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Chapter  I 
INTRODUCTION 

Two  very  fundamental  functional  properties  of  an  associating 
system  of  ma c romol ecnl es  are:(l)  the  s t o i chiome tr i e s  of  the 
interactions  that  occur  (i.e.  the  compositions  of  the  molec- 
ular species  that  exist)  in  solution  and  (2)  the  equilib- 
rium constants  that  govern  these  interactions. 

To  determine  these  properties  of  a  system,  thermodynamic 
equilibrium  techniques  such  as  osmometry,  light  scattering, 
equilibrium  sedimentation  and  small  angle  X-ray  scattering 
are  usually  employed.    The  variations  of  molecular  weight 


The  Law  of  Mass  Action  states  that  the  rate  of  a  reaction 
at  a  given  time  is  proportional  to  the  active  masses  of 
reacting  substances  present  at  the  time.  Thus,  for  the 
reversible   reaction: 


A  +  B  <=>  AB  , 


the  speed  with  which  A  and  B  react  is  proportional  to  the 
product  of  their  activities.  The  rate  of  the  forward  reac- 
tion, v^_  is  given  by  ki(A)(B),  and  that  of  the  reverse 
reaction,  v,  ,  by  k,(AB).  At  equilibrium.  Vi  =  V2  and 
hence,  also  at  equilibrium. 


ti(A)(B)  =  k,(AB) , 


and  therefore, 

l^i/k:,  =   (AB)/(A)(B) 

-  1  - 


averages  with  the  concentrations  of  the  interacting  constit- 
uents,  obtained  by   these   methods,   provide   diagnostic  data 

which,   ideally,   would  fit   only  one   correct  reaction  model 

2 
for  the  system.   Frequently,  however,  more  than  one  reaction 

scheme   can   fit   all   of   th«   molecular   weight   data   equally 

well,   and  additional  experiments  are  required  to  eliminate 

the  remaining  incorrect  reaction  models. 

One  potential  source  of  complementary  distinguishing  data 
is  the  behavior  of  the  system  during  a  mass  migration  exper- 
iment such  as  velocity  sedimentation.  If  the  shape  of  a 
system's  migrating  boundary  is  a  distinctive  indicator  of 
its  reaction  pattern,  as  it  is  generally  thought  to  be,  then 
a  comparison  of  the  real  system's  boundary  shapes  with 
boundary  shapes  predicted  for  different  possible  model  sys- 
tems should  help  to  eliminate  at  least  some  of  the  ambigui- 
ty. 

To  predict  the  shapes  of  migrating  boundaries  for  inter- 
acting systems  during  velocity  sedimentation  experiments  it 
is  necessary  to  invoke  numerical  methods,  since  analytical 
solutions  do  not  exist  for  the  differential  (continuity) 
equations  that  describe  the  behavior  of  interacting  systems 
during  transport.  Several  such  numerical  methods  have  been 
developed  and  computerized  for  simulating  velocity  sedimen- 
tationintheultracentrifuge  (2-38). 


This  last  quantity  defines  the  equilibrium  constant,  K 
for  the  reaction  with  the  s t o i chiome t ry  given  above. 

An  overview  of  equilibrium  methods  can  be  found  in  (1). 


If  these  techniques  are  to  be  useful  in  distinguishing 
between  various  possible  models  for  an  interacting  system 
then  the  sediment ing  boundaries  of  significantly  different 
model  systems  must  not  be  alike.  We  are  thus  interested  in 
knowing  the  extent  to  which  boundaries  of  different  systems 
may  be  distinguished  from  one  another. 

In  the  work  described  herein  the  distorted-grid  method  of 
Cox  (29-35)  was  used  to  simulate  velocity  sedimentation  of 
various  systems  undergoing  mixed  association  reactions  (i.e. 
reactions  between  dissimilar  protein  molecules)  and  an  at- 
tempt was  made  to  assess  the  diversity  in  the  resulting  sim- 
ulated gradient  profiles.   Only  AB   svstems  were  considered, 

n   ■' 

and,  within  this  class  of  systems,  the  focus  was  almost  ex- 
clusively on  AB^  systems.   ABj  systems  are,  conceptually  and 

computationally,   the   simplest   systems   among   AB    svstems 

n    ■' 

(which  are   the   simplest   among  A  b    systems)   and  therefore 

n  m    ^ 

represent  a  convenient  and  logical  starting  point  for  such 
studies.  These  studies  are  described  and  discussed  in  chap- 
ters three  and  four. 

The  theory  of  the  Cox  distorted-grid  model  for  simulating 
velocity  sedimentation  is  presented  in  chapter  two.  This 
was  written  primarily  as  an  exercise  and  is  presented  here 
for  the  sake  of  completeness.  More  elegant  treatments  of 
this  theory  may  be  found  elsewhere  (29-35)  . 


Chapter  II 
THE  DISTORTED  GRID  METHOD 

In  the  distorted  grid  model  the  ultracentrifuge  cell  is 

treated  as  an  array  of  narrow  concentric  cylindrical  sectors 

(boxes).   The  distribution  of  solute   in  the  model  cell   is 

described  by  three  arrays:  r,  r,  and  c.   Arrays  r  and  r  c on— 

tain  the  values  of  the  distances  from  the  axis  of  rotation 

to  the  box  b  oundar  ies  and  box  centers,  respectively,  and  the 

array  c  contains  the  values  of  solute  concentrations  at  the 

box  centers.   Thus  the  j th  box  in  the  model  would  have  upper 

(nearer  to  the  axis)  and  lower  boundaries  at  r.   and  r.,-, 

respectively  and  a  solute  concentration,  c.,  at  its  center, 

r  . 
J 

When  two  different  constituents  (ie.  types  of  monomeric 
species)  are  present,  the  arrays  r,  r  and  c  become  two  di- 
mensional so  that  the  distribution  of  each  constituent  is 
described  independently  by  its  corresponding  track  of  r,  r 
and  c  values.  For  example  the  j th  boundary  position,  the 
box  center  position,   and  the  box  center   concentration  of 

constituent  A  would  be  given  by  r    .^  ^a      and  c* 

A,  J  '   A, J  '        A, J  * 

Velocity  sedimentation  is  simulated  in  the  model  cell  by 
making  the  appropriate  changes  in  the  three  simulation  ar- 
rays (r,  r  and  c)  to  describe  sedimentation  and  diffusion  in 
short  time  intervals. 
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An   ultracentrifuge   cell   at   time   zero   in   a   velocity 

sedimentation  experiment   is  represented   in  figure  1.    The 

length  of  the  cell  is  r,  -  r^  where  r,   and  r.   are  the  dis- 

D      t  b        t 

tances  from  the  axis  of  rotation  to  the  bottom  and  top  of 

the  cell.   The  initial   sharp  solute  boundary  (meniscus)   is 

at  r   and  the  initial  plateau  region  is  between  r   and  r.  . 
o  IT  b  o       b 

Typical  values  for  r.,  r,  and  r   are  6.0,  7.0  and  6.1  cm. 
t    b       0 


Figure  1:   Ultracentrifuge  Cell  at  t 


2  . 1    INITIALIZATION  OjP  THE  SIMULATION  ARRAYS 

To  simulate  a  velocity  sedimentation  experiment^  using 
the  distorted  grid  model,  the  simulation  arrays  (r,  r  and  c) 
must  first  be  initialized  so  that  they  describe  the  initial 
conditions  in  the  cell.  Array  r  (box  boundary  positions)  is 
initialized  so  that  the  boundaries  are  evenly  spaced^  or 
otherwise  systematically  placed,  between  positions  about  0.1 


cm  above  r^  and  the  cell  bottom^. 

To  avoid  ambiguity  in  assigning  boj:  center  concentra- 
tions, one  of  the  box  boundaries  is  positioned  at  the  menis- 
cus, r  jhe  initial  position  of  this  boundary  is  then  used 
as  a  reference  point  for  assigning  positions  to  the  bounda- 
ries above  and  below  r^.  por  instance  if  boundaries  sepa- 
rated by  a  constant  distance,  Ar,  are  desired  and  the  nth 
boundary  is  given  the  position  of  r^  ( r  ^^  =  r  )  then  the  po- 
sition of  the  boundary  immediately  below  r  is  r  +Ar  and  the 
position  of  the  boundary  immediately  above  r   is   r  -Ar, 

The  following  algorithm  illustrates  equal  box  width  ini- 
tialization of  boundaries  for  a  100  box  model  of  a  1.0  cm. 
cell.  Note  that  the  number  of  box  boundaries  must  always  be 
one  greater  than  the  number  of  boxes  and  that  giving  the 
eleventh  boundary  the  position  r  places  the  first  ten  boxes 
above  r 


The  upper  and  lower  limits  of  box  boundary  positions  are 
not  necessarily  the  same  as  the  positions  of  the  top  and 
bottom  of  the  real  cell. 


Al  gor  i  tlim 

I)  number  of  boxes  =  100 

II)  cell  length  =  1.0 

III)  Ar  =  (.cell  length)  /(num.  boxes) 


IV) 

V) 

VI) 


r(ll)  =  r. 


J  =  11  ■    .. 

Repeat  from  i  =  1  to  10 

A)  j  =  j  -  1 

B)  r(j)  =  r(j+l)  -  Ar 
End  Repeat 

VII)   Repeat  from  j  =  12  to  (num.  boxes 
A)  r(j)  =  r(j-l)  +  Ar 
End  Repe  a  t 
End  Al gor  i  thm 


+  1) 


In  some  simulation  programs  the  bo unda ries  are  not  evenly 
spaced  in  the  initialization  routine.  Instead  they  are  po- 
sitioned so  that  the  distance  to  each  boundary  from  the  axis 
of  rotation  is  greater  than  that  of  the  boundary  immediately 
above  it  by  a  constant  factor  y. 


r«  +  Ar 


(1) 


where  Ar  is  the  average  box  width.  Since  r  is  always  much 
larger  than  Ar ,  y  is  just  slightly  greater  than  one.  The 
following  algorithm  illustrates  initialization  of  box  bound- 
aries using  a  factor  of  y  for  spacing. 


Al gor  i  thm  ^ 

I)  number  of  boxes  =100 

II)  length  of  cell  =  1.0  '  ■'   ''■ 

III)  Ar  =  (length  of  cell)  /(num.  boxes) 

IV)  y  =  (r^  +  Ar)/r^ 

V)  r(ll)  =  r 

o 

VI)  j  =  11 

VII)  Repeat  from  i  =  1  to  10 

A)  j  =  j  -  1 

B)  r(j)  =  r(j+l)/Y 
End  Repeat. 

VIII)  Repeat  from  j  =  12  to  (num.  boxes  +  1) 
A)  r(j)  =  r(j-l)Y 

End  Repeat. 
End  Al gor  i  thm 


The  number  of  boxes  used  in  a  simulation  is  not  always  100 
but  it  usually  is.  Fewer  boxes  may  be  used  to  conserve  com- 
puter time,  but,  since  a  decrease  in  the  number  of  boxes 
will  result  in  an  increase  in  the  box  width,  the  interpola- 
tions of  concentrations  between  box  centers  that  occur  dur- 


ing  the  simulation  will  become  less  accurate  and  the  simu- 
lated result  less  reliable.  On  the  other  hand,  very  little 
is  usually  gained,  in  terms  of  precision,  by  using  more  than 
100  boxes . 

With  box  boundary  positions  initialized  as  described 
above  the  box  center  concentrations  are  initialized  to  match 
the  situation  in  the  real  cell.  Boxes  below  r  are  given 
concentrations  equal  to  the  initial  plateau  concentration, 
C  ,  and  boxes  above  r  are  assigned  initial  concentrations 
of  zero.  The  positions  of  the  box  centers  are  computed  from 
the  boundary  positions 


r  .  = 
J 


^j  ^  ^j-i 


(2) 


The  algorithm  that  follows  illustrates  the  initialization  of 
arrays  c  and  r . 
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Al gor  i  thm 

I)  Repeat  from  j  =  1  to  10 

A)  c  .    =  0.0 

B)  i.  =  (,.  ^  ,,^^)/2   ^    .^.-  , 
End  Repeat. 

II)  Repeat  from  j  =  11  to  (num.  boxes) 

A)  H.  =  c„ 

B)  ij  -  (Xj  +  r  j^i)/2 
End  Repe  a  t . 

End  Al gor  i  thm 


In  a  two  constituent  model  both  constituent  tracks  are 
initialized  the  same  way  that  a  single  constituent  array  is 
initialized.  Box  boundary  positions  are  given  the  same  val- 
ues in  both  tracks  < ^^  =  r^  - )  .  The  only  difference  ini- 
tially between  constituent  tracks  is  in  their  plateau  con- 
centrations C^^^  and  C^Q, 

The  algorithm  above  creates  a  sharp  initial  boundary, 
which  is  the  usual  choice.  It  is  also  possible,  and  some- 
times desirable,  to  begin  with  an  initial  situation  in  the 
model  cell  describing  a  diffused  initial  boundary.  In  this 
case  the  initial  concentration  in  each  box  would  either  be 
read  directly  by  the  program,  or  arrived  at  by  subjecting  a 
sharp  boundary  to  a  few  rounds  of  simulated  diffusion  wit fa- 
out  the  corresponding  rounds  of  simulated  sedimentation. 
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Once  the  arrays  have  been  initialized,  the  simulation  be- 
gins,  with  alternating  rounds   of   sedimentation  and  diffu- 


2 . 2    SIMULATION  OF  DIFFUSION 

2.2.1    Sj s t em^  not  Un d e r^ oi n^  Mixed-Association 

A  short  interval,  Atp,  of  diffusion  (without  sedimenta- 
tion) is  simulated  by  computing,  for  each  constituent,  the 
mass  of  that  constituent  that  would  flow  upward  across  each 
boundary  in  its  track  during  At_  seconds  and  subsequently 
computing  new  constituent  concentrations  for  each  box. 

In  a  single  constituent  system  the  rate  at  which  constit- 
uent mass  moves  upward  in  the  cell  per  unit  cross  sectional 
area  at  boundary  j  is  given  by  Ficks  law  as: 


r  3c  1 
-J  =  D   (3) 

L  «'  Jj 


where  D   and  Oc/3r).  are  the  local  average  diffusion  coef- 
ficient and  concentration  gradient  at  r.. 

The  total  flow  rate  of  mass  upward  at  boundary  j,  F.,  is 

-J  times  the  cross  sectional  area  of  the  cell  at  r. 

J  ■ 


r  3c  1 

L  ''    Jj 


''j  =  "^-^jBj  '"-  (4) 


where  ber^  is  the  cros 
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sectional  area  of  the  cell  at  r.    b 


and  0  are  the  height  and  sector  angle  of  the  cell. 

The  mass  that  would  flow  past  boundary  j  during  the  time 
from  t^  to  t-^  ,    m-,     is  given  by  the  following  expression. 


"j  =  "«- 


'l>' 


r  ac   1 


arr  Be  1     11 


D.       dt 

J 


L    ^"^    Jo,j        3HL   ^'^    Jj      JJ 


=    ber. 


r     3o     1 


r   ac   1     1 


L-* 


I       0  c       I  pO  I       DC       I  I 

D^     . dt    +       — n  . dt 

L   ''   Jo.j  Jat   \   3-   Jj    J 


(5) 


The  second  integral  in  equation  5  cannot  be  evaluated  di- 
rectly.  However,  since 


r   ac  1 


r   ac  1 


)     ---       >>  — -D   — -  , 

L  ''   Jo.j     3t   L  3'  Jj 


(6) 


the  contribution  of  the  second  integral  can  be  considered 
negligible  as  long  as  the  interval  from  t  to  t,  is  very 
shor  t . 


1  imit       in   =  bGr. 
(At  ->  0)   J        J 


r  30  1 


0.  J 


dt 


3r 


L  "  Jo.j 
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(7) 


Since  D^  ;  and  (3c/3r)^  .,   the  average  diffusion  coeffi- 


cient and  local  gradient  at  t 
the  last  expression  gives: 


ire  constant,   evaluation  of 


where 


r  3c  1 

L  "^    Jo.j 


^'d  =   ti 


(8) 


The  mass  of  solute  (constituent)  that  a oc nmul a t e s  in  box 
j  during  Atp,  Am,,  is  equal  to  the  mass  that  flows  upward 
across  boundary  j+1,  into  the  box,  minus  the  mass  that  flows 
upward,  across  boundary  j,  out  of  the  box. 


-  -"j+l  -  -nj 

r  r  3o  1  r  r  ac  i      i      i 

^  ^^    '^j+lPo,j  +  l  -:--  -   "^jD^.j Atp 

L         L  «'  Jo.j.i   L    '\  3r  j„,jj   J  ^^^ 


Dividing  the  change  in  mass.  Am 


j,  by  the  volume  of  box  j. 


i-vX 


Vj  =  berj(rj,i  -  rj), 
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(10) 


results   in  the  following  expression  for   the   change   in  the 
concentration  of  box  j  during  At,,, 

,-  •  "  I .  ■<•  •  . 

Am  . 
AC.  = i-      -  .   ,.'  ,   -  '   '■  - 


r  3o  1 


r  3c  1      1 


"j+iOj+i 


L  ''    Jj  +  i   ''  \   '^    Jo.jJ 


At, 


'j'"^j+l  -  "^j' 


(11) 


Finally,   the  new  concentration  in  box  j. 


;,  at  the 


new , J ' 
end  of  diffusion  is  equal  to  the  original  concentration  in 

box  j  plus  the  change  in  its  concentration. 


new,  J 


=  o,j  ^  ^Cj 


(12) 


Thus,  to  simulate  Atp  seconds  of  diffusion  (in  a  one  con- 
stituent system),  equations  (11)  and  (12)  are  applied  to 
each  box  in  the  array. 

Expressions  analogous  to  (11)  and  (12)  for  the  diffusion 
of  constituents  A  and  B  in  a  noni nt er a c t i ng  two  constituent 
system  are  the  following. 
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A,  j  +  l"A.  j  +  1 


Ac 


A,  j 


L  "  JA,j+l 


^a,j5a,j  ----     Ati, 

L  ^^  JA.jJ 


'A,j<'^A,j  +  l 


^A,j' 


(13a) 


B, j+l"B, j+1 


Ac 


B,  j 


.__J:____J?iJil 

'B,j('^B,j  +  l  -  -^B.j' 


"B,j"B,j 


Tacgjl   1 


At, 


(13b) 


new.A.j  =  "Ao.j  +  '^^K.  j 


(14a) 


"new.B.j  ~  °Bo,j  ■*■  ^*^B,j 


(14b) 


To  simulate  diffusion  in  a  two  constituent  system  (not 
undergoing  mixed  association  reactions)  equations  (13a)  and 
(14a)  are  applied  to  each  boi  in  the  constituent  A  track  and 
equations  (13b)  and  (14b)  are  applied  to  each  boi  in  the 
constituent  B  track. 


2.2.2    Ml  led  A^i  o^i  a^l  n^  Ji^^^m^ 

A  mixed  association  reaction  is  one  in  which  constituent 
monomers  of  different  types  associate  to  produce  a  heteroge- 
neous (mixed)  aggregate. 


mA  +  nB  <=>  A 


B„ 
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A  system  undergoing  mixed  association  may  contain  a  num- 
ber of  different  mixed  aggregates  in  addition  to  constituent 

monomers   and  possibly   some   self   association  products   (A 

n' 

Bjjj)  .  The  distribution  of  constituents  among  the  various 
species  present  is  governed  by  the  set  of  association  con- 
stants for  all  aggregate  species  present. 


m    C. 


A  B 
m^n 


"A.B. 


pm  pu 


(15) 


and,  for  a  given  set  of  association  constants,  will  depend 
on  the  total  concentrations  of  each  of  the  two  constituents. 

The  flow  of  a  constituent  in  response  to  its  gradient  in 
such  a  system  is  the  sum  of  constituent  weight  fractions  of 
the  flows  of  the  individual  species  in  response  to  their  in- 
dividual species  gradients.  The  flow  of  constituent  A,  for 
example,  is  the  sum  of  constituent  A  weight  fractions  of  the 
mass  flows  of  all  species  containing  constituent  A. 

The  constituent  A  gradient  is  the  sum  of  constituent  A 
weight  fractions  of  the  individual  species  gradients.  Since 
the  individual  species  have  different  diffusion  coeffi- 
cients, the  flow  of  A  depends,  not  only  on  the  magnitude  of 
the  A  gradient,  but  on  its  composition  as  well.  Whereas  the 
distribution  of  constituent  A  among  species  at  any  point  in 
the  cell  depends  only  on  the  two  constituent  concentrations. 


17 
the  distribution  of  the  A  gradient  among  species  gradients 
at  that  point,  which  will  determine  the  average  rate  of  dif- 
fusion of  constituent  A  at  that  point,  depends  on  both  of 
the  constituent  gradients  in  addition  to  the  constituent 
concentrations.  Since  the  gradient  of  constituent  B  affects 
the  distribution  of  the  A  gradient  among  gradients  in  molec- 
ular species  that  have  different  diffusion  coefficients,  it 
must  also  affect  the  flow  of  constituent  A.  (i.e.  a  cross 
diffusion  effect)  The  total  flow  of  constituent  A  can  be  ex- 
pressed as  the  sum  of  a  self  diffusion  term  (diffusion  of 
constituent  A  in  response  to  the  A  gradient)  and  a  cross 
diffusion  term  (diffusion  of  A  in  response  to  the  B  gradi- 
ent )  . 


■^A  =  "aa  ----  +  Dab  -:-- 
L  ar  J       L  3r  J 


(16) 


^AA  ^^^    ^AB  ^^  this  expression  are  the  direct  and  cross  dif- 
fusion coefficients  of  constituent  A. 

Similarly,  the  flow  of  constituent  B  would  be  affected  by 
the  A  gradient,  and  its  flow  is  likewise  expressed  as  the 
sum  of  a  direct  and  cross  diffusion  term. 


-J„ 


^BB 


L  ^"^  J 


'"scat"' 


L  ^"  J 


(17) 
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Expressions  analogous  to  equations  (11)  and  (12)  for  the 
change  in  box  constituent  concentrations  and  resulting  new 
box  concentrations  after  a  time.  At  of  diffusion  in  a  two 
constituent  (A  and  B)  system  undergoing  mixed  association 
are  obtained  by  replacing  each  of  the  diffusion  terms, 

r  3c  1 

D. ,  ■         ' 

.        'l   ^^    Jj 

with    tlie    appropriate    cross    and    direct    diffusion    terms. 


'^"A.j    =    ^H    ^ 


r^^AT^ 


rac3,i      n 


A,j+1    "AA,j+l       , 

L  L  >    ^'^ 


AB, j+1 

LlLhti L!!JiiiJJ 


'A,j('^A,j  +  l 


^A,j' 


^A,j    "AA.j 


L  L  L   ""^   J 


3r 


D 


Tacgjl   n 


AB,  j 


A,j<'"A,j  +  l    -    "^A,:' 


(18a) 


'A,new,j    -    <=A,j    +    AC^^j 


(19a) 
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^°B,j    =    AtD 


B,j+1    "BB,j+l 


rar_i     n 


D 


AT 


J=_:!_Jiil. 


BA, j+1 


i-!f_JiilJJ 


"b,  j<'^B,  j  +  1    -    "^B,  j' 


r       r        rac3,i 

^B.j    °BB,j    -:--       ^    "BA.j 

L_ L L!!_Ji k_!!JiJJ 

^B.jf'^A.j  +  l    -    -^B.j) 


3r 


(18b) 


=  B,iie*,j    -    =B.j    ^    '^Cgj 


(19b) 


2.2.3         SoEZSliSA    Gradients    a_t    Box    B o un d a rj, ei 

The  boundary  gradients,  (dc/dr),  or  simply  g,  used  in  the 
diffusion  routine  (equations  13  and  14  and  equations  18  and 
19)  are  computed  immediately  prior  to  each  round  of  simulat- 
eddiffusion. 

The     gradients,      g  at     the     box     boundaries,      r.,      in     a     one 

constituent  systems  are  computed  in  two  steps.  In  the  first 
step,  gradients,  s';,  at  positions  midway  between  adjacent 
box  centers,  are  computed  along  with  the  positions,  r'  .  ^o 
which    they    correspond. 


g  ' 


T  .     _ 


(20) 


■^j-l 
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r  .  -  r  .  , 

r'   =  -i i=i 

J        2 


(21) 


Because  the  box  widths  are  generally  not  uniform  (sedi- 
mentation operations  distort  the  spacing  of  the  grid)  the 
position,  r'  .,  of  the  midpoint  between 
necessarily  the  same  as  the  position, 
between  the  boxes  (figure  2). 


r .  ^  and  r.   is  not 
r  ,   of   the  boundary 


'  ^^      » 


Figure  2;   Gradient  Interpolation  Scheme 


The  gradients,  g.,  at  the  box  boundary  positions,  r.,  are 
J  J 

computed  by  linear  interpolation  between  the  gradients,   g', 

at  the  box  center  midpoints.   For  instance,   the  gradient  g. 

at  r   is  interpolated  between  g'.  at  r'.  and  g'.    at  r'. 

J  J        J         J  "*"1        j  '*"1 


^j  = 


'«'j+l  -  8'j'<>^j  -  ''j' 
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(22) 


Analogous  expressions  for  computing  the  gradients  of  each 
constituent,  g^  and  gp  at  its  own  boundary  positions  in  a 
two  constituent  model  are  given  below. 


g  ' 


A,  j 


'a.j  -  "^A.j-l 


(23a) 


■^A.j  -  'A,j-1 


A.  j 


(24a) 


S  ' 


'B,j 


'B, j-1 


B,  j 


B,j  -  "^B.j-l 


(23b) 


■^B,  j  ~  "^B,  j-1 


B,  j 


(24b) 


'a,j  -  8'a, 


*»'a,j  +  1  -  8'A,j>(>^A,j 


A,j' 


A,  j+1 


A,  j 


(25a) 


(g' 


'b.j  -  8  B,j 


B,  j  +  1 


'b,  j)("^B,  j 


B,j' 


B,  j  +  1 


B,  j 


(25b) 


The  cross  gradients  required  for  diffusion  in  mixed  asso- 
ciating systems  (ie.  the  gradients  in  constituent  B  at  the 
constituent  A  boundaries  and  the  gradients  in  constituent  A 
at  the  constituent  B  boundaries)  are  computed  in  a  similar 


vi- 
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fashion.  The  first  step  is  again  to  compute  gradients  mid- 
way between  box  centers.  For  example,  in  determining  gradi- 
ents in  constituent  B  at  A  boundaries,  the  first  step  is  to 
calculate  S'g  j  and  r'fl.j  between  each  adjacent  pair  of  B 
box  centers  (equations  23b  and  24b).  A  complication  arises 
at  this  point  because  box  boundaries  in  the  two  constituent 
tracks  are  usually  not  in  register.  That  is,  the  position 
of  boundary  j  in  constituent  track  A  is  not  necessarily  the 
same  as  the  position  of  the  j  th  boundary  in  constituent 
track  B.  In  fact  r^  and  ig  .  may  be  separated  by  several 
boxes.  This  lack  of  correspondence  in  position  between  the 
boxes  of  the  two  tracks  develops  durring  the  sedimentation 
routine,  in  which  the  boundaries  of  the  two  tracks  are 
shifted  independently. 

The  second  step  in  finding  the  B  gradient  at  an  A  bounda- 


ry. 


'A,j 


p       is   to   locate   the   pair   of   gradients,   g' 


B,  n 


and 


8  g  n-1 *   located   at   positions,   r'p     and   r' 


B,n-1, 


bracket  r 


^  1*        This   is   accomplished  by   searching  from   the 
top  of   the   array   through   successive  values   of   r'    .   until 


the  first  occurrence  of  r' 


g  j  greater  than  r^  ■   is  located. 


The  gradient  in  B  at  r^^  ^  is  then  interpolated  between  the 
gradient,   s'g  ^,   at  this  position,  which  is  r'g  ^,     and  the 

gradient,  g'       at  the  position  r'      -, 
D,  n  -^  B,  n— 1  • 


(g' 


^BA,j  -  8  B,n-1 


B,n  "  8  'B,n-l'  <  "^A, 


B,n-1' 
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(26) 


B.n 


B,  n-1 


2.2.4    Com£ ul a_tl on  ot,    Dlf  f  ns  ion  Co eHi o J, e n^^ 

The  difficulty  involved  in  computing  the  appropriate  diffu- 
sion coefficients  for  simulated  diffusion  operations  depends 
on  the  compleiity  of  the  system  being  studied. 

2.2.4.1    Non-Interacting  Systeas 

The  simplest  systems  are  those  which  undergo  no  associa- 
tion whatsoever.  Since  the  constituent  molecules  in  systems 
of  this  type  spend  all  of  their  time  as  monomers,  the  con- 
stituent diffusion  coefficients  are  the  same  as  the  monomer 
diffusion  coefficients  and  are  concentration  independent. 
For  non- as  so c i a t i ng  systems,  then,  each  constituent  diffu- 
sion coefficient  is  computed  only  once. 


RT  f 

D  =  -- ,  where  f  =  f  __ 
Nf  "f 


(27) 


f/fg  is  the  ratio  of  the  frictional   coefficient  of  the 
molecule,   f,   to  that  of  an  unhydrated  sphere  of  equal  vol- 

"°"'  '^o-   fo  's  given  by  Stokes  law. 


fj,  =  SnnRo'  where  R^ 


I  4nN  I 


24 


(28) 


and  11  is  the  solvent  viscosity.  M  and  v  are  the  molecular 
weight  and  partial  specific  volume  of  the  ma cromol ecul ar  so- 
lute . 

Alternatively,  D  may  be  obtained  from  the  sedimentation 
coefficient  at  infinite  dilution  and  the  molecular  weight, 
eliminating  the  need  to  know  f. 


S   = 


M(l  -  vp) 


Nf 


(29) 


D  = 


RT  ^    S^RT 

Nf    M(l  -  vp) 


(30) 


2.2.4.2    Self-Associating  Systens 

In  self  associating  systems  a  constituent  molecule  spends 
part  of  its  time  as  a  monomer  and  part  of  its  time  as  part 
of  any  number  of  aggregates,  diffusing,  at  any  instant,  at  a 
rate  determined  by  the  diffusion  coefficient  and  local  gra- 
dient of  whichever  species  it  is  existing  in  at  that  in- 
stant. Since  the  diffusion  coefficients  of  the  monomer  and 
aggregate   species   are  not   the   same,   the  diffusion   coeffi- 
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cient  of  a  constituent  molecule  that  spends  part  of  its  time 
as  one  or  part  of  each  different  species  must  be  an  average 
of  some  sort  over  all  species  present. 

The  total  mass  flux  of  a  self  associating  system  is  the 
sum  of  monomer  and  aggregate  species  fluxes  in  response  to 
their  respective  species  gradients. 


ac 


di 


(1  <  j  <  n)  (31) 


Thus  the  contribution,  to  the  total  flux,  of  each  species 
is  proportional  to  its  species  gradient  and  the  average  dif- 
fusion coefficient,  D,  is  a  species  gradient  weighted  aver- 
age . 


rac^l 

D 

L  ^'J 


(32) 


rad 

^ar  J 


(1  <  j  <  n)  (33) 
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The   species   gradients   Oc^/ar),   in   equation  33   can  be 
re-eipressed  in  terms  of  the  monomer  concentration  C  ,  and 


the  species 


sociation  constants,  K.,  as  follows. 


K.  =  -i 


<34) 


Cj  =  Kj  Ci 


(35) 


'(c)    a(K^cJ) 

)L_  = J.___ 

3t        ar 


(36) 


If   the   association   constants   are   pre s s ur e- independent 
then  equation  36  may  be  rewritten  as: 


3(c)  facJl 

i-    =    K.     — i 

ar  "'l    ^"^J 


(37) 


8^«j)  _      racjirac^i 


ar 


L^'^iJL   ^n 


(38) 


Partially  differentiating  with  respect  to  C  and  substi- 
tuting the  result  for  the  species  gradients  in  equation  33 
gives 


D   = 


L   ^"^J 


^JKjCf 


-1 1 

L   ^'J 
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(39) 


which    is    equivalent     to 


h4 


(40) 


2.2,4.3    Mixed-Associating  Systems   ' 

The  derivation  of  the  following  expressions  for  the  average 
diffusion  coefficients  in  a  mixed  associating  system,  in 
terms  of  the  monomer  concentrations,  C  and  Cd  ,  though 
St rai ght- forward,  is  quite  cumbersome  and  will  not  be  pre- 
sented here.  A  detailed  derivation  of  these  expressions  may 
be  found  elsewhere,   (35) 


"aA  -  '"^AA^BB  "  i^ABlBA)  ^  ^l 
"ab  =  '"^ABlAA  -  '*AA«AB)  /  ^^ 


(41) 


BB 


BA 


'^BB^AA  ~  "^BAlAB)  '  ^1 
'"^BA^BB  "  "^BB^BA)  /  ^4 
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wher  e 


^'^         "^BBlAA  -  lABlBA 

''aA  =  ^;^iDijfA,ijKijCr'c^ 

(0  <  i  i  m)  (1  i  j  <  n) 

"aB  =  &'>ijfA,ijKijCic^-l 

(0  i  i  i  m)  (1  <  j  <  n) 

%B  =  ^^JDijfA.ij^ijCicJ-' 

(0  <.  i  <  m)  (1  i  j  <.  n) 

Sa  =  ^^iDijfA.ij^ijCi-^CJ 

(1  <  i  i  m)  (0  <  j  <  n) 


(42) 


and 


'aA  =  ^^fA.ijiKijCr'^i- 

(0  <  i  <  m,  1  <  j  <  n) 

'aB  =  ^^fA,ijJKijCicJ-l. 

(1  <  i  1  m,  0  i  j  i  n) 
*BB  =  ^^fB,ijJKijCicJ-l, 

(1  <  i  1  m,  0  <.  j  <.  n) 
*BA  =  ^^fB.ijiKijCi-^cJ, 

(0  <  i  <  m,  1  <  j  <  n) 


(43) 


The  summations  in  these  expressions  are  over  all  species 
"^Bj  that  are  present  in  the  system  with  i  and  j  equal  to  or 
greater  than  the  values  given  below  each  expression.  The 
summation  limits,  m  and  n,  are  the  largest  numbers  of  A  and 
B  monomers   occuring  in  any   of   the  aggregates.    f. 


A,  ij 


and 
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fg  J-  are  the  weight  fractions,  in  species  A.fi-  of  constit- 
nent  A  and  constitnent  B,  respectively,  D..  is  the  diffusion 
coefficient  of  species  A.b^  and  Kjj  is  the  association  con- 
stant for  A .B . . 


^A,  ij  = •   fB,ij  = (44) 

'"a  +  JMb  i"A  "■    JMb 


^ii     =    — T-T  (45) 

a'-b 


2.3    SIMULATION  OF  SEDIMENTATION  '   - 

Sedimentation  is  simulated  by  shifting  box  boundaries 
downward  in  the  array  and  away  from  the  axis  of  rotation. 
In  the  time  interval  for  one  simulated  sedimentation  shift, 
Atg,  which  corresponds  to  At„  seconds  of  real  experimental 
time,  a  boundary  initially  at  r  will  move  to  a  new  posi- 
tion, r  so  that 
new, J '  *"  tnai 


^ew,j  =  rjefSj"'^*s'  (46) 


wher  e 
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u  is  the  rotor  speed  in  radians  per  second  and  B.  is  the 

local  average  sedimentation  coefficient  at  r 

J  • 

Solute  does  not  cross  boundaries  in  the  model  during  a 
sedimentation  shift;  all  solute  initially  contained  between 
boundaries  r^^^  and  r.  will  be  contained  between  boundaries 

^new,j  +  l  ^""^  "^new.j  **  **"  ^""^  °^  ""^  sedimentation  shift. 
Because  the  boundaries  move  at  different  rates,  however,  the 
box  volumes  and  box  center  concentrations  will  be  changed'*. 
The  initial  volume  of  box  j,  Vj ,  and  the  volume  of  box  j  at 
the  end  of  a  sedimentation  shift,  V^^^^^,  are  given  by  the 
following  expressions. 


"^'i^i 


(47) 


be(r2 

new, j+1 


new, J ' 


(48) 


Thus  if  the  initial 
center  concentration)  then  the  concentration  at  the  end  of 


oncentration  in  box  j  was  c.   (box 


the  sedimentation  shift. 


new, J ' 


would  be  as  follows 


Because  the  cell  is  sector  shaped,  the  movement  of  two  ad- 
jacent boundaries  away  from  the  axis  of  rotation  will  re- 
sult in  a  change  in  the  volume  and  concentration  of  the 
intervening  box  even  when  the  boundaries  move  at  the  same 
rate. 


"iCr 


J+1 


) 


new,  J 


new, j+1 


.11. 

2 
"^new,  j 
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(49) 


For  a  two  constituent  system  the  equations  for  sedimenta- 
tion, analogous  to  equations  46  and  49,  are  written  as  fol- 
lows. 


'^new,A,j  =  '^A,j"P(SA,j"^^ts) 


(50a) 


new , A, j 


'a,  j<'^A,j  +  l 


^A,jJ 


new, A, j+1     new,A,j 


(51a) 


"^new.B,  j  =  "^B,  j=^p(Sb_  jO,  At^) 


(50b) 


new , B, j 


=  B,  j('^B,j  +  l 


^B,j' 


new.B,  j  +  1    "^new.B,  j 


(51b) 


where  '^^  j   and   Sg  j   are   the   local   average   sedimentation 
coefficients  of  constituents  A  and  B. 


2.3.1    Com£ut a tlon  of  S^dimeni a_tlon  Co£ f f _i olen_t^ 
2.3.1.1    Non-Associating  Systems 

The  local  average  sedimentation  coefficients  are  general- 
ly concentration  dependent.   In  non-associating  systems  only 
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the  hydrodynam  i  c  dependence  of  'S  needs  to  be  considered.  In 
this  case,  'S  is  obtained  directly  from  S  ,  the  sedimenta- 
tion coefficient  at  infinite  dilution. 


5.  = 
J 


1  +  kC. 


(52) 


where  C.  is  the  concentration  at  boundary  j  and  k  is  the 
hydrodynamic  constant.  For  a  two  constituent  non- as  so c i a t- 
ing  system 


A,  o 


A,  j 


1  +  ^aC 


aa'^at,  j  ■^  ^ab'^bt,  j 


(53a) 


B,  0 


B,  j 


1  +  K^C, 


bb'^bt,  j  "^  ''ba^at,  j 


(53b) 


where  C^^  and  Cd't-  ;  are  the  concentrations  of  constituents 
A  and  B  at  boundary  j  and  k^^,  t^^^  t^g,  and  t^^  are  the  ap- 
propriate hydrodynamic  constants. 

2.3.1.2    Self-Associating  Systens 

If  constituent  A  self  associates  then  its  average  sedi- 
mentation  coefficient   will   depend   upon   its   distribution 


;..-.  ..-.1*^ 
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among  aggregate  species  and  monomer.  The  average  sedimenta- 
tion coefficient  of  constituent  A  is  a  weight  average  over 
all  species  in  which  constituent  A  participates.  An  appro- 
priate expression  for  the  average  sedimentation  coefficients 
of  a  self  associating  constituent  A  at  indefinite  dilution 
(not  corrected  for  hydrodynamic  dependence),  S  _  js  the 
f ol 1 ow  ing . 


A,  0 


(54) 


where   S,      jj   ^jj^   ideal   sedimentation   coefficient   of 
Aj.o 

species  A^  (ie.,  in  the  absence  of  hydrodynamic  effects)  and 

c  •  *  ' 

^^   is  Its  concentration. 

i 
Since 


Ka.Ca 


(55) 


equation  54  can  be  reexpressed  in  terms  of  the  monom 
centration,  C^,  and  the  association  constants,  K, 


er  con- 


A,  o  = 


(56) 
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2.3.1,3    Mixed-Associating  Svste 


The  appropriate  expression  for  'S^  ^  i  n  a  mixed  associat- 
ing system  is  similar  to  equation  54  but  in  this  case  the 
summations  are  two  dimensional  and  each  term  in  the  summa- 
tion includes  an  additional  factor,  f,  ..,  the  weight  frac- 
tion of  A  in  species  A.B-, 


A,  o      7  r 

22fA.ijCij 


o^A,  ij  '^ij 


(57) 


(0  i  i  i  m,  1  <.  j  <  n) 


wher  e  S  . 


nd  C^i      are   the   sedimentation  coefficient   and 


ij  .0 

concentration  of  A.R. 
1  J 

Equation  57  can  also  be  re-expressed  in  terms  of  the  mo- 
nomer concentrations,  C  and  C„ ,  and  association  constants 
K^j  as  follows. 


A,  o 


IJSij,<,fA.ij4cJ 


(1  i  i  i  m,  1  i  j  i  n) 


'A.ijCiC^ 


(58) 


The  average  sedimentation  coefficients,  "S.  ^  jn  equa- 
tions 54-58  are  infinite  dilution  coefficients.  The  correc- 
tion for  hydrodynamic  dependence  is  the  same  as  in  nonasso- 
ciating  systems  (equation  52-53).   Thus,  for  instance. 


\  = 


A,  0 


^  '^  ^AA^AT  "^  '^AB^BT 
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(59) 


2.3.2    Co m£ u_t^ n^  L o£i j.  Cons 1 1  tnent  Concentrations 

The  boundary  monomer  concentrations  needed  to  compute  lo- 
cal average  transport  coefficients  are  extracted  from  bound- 
ary constituent  concentrations.  Boundary  constituent  con- 
centrations which  are  also  needed  for  computing  the 
hydrodynamic  dependence  of  the  local  average  sedimentation 
coefficients,  are  interpolated  from  box  center  constituent 
concentrations.  The  concentrations  of  a  constituent  at  its 
own  boundaries  are  computed  directly  by  linear  interpola- 
tion. 


'-AT.j  -  "AT.j-l  ■*■ 


^^AT,  j  -  <=AT,  j-lHrA,  j 


■^A.j-l' 


(r 


(60) 


A,j  "  "'A.j-l' 


The  cone  en t ra t i on  J  at  a  boundary  in  one  constituent 
track,  of  the  other  constituent  (eg.  C^^  at  r^  ■)  is  inter- 
polated between  the  two  box  centers  in  the  other  constituent 
track  that  bracket  the  position  of  the  boundary  being  con- 
sidered. 


.■'.rii'4'--v:\v,^  'K 
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r   ,,   .,    -  *°BT,ii  -  <=BT,n-l'  ("^A,  j  "  "^B,!!-!) 


wher  e 


*^B,n  "  ^B.n-l) 


'B,a-1  >  ^A.j  >  ^B,ii 


2 . 4    TABLE  ASSEMBLIES 

In  all  but  a  few  simple  cases  the  calculation  of  trans- 
port coefficients  from  local  constituent  concentrations, 
which  is  essentially  a  problem  of  extracting  monomer  concen- 
trations from  constituent  concentrations,  is  fairly  complex 
and  expensive  in  terms  of  computer  time.  A  typical  two  con- 
stituent 100  box  model  simulation  consisting  of  100  sedimen- 
tation transfers^  and  500  diffusion  transfers  would  require 
100  X  600  X  2  =  120,000  such  calculations  which  would  con- 
sume a  considerable  block  of  of  computer  time. 

An  alternative  approach  which  greatly  reduces  the  number 
of  calculations  is  to  assemble  a  table,  prior  to  the  simula- 
tion, containing  transport  coefficients  at  a  number  of  con- 
stituent  concentrations  over  an  appropriate  range,  and  to 
then  interpolate  transport  coefficients  from  the  table  as 
they  are  needed  during  the  simulation.  The  tables  are  pro- 
duced by  table  assembly  programs,  of  which  there  are  several 
types,  each  designed  to  create  tables  as  efficiently  as  pos- 


The  maximum  all ow able  diffusion  transfer  time  (At)  js  us- 
ually  about  1/5  as  long  as  the  sedimentation  transfer 
time.  Five  diffusion  transfers  are  thus  required  for  each 
sedimentation  transfer  (see  p.  94). 
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s  ibl e  for 


particular  class  of  systems.   (See  table  1.) 
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TABLE  1 
Varieties  of  Table  Assemblies 

I)  Self  Association  (nA  <=>  A  ) 

n' 

A)  A^  only  (n  =  2)  . 

B)  One  A^  aggregate  (n>2). 

C)  Multiple  A^  aggregates   (n  =  any  positive  inte- 
ger). 

II)  Miied  Association  (mA  +  nB  <=>  A  R  ) 

m  n'  • 

A)  Aj^Bj  only  ( n  =  m  =  1) 

B)  One   or  more  ^j^B^     aggregates   (n  =  any  positive 
i  nt e  ger  )  . 


integer) 


"2°n 


D)  One   or  more  A^B^  aggregates   (m  =  any  positive 
integer,  n  =  any  positive  integer). 
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2.4.1    Sel  f-As  soc  lat  Inn  sts  terns 

In  preparing  tables  for  self  associating  systems  entries 
are  made  at  regular  intervals  in  either  c  D.  or  S  ,  all  of 
which  are  functions  of  the  monomer  concentrations: 


=  1   +  ^K_jci 


(1  i  j  <.  n)  (62) 


?   = 


^o  -    5Sj,<,Kjci-l 


1  + 


h^r 


(2  <  j  <  n)  (63) 


0   = 


1    + 


(2  i  j  <  n)  (64) 


With  one  of  the  three  variables,  o  g,  or  'S  ,  fixed  at 
each  point  in  the  table,  the  other  two  are  computed  using 
equations  62-64.  This  involves  finding  the  value  of  c  that 
corresponds  to  the  value  of  the  fixed  variable  and  then  com- 
puting the  other  two  variables  from  c  Equation  62  can  be 
solved  directly  for  c^,  as  long  as  no  aggregates  higher  than 
dimer  are  present,  in  which  case  equation  62  is  a  quadratic. 


■1  Ml  -  4K,cJl/2 


1  = 


2K, 


(«5) 


Equations  (63)  and  (64)  can  be  solved  for  c.  as  long 
one  aggregate  (A   where  j  is  any  integer)  is  present. 
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is  only 


r  s.     -  s      1  ii/(j-i) 

A, o    "o    ^ 
"l  = (j  >  2)  (66) 

L^o  -  ^A.  0   Kj  J 


L"  -  °A.   JK  J 


(j  >  2)  (67) 


Finding  c^^  from  c^  when  aggregates  beyond  dimer  are  pres- 
ent, or  from  ?  or  D  when  two  or  more  aggregates  are  pres- 
ent, requires  a  binary  search,  whereby  the  correct  value  of 
c^  is  approached  through  a  logical  trial  and  error  proce- 
dure. The  general  scheme  of  the  binary  search  is  illustrat- 
ed in  the  algorithm  below.  The  first  test  value  of  c  is 
that  midway  between  its  maximum  and  minimum  values.  The 
maximum  possible  value  of  c^  is  c^  when  the  search  is  for  Cj^ 
at  some  value  of  c^,  and  c^,  when  the  search  is  for  c^  at 
some  value  of  'S  or  D.  The  scheme  shown  below  is  general;  y 
may  represent  either  c^,  s^  or  5.  Ytest  '*  ^^^  value  of  ei- 
ther c  5  or  D  computed  from  the  test  value  of  c,  and  v 
to  1  '       ^  o 


The  exact  value  of  c^^  corresponding  to  y^  is  usually  not 
found.  What  is  found  is  a  value  of  c^^  that  corresponds  to  a 
value  of  Yjgjj  that  is  acceptably  close  to  y^.    The  values 
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of  y  and   o^   stored   in  the   table   are   the   final   values  of 


test 


of  c. 


and  c 


(yteet-yo'-   yo  ^'  °°'  stored.    The  pairs 


1 ,  te  St  '■  J^  test-^o 


-J  and  y,  (y  =  c^,  5.  or  ^g),  stored  in  the  table  are 
thus  in  exact  agreement.  The  acceptable  difference  between 
^test  ^^^  -^o  ^^  arbitrary  but  should  be  small.  For  most  of 
our   programs   an   acceptable   difference   of   C    /10,000   is 

used. 


Al gor  i  thm 
I)      C 
II) 
III) 
IV) 
V) 


=  c 

top     max 


bot 

l.test  =  (Ctop  +  Ci,„t)  /  2 


test  "  f''=l,test 


^t 

"  '^test  -  yo'  <  C„„  /  10,000 

A)  Then  begin 

^-   "l  =  <=l,test 

2-     y   =   f(C^)   =   y^^^^ 

End. 

B)  Else  If  y^   ^ >y 

test^'o  ,■■ 

1 .  Then  begin 

"'  '^top  =  <=l,test 
b)  Go  to  III 

End. 

2.  Else  begin 

*'  ^bot  =  =l,test 
b)  Go  to  III 


ye 


End. 
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End  Algorithm 


The  operations  involved  in  creating  a  table  with  regular 
intervals  in  S   are  outlined  below. 


Al gor  i  thm 


I) 


De  t  erm  i  ne  5 


o.mai  '^o  "'  *^max'   '^^^     largest  value  of 


S^  that  will  occur  during  a  simulation  is  the  value 

of  ?„  when  c^  =  C„^^ 

A)  Find  c,  at  C 

1      max 

1.   Only  one  aggregate  (A.), 
a)  If  j  =  2  Then  Begin 
i)      c^  =  c,  +  K^o? 


ii) 


-1  +  (1  +  4K2Cj)l/2 
2K» 


End 
b)  If  j )2  Then  Begin 

i)      Requires  a  binary  search. 


=  t  =  "1  +  KjcJ 
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ii)     y  =   0^  =   f(cj^)   =  ^Kjcj 


End 


2.   More  than  one  aggregate. 
Requires  a  binary  search. 


O  S  =  Ci  +  ^KjcJ 


b)  y  =  0^  =  f(c^)  =  ^Kjc| 

B)  Compute  S^_^^^  from  c^  at  C„„. 
5        ^Ao  ^  ^SAj.oKjci-' 


II)    Compute  AS^,   the  increment   in  B^  between  adjacent 
table  entries. 


AS 


o,  max     o,  A 
n-1 


where  n  is  the  number  of  table  entries.    (n  =  100 
is  usually  sufficient.) 
Ill)   Compute   Cj^(i)   at   each   S^(i)   in   the   table   for 
(l<i<n) . 

^o(i'  =  Sao  +  (i-l)AS„ 


A)  Only    one    aggregate. 
1  . 

S  (i)  = 

1  .  K^jci-l(i) 

Cl(  i)  =   

L^o(i)  -  Sa.j.o  Kjl 

B)  More  than  one  aggregate. 
Requires  binary  search. 

1  +  ^KjcJ  Mi) 

IV)     Compute  c   and  D(i)  for  each  c  (i). 
A)  c^(i)  =  c^(  i)  +  ^KjC^(  i) 
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B> 


^jK.DicJd) 
^JKjci(i) 


End  Al gor  i  thm 
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2.4.2    Mi led  A^^ o ci a^^n^  Systems 

Tables  for   systems  undergoing  mixed  association  are  two 

dimensional.    Within  a  given  row  o{     the  table   C._  is  held 

Al 

constant  and  C^^  is  incremented.   Within  a  given  column  C, 
is  constant  and  C    is  incremented.    Each  position  (row  i, 
column  j)   in  the  table  corresponds  to  a  constituent  A  con- 


BT 


centr  a  t  i  on  C, 


constituent  B  concentration  Cr,rr     •■     and 


AT,  ij' 
the   average   transport  coefficients   of   the   system  at   those 

concentrations   D..        n„.,  .      n.         fi  "^ 

AA,  ij  '    "BB.ij'    "AB,ij'    "BA,  ij  •    ='Ao.ij' 

Bo  ii'  Surface  diagrams  of  some  sample  sedimentation  coef- 
ficient tables  for  mixed  associating  systems  are  shown  in 
figures  3  and  4 

The  constituent   concentrations  and  all  of  the  transport 

coefficients  are  functions  of  the  monomer  concentrations  C 

A 

and    Cg. 


'^AT    =    25fA,  ijKijCicJ 

"^BT  =  IhB.a^ijciH 

(0    <     i    <    m)     (0    <    j    i    n) 


(68a) 


(68b) 


Ao 


?L*A,^jfijf^j,o_CA_ci 
i^A.ijKijCicJ 

(0    i    i    i    m)     (0    i    j     <    n) 


(69a) 
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Figure  3:   Topologies  of  Sedimentation  Coefficient  Tables 
for  an  Uncooperative  AB^  System  (Kj=10*) 

See  p. 97  for  description  of  Uncooperative  AB 
systems.  ^ 


a.  ^A  vs.  Cf^j    and  Cbt 
'' •  ^B  vs.  C^j    and  Cgj 
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Figure  4:   Topologies  of  Sedimentation  Coefficient  Tables 
for  a  Completely  Cooperative  AB.  System 
(K,„=10l<') 

See  p. 148  fox  description  of  Cooperative  AB 

sy  s t  em  s . 


^^  vs.  C.J  and  Cp-T- 


'     '\    'y'y\\       ''•  ^B  ^''  ^AT  ^'^^    •^ 


BT 


:orrc- 
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wher  e 


and 
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*Bo    " vv 7"" —  (69a) 


'J     ij 
(0    i    i    i    m)     (0    <    j     <    n) 


AA    =  ("IaaIBB  ~  ^ABlBA)  ''  ^1 

"aB    =  ('^ABlAA  -  "^AAIaB)  '  ^1 

"bB    =  '"^BB^AA  "  "^BA^AB)  '  ^"i 

°BA    "  '"^BAIbB  "  '^BBlBA)  '  "^^ 


■^l    ~    "bbIAA    -    lAB^BA 


•■aA    =    ^^OijfA.ijKijCi    ^CJ 

(0    i    i     <    m)     (1    <.    j     <    n) 

''aB    =    &DijfA,ij'fijCicr' 

(0    <     i     <    m)     (1    i    j    i    n) 

Sb    =    ^^J^ijfA.ijKijCicJ-l 

(0    <     i    1    m)     (1    <    j    i    n) 

Sa    =    ^^iOljfA.ijKijCr^CJ 

(1    <     i    i    m)     (0    <    j    i.    n) 


(70) 


(71) 
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and 


(0  i  i  i  m,  1  <  j  <  n) 

Iab  =  l^fA.ijJKijCicr'- 

(1  <  i  1  m,  0  i  j  <  n) 

Ibb  =  I^fB.ijJKijCicJ-l. 

(1  <  i  1  m,  0  <  j  <  n) 
IbA  =  ^^B.iji^ijCi-lcJ. 

(0  i  i  i  m,  1  <  j  i  n) 


(72) 


f .  ■ 


i_  jj  =  iMW^  /  (iMW^  +  jMWg) 


B,  ij 


jMWg  /  (  iMW .  +  jMWp) 


(73) 


fi(o,i)  =  fgd-o)  =  0 

f^CLO)  =  fgCO,!)  =  1 

^ij  =  Ca.b.  /  cicJ 
1  J 

K(0,1)  =  K(1,0)  =  1 


(74) 


Here,  i  and  j  refer  to  the  numbers  of  A  monomers  and  B 
monomers,  respectively,  in  each  aggregate  species,  and 
should  not  be  confused  with  indices  of  table  entries  or  sim- 
ulation array  elements. 

The  summations  in  these  expressions  are  over  all  species 
^jB;  present  in  the  system  with  i  and  j  equal  to  or  greater 
than  the  values   given  below   each   equation.    For   example. 
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suppose  we  are  considering  a  system  which  contains,  in  addi- 
tion to  the  monomers,   A  and  B,   the  two  associated  species, 

"■■]_^2      ^^^     ^2^4*    "^^^   summations   in  this   case   include   only 
four  te  rm  s : 

1.  i  =  0,  j  =  1,  corresponding  to  monomer  B  (A„Bi). 

2.  i  =  1,  j  =  0,  corresponding  to  monomer  A  (A-Bn). 

3.  i  =  1,  j  =  2,  corresponding  to  species  A. b^  , 

4.  i  =  2,  j  =  4,  corresponding  to  species  A  B,  . 

In  some  of  the  summations,  one  of  the  monomer  species  is 
excluded  by  the  lower  limits  indicated.  For  instance,  (1  <. 
i  i.  m)  (0  <.  j  i  n)  indicates  that  i  is  greater  than  one  in 
all  of  the  species  to  be  included  in  the  summation,  which, 
means  that  the  monomer  B  term  ( i=0 ,  j=l)  is  not  included  in 
the  summa  t  i  on  . 

Thus,  for  example,  d^^  for  an  A+B+AB^+AjB^  system  could 
be  written  as  follows. 


''aA  =  ^liDA.B/A.  ijKijCi-^C^  = 

(1  i  i  i.  m)  (0  <  j  <.  n) 

°AfA(l'0)K(l,O)  +  D^B  fA(l,2)K(l,2)C^ 
\b/a(2''KK(2,4)CaC^ 


2D, 


(75) 


In  constructing  these   tables,   C 


AT 


and   C 


BT 


ire   fixed  at 


each  table  entry,   C^  and  Cg   are  computed  from  C^^     and  Cg^ 
using  equation  (68a)  and  (68b),  (this  usually  involves  bina- 
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ry  search)  and  S^^,  Sb,,,  D^a.  0^3,  ^^b-  ^-"i  ^BA  »"  computed 
from  the  monomer  concentrations  using  (69-74). 

The  complexity  of  a  mixed  association  table  assembly, 
which  is  related  to  the  number  of  binary  searches  involved, 
depends  upon  the  types  of  species  present.  When  only  one  ag- 
gregate, AB,  is  present,  (68a)  and  (68b)  become  simply: 


AT 


Ca    +     fA(AB)KABCACB 
=     C^     (1    +    fA(AB)K^BCB) 


(76) 


and 


Cgd    +    fB(AB)KABCA) 


(77) 


Solving      (76)       for      C^      and      substituting      for      C^      in      (77) 


AT 


1    +    fA(AB)KABCB 


{78) 


St  -   ^B 


^AT_fB^AB)KABCB 
1    +    fA(AB)KABCB 


(79) 
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(79)  can  be  written  as  a  quadratic  in  terms  of  C„  as  fol- 

B 

low  s  : 


Cbt(1    +    fA^ABCB'     = 


S<1    +    ^aKabCb)     +    CATfBKABCfl 


(80) 


BT  ■^  ^BT^^A^AB^b' 


'^B  "^  '^b^a^ab^b  *   '^at^b^ab'^b 


(81) 


^A^Ab'^B  *  L^  '^  '^AT^B^AB  ~  ^BT^A^AbJ^B  "  '^BT    '  ^^ ' 


where  f^  =  f^(AB)  =  M^/M^g  "-"^  ^B  =  ^flfAB)  =  Mg/M^g. 

Thus,  Cg  can  be  determined  directly  from  C,_  and  C„_, 
ing  the  quadratic  formula. 


_  !Ll!Al!?!AB_  +  _CBTfAKAB^  ^ 

"a^ab 
[(1  +  c^t^b^ab  "^  ^^bt^a^ab*  ~  '*^a^ab'^bt] 


"a^ab 


(82) 


C.  can  then  be  computed  directly  from  equation  (78). 

The  following  algorithm  outlines  the  table  assembly  pro- 
gram for  mixed  association  with  the  only  aggregate  species 
present  bei  ng  A.  fi^  , 


a''> 
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Al gor  i  thm 

I )     n    =    numb er  of  rows  and  col amns . 

"'  "at    =     C^Tmax  '     (■'"D 

"^'       -^^BT    =     CBx„a^  /     (m-1) 


IV) 
V) 


Repeat    from    i    =    1     to     (n+l) , 


A)     C 


AT 


■-AT 


+     (i-l)AC 


AT 


B)  Cgj    =    0 

C)  Repeat    from    j    =    1     to     (n+l). 

^-  St    =    Cbx    +    (j-DACgx 

2.  Compute  Cg  _  using  eq.-83. 

3.  Compute  C^jj  using  eq.  78. 

4 .  Comput  e   ^ 


A,o,ij'   ^B.o.ij'   Bj^j^jj,   D^^Bij, 
BB, ij '  ^BA  ii'  nsing  equations  69-74. 
End  repe  at . 

End  repeat. 
End  Algorithm 

A  binary  search  for  monomer  concentrations  is  needed  to 
construct  tables  for  systems  with  two  or  more  mixed  associa- 
tion products.  If  all  of  the  mixed  aggregates  present  con- 
tain only  one  subunit  of  one  of  the  constituents  (ie.,  if 
all  mixed  aggregates  are  from  A-^B^,  where  n  is  any  integer), 
and  if  there  are  no  self  association  products  in  constituent 
A  then  i  =  l  in  every  term  in  equation  (68a),  which  can  then 
be  rewritten  as  follows: 


'^A  = 


AT 


^f^d,  j)K(l,  j)C^ 
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(0  <  j  i  n)  (84) 


Substituting  the  right  hand  side  of  (84)  for  C,  in  equa- 
tion (  68b  )  gives: 


r 


AT 


St  =  }ltB^i.i)^(i.i)ci  -  . 

L2fA(l.J)K<l.J)cJj 
(0  i  i  i  1)  (0  i  j  i  n) 


Ij 


(85) 


If  there  are  no  self  association  species  in  B,   equation 
(85)  can  be  simplified  as  follows: 


■^BT  =  Cg  +  ^fgd-j) 


AT 


5fj^(l,j)K(l,j)CJ 
(0  i  j  i  n) 


(86) 


Whether  using  equations  (85)  or  (86),  a  binary  search  is 
required  to  determine  Cg  at  a  given  C^^^  and  Cg^.  The  object 
of  the  binary  search  is  to  find  a  value  of  C-  which  corre- 
sponds to  a  value  of  C^^  that  is  very  close  to,  but  not  nec- 
essarily equal  to  a  target  value  of  Cg^.  The  target  values 
°f  Cg^  are  not  stored  in  the  table,  but  serve  as  guidelines 
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for  constructing  a  table  that  is  nearly  evenly  spaced  in 
Cgj.  The  values  of  Cg^  that  are  finally  stored  in  the  table 
are  those  which  correspond  exactly  to  the  final  values  of 

Once  Cg  is  found,  C^  is  computed  directly  from  (84).  The 
general  scheme  of  an  A^^B^  table  assembly  with  a  binary 
search  is  illustrated  below. 


Al gor  i  thm 

I)     Compute  increments  for  C,„  and  C„t  .. 

Al       Bl, target' 


A) 


AC 


AT, max 


AT 


number  of  rows 


B) 


AC, 


BT,  max 


BT 


number  of  columns  -  1 


II)    Repeat  from  k  =  0  to  (number  of  rows). 

A)  C^T.k    =    kAC^T 

B)  Repeat    from    m    =    0    to     (number    of    columns), 

^  ■       St,  target    =    mACgT 

2.       Compute    initial     test    value    of    C 

^'     S,bot    =    0 

b  )     C„  -    r 

B,  top    -    l-BT,  target 


c)     C 


B,test    -    <<^B,top    ""■Cg^bot)     /    2 
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3.       Compute    test    value     of    Cg^    from    the    test    val- 
ue   of    Cg    using    equations     (85)     or     (86). 

Sl.test    =    5$fB<i'J)K(i.J)C^^ --—7 

2fA(l.J)K(l,j)cJ_t,^, 

(0    <     i     <    m)     (0    <    j     <.    n) 


Compare       C,, 


'■"^       '^BT.tETgef  " 


BT,  test 

St, test      '*      ^l"'"      enoagi      to      CBT.tngef 
store     Cg^^^^^^     and     Cg^g^^     as     Cg^d^.m)     and 
Cg(k,m),     respectively,     and    compute    and    store 

^^  'St,  target'^^BT,  test  I  <^BT,  targe  t'lO 

Then    begin 

a)     Cg^(k,m)     =    C^j^,,^, 


b)     Cg(t,m) 


B, test 


"AT 


c)     C^(k,m)     = 


5fA(l,j)K(l,j)cJ(k, j) 
(0    i    j     <.    m) 


End 
5  .       I  f     C, 


BT.test      1*     '°°     large,      set      Cg    ^.^^^     equal 
'°    S.test    ""^    compute    a    new    Cg    ,.^^^. 


Else    If    C„ 


BT,  test^^BT,  target 


Then    Begin 


•   ,  ,       ■  5  9 

*'  S,  top  =  ^B,  test 
b )  Go  to  2o . 

End. 

*•   ^^     St, test   *«  '°°  small,   set   Cg  ^^t  ei"al 
*°  ''B.test  """^  compute  a  new  Cg  test- 
Then  Begin 

*'  '^B,  bot  "  '^B,  test 
b)  Go  to  2c. 

End. 
7.   Compute   S^^  j^___,    Sg^  ,^„,    S^a.^^.    SBB.km- 

^AB.km'  SBA.km  using  equations  69-74. 
End  Repe  at . 
End  Repe  at. 
End  Algorithm 


^2"n' 


present,  equation  (68a)  is  a  quadratic. 


'^AT  =  cJf^(l,j)K(l,j)C^  +  c25f^(2,j)K(2,j)C^ 

(0  i  j  i  n) 


(87) 


:j 
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C.  = 


^f^(l, j)K(l, j)C^ 


2^fA(l,j)K(2,j)CJ 
(0  i  j  i  n) 


(88) 


Substituting  the  right  hand  side  of  this  expression  for 
Cj^  in  equation  (68b)  gives. 


St  =  ^2fB(i'J'K(i'J)cJ 


r-  5f^(l,j)K(l.j)cJ 


L  2^f^(l,j)K(2,j)CJ 


(0  i  i  1  m) 
(1  i  J  <  n) 


+  !*lfA'^'i!^'^'i!5B>^  +_4^fj^(l,j)K(2,j)CJC^j]l'2T 
25fj^(l,j)K(2,  j)CJ  J 

(0  <  j  i  n) 


(89) 


A  table  assembly  program  for  this  type  of  system  differs 
from  that  for  an  A^b^  system  only  in  that  equations  (88)  and 

(89)  are  used  in  place  of  equations  (84)  and  (85).   The  A.  R 

2  n 

assembly   program   is   suitable   for   systems   including   A- B 
species  and  self  associating  species  in  constituent  B    as 
well  as  A^  and  A2Bn  species,  but  cannot  be  used  for  systems 
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that  include  self  association  products  in  A  beyond  the  di- 


In  dealing  with  more  complex  systems  (A  b  ,  m>2)  it  is 
not  possible,  as  it  is  with  simpler  systems,  to  simplify  or 
combine  equations  (68a)  and  (68b)  in  any  way  that  will  allow 
a  direct  computation  of  one  of  the  monomer  concentrations. 
A  binary  search  must  therefore  be  employed  to  find  C,  and  C„ 
that  simultaneously  satisfy  equations  (68a)  and  (68b)  at  a 
given  C^^  and  Cg^.  In  this  case  target  values  are  used  for 
both  C,_  and  Cot  at  each  position  in  the  table.   Values  of 


C^  and  Cg  that  correspond  to  values  of  C,„  and  Cp^.  that  are 
acceptably  close  to  the  target  values  are  located  and  stored 
along  with  the  constituent  concentrations  to  which  they  cor- 
respond directly.  The  target  constituent  concentrations  are 
not  stored  in  the  table. 

The  organization  of  the  search  is  such  that  for  each  tri- 
al value  of  C^  a  value  of  Cg  is  found  (though  a  binary 
search  like  that  used  in  simpler  assembly  programs)  that 
satisfies  equation  (68b).  This  value  of  C_  and  the  trial 
value  of  C,  are  then  tested  in  equation  (68a).  If  equation 
(68a)  is  not  satisfied  the  trial  value  of  C,  is  adjusted  ap- 
propriately and  the  preceding  steps  are  repeated.  This  pro- 
cess continues  until  values  for  C^  and  Cg  that  satisfy  both 
equations  (68a)  and  (68b)  have  been  located. 

It  should  be  noted  that  this  table  assembly  is  capable  of 


dealing  with  A 


j^Bjj  and  A2Bjj  systems  as  well  as  more  complej 
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cases.  Simpler  assembly  programs,  however,  are  preferred 
whenever  possible  to  save  computing  time, 

2 . 5    TABLE  LOOKHaP 

In  order  for  a  simulation  program  to  extract  transport 
coefficients  from  a  table,  it  must  be  equipped  with  a  table 
look-up  routine.  The  lookup  routine  consists  of  (1)  a 
search  routine  which  locates  the  table  entries  with  constit- 
uent concentrations  that  bracket  the  boundary  constituent 
concent  rat i on ( s ) ,  and  (2)  an  interpolation  routine  which 
computes  the  local  average  boundary  transport  coefficients 
by  linear  int erpo 1  a t i on ( s )  between  the  bracketing  table  val- 
ues. 

For  single  constituent  (self  associating)  systems  the  ta- 
ble search  and  interpolation  routines  are  fairly  simple.  To 
obtain  5^,  g  at  a  boundary  where  C  =  Cj(j),  the  search  rou- 
tine inspects  each  successive  value  of  C  in  the  table  until 
it  finds  a  value,  C^(h),  greater  than  or  equal  to  C  (j). 
The  local  average  transport  coefficients  at  boundary  j  are 
then  computed  by  linear  interpolation  between  the  values  of 
the  transport  coefficients  at  C^(ii-i)  and  Cj(h)  which  brack- 
«t  C^(j). 

The  following  algorithm  illustrates  the  table  search  and 
interpolation  procedure  for  single  constituent  (self  associ- 
ating )  sy  st  em  s  . 
Al  gor  i  thm  .'!;■..'" 
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I )  Search  Rout  ine 

A)  h  =  2 

B)  If  C^(h)  >  Cj(j)  go  to  II 

C)  li  =  h  +  1 

D)  go  to  B 

II)  Int erpol at i on (T  ;  ^   or  B) 
T(j)  =  T(li-l)  + 

<T_(M  -  Kh-D)  (C^(j)  -  Cjdi-l)) 
(C^di)  -  Ct(h-l)) 


End  Al gor  i  thm 

The  table  look-up  for  a  two  constituent  system  not  under- 
going mixed  association  is  the  same,  for  each  constituent, 
as  the  single  constituent  system  look-up  described  above. 
For  example,  'S^  qCA.j),  the  local  average  sedimentation 
coefficient  of  constituent  A  at  boundary  j  of  constituent 
track  A,  is  obtained  by  interpolating  between  table  values 
A  o^^^  ^^^  ^A  o^^~^^'  which  correspond  to  constituent  A 
concentrations    C   (h)    and   C,_(h-1),    where    C._(h-1)    < 


"AT^ 


Cat(A,J)  <  Cj^t(^> 


^A,o(A,j)  =  S^  „(h-l)  + 


<^A  „(t) 


_Sa^o<''-i^1_(<Cax(a^j)  -_5at(5-i>> 

(Cj^T(h))  -  (Cj^j(h-l)) 


(90) 


C^(A, j)  =  D^(h-l)  + 
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(91) 


(Cg^(h-l))  -  (Cb^(B,j)) 


(92) 


BgCB.j)  =  DgCh-l)  + 

(Cgdi)  -  Dgdi-D)  (Cb^(B,j)  -  CB^(h-l)) 


(Cg^di) 


CB^(h-l)  ) 


(93) 


When  mixed  associations  are  involved  the  table  look-up 
procedures  become  more  complex.  The  table  transport  coeffi- 
cient  arrays,  ?^,  Sg,  5^^,  D^b.  633  and  034,  are  all  two  di- 
mensional.  Each  member,  (i,j),  of  these  arrays  corresponds 
to  a  pair  of  constituent  concentrations  in  the  C,_  and  Cox 
arrays. 

For  the  simplest  mixed  associating  system,  ( AB ) ,  it  may 
be  recalled  that  the  construction  of  the  table  does  not  re- 
quire a  binary  search,  and  that  the  values  of  constituent 
concentrations (C^^  and  Cg^)  for  each  table  entry  may  there- 
fore be  selected  to  create  a  table  with  rows  varying  by  ex- 
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act  increments  in  C^^  and  columns  varying  by  exact  incre- 
ments in  Cg^,  Cat  is  constant  within  a  given  row  and  C„-,  is 
constant  within  a  given  column.  The  C^^  and  Cgj  arrays  in  a 
table  of  this  type  are  one  dimensional. 

The  following  algorithm  outlines,  in  general  form,  the 
procedure  used  by  the  table  look-up  routine  to  extract  a  lo- 
cal average  transport  coefficient,  T.  .  or  T„  .  (T.  .  de- 
notes  a  transport  coefficient  at  boundary  j  in  the  constitu- 
ent A  track  and  Tg,j  denotes  a  transport  coefficient  at 
boundary  j  in  the  constituent  B  track  (ie.  T, ,  j  =  Sa(A,j), 
°AA<A'J'  »'■  "ab'^'J'  '""l  Tg,j  =  Sb(B,j),  5gB(B,j)  or 
BA^^'J^^^'  corresponding  to  boundary  constituent  concentra- 
tions, C^j(A,j)  and  C^jiA.j)  or  Cf^j(B.  j)  and  Cbj(B,j),  from 
a  table  of  the  type  described  above. 


Al gor i  thm 

I)     Table  search  routine 

A)  Locate  the  two  members  of  the  C,_,  table  array 
'C^I'(h-l)  and  Cj^j(h))  whose  values  bracket  the 
boundary  j  constituent  A  concentration  (C,_(j)). 
1  .   h  =  2 

2.   If  C^j,(h)  >  C^x<J'  SO  to  B 
3  .   h  =  h  +  1 
4 .   go  to  2 

B)  Locate  members  of  the  table  C^^  array  (CB^Ck-l) 
^°"  *-B'j'(k))  whose  values  bracket  the  boundary 
constituent  B  concentration  C„_(j). 


II) 
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1.   k  =  2 

2-   if  %j(.k)    i_    CgjCj)  go  to  II 

3  .   k  =  k  +  1 

4  .   go  to  2 
Interpolation  routine 

A)  Interpolate   the  value   of   the   transport   coeffi- 


c lent   at 


AT 


C^j(h-l)   and  C^j 


^BT 


(j). 


T(h-l,j),  between  T(h-l,k-I)  and  T(h-l,k). 


T(h-l,j)  =  T(h-l,k-l)  + 

(T(h-l,k)  -  T(h-l,k-l))  (Cgj(j)  -  CBT(k-l)) 


<Cg^(k) 


Cg^(k-l)) 


B)  Interpolate  the  value  of  the  transport  coeffi- 
cient at  C^^  =  C^^(h)  and  Cg^  -  Cgj(j).  T(h,j), 
between  T(h,k-1)  and  T(h,k). 

T(h, j)  =  T(h,k-1)  + 

I!!'':'''  "_^'^:!!I^!!_!fBT(Ji_:_5BTik:i>' 

(Cg^(k)  -  Cgj(k-l)) 


C)  Interpolate  the  value  of  the  transport  coeffi- 
cient  at  C^^  =.  c^^  ( j  )  ^nd  Cg^  =  CgxtJ).  T(J)- 
between  T(h-l,j)  and  T(h,j). 
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T(j)  =  T(li-l,j)  + 


End  Al gor  i  ttm 


It  will  be  recalled  that  tlie  construction  of  tables  for 

systems  with  mixed  aggregates  of  the  form  AB  .  where  n  is  an 

n 

integer  greater  than  one,   requires  a  binary  search  to  find 

corresponding  pairs  of  Cg^  and  Cg ,  and  that  therefore,  while 

the  target  values  of  Cg^  are  exactly  incremented  in  each  row 

and  constant  within  each  column,   the  values  of  C__,  that  are 

finally  stored  in  the  table,   those   that  correspond  exactly 

to  the   stored  transport   coefficients,   being  only  close   to 

the  target  values,   are  not  exactly  incremented  within  rows 

or  constant  within  columns.   The  value  of  C    is  thus  unique 

BT  ^ 

at  each  position,  (h,k),  in  the  table  so  that  C_-,  requires  a 
two  dimensional  array. 

As  long  as  no  mixed  aggregates  containing  more  than  one 
molecule  of  constituent  A  (and  no  self  aggregates  of  A  larg- 
er than  A^ )  are  present,  a  binary  search  is  not  required  for 
^j^j    and  the  C^j    array  can  be  kept  one  dimensional. 

The  following  algorithm  illustrates  the  general  procedure 
for  extracting  transport  coefficients  from  tables  of  this 
type  (ie.  a  two  dimensional  Cg^  array  and  a  one  dimensional 
''AT  array)  . 
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Algorithin  , '- ' 

I)      Search  routine 

A)  Locate  members  of  the  table  C^^  array  (C.jCh-l) 
^^d  ^AT^^^^  whose  values  bracket  the  boundary 
constituent  A  concentration,  ^at'J'- 

1  .   h  =  2 

2.   if  C^^(h)  >  C^x(j)  go  to  B 

3  .   h  =  h  +  1 

4 .   go  to  2 

B)  Locate   members   of   the 


^gj  array  in  row  h, 
(Cg^(h,a-1)  and  Cg™(h,a),  whose  values  bracket 
the  boundary  j  constituent  B  concentration 
(CgjCj))  . 

1  .   a  =  2 

2  . 

3  , 

4  , 


If  Cg^(h,  a)  >  Cg^(j)  go  to  C 
a  =  a  +  1 
go  to  2 


C)  Locate   members   of   th 
and 
bracket  Cg^(j)  . 


(Cg^(h-l,b-l) 


■nj     array      in      row      h-1, 
Cg™(h-l,b),         whose        values 


1 
2 
3 

4  .   go  to  2 
II)    Interpolation  routine 


b  =  2 

if  Cg^(h-l,b)  >  Cgj(j)  go  to  II 

b  =  b  +  1 
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A)  Interpolate  the  value  of  the  transport  coeffi- 


cient  at 


^BT 


'-BT 


(j). 


T(h-l,j),  between  T(h-l,b-l)  and  T(h-l,b). 


T(h-l,j)  =  T(h-l,b-l)  + 
(T(h-l^b)  -  T(h-l,b-l))  (C3j,(j)  -  CBT(h-l,b-l)) 
(Cg^(h-l,b)  -  CgjCh-l.b-l)) 

B)  Interpolate  the  value  of  the  transport  coeffi- 
cient at  C^^  =  Cf^.j.(h)  and  Cg^  =  Cg^(j),  T(h,j), 
between  T(h,a-1)  and  T(h,a). 

T(h, j)  =  T(h, a-1)  + 

(TU,a)  -  T(h,a-1))  (Cg^Cj)  -  Cg^(h,a-1)) 
(Cg^Ch.a)  -  Cg^Ch.a-D) 

C)  Interpolate  the  value  of  the  transport  coeffi- 
cient at  C^j  =  C^T<J)  »■"»  ^Bx  =  Cbx(J).  T(j), 
between  T(h-l,j)  and  T(h,j). 

T(j)  =  T(h-l,j)  +      .  .,        ,  j 

L'^-t'O.    "_!!^'^:i'l_'^AT(J'  -  CAx(t-l)) 


End  Al gor  i  thm 

Generating  tables  for  more  complex  mixed  associating  sys- 
tems  (ie.   systems  with  mixed  aggregates  of  the  form  A  n 

m  n' 
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where  m  >  1  and  n  >  1,  or  systems  with  only  AB  mixed  aggre- 
gates but  with  one  or  more  A^  aggregates  with  m  >  2)  in- 
volves a  binary  search  for  C^^  -  c^  pairs  as  well  as  for  Cgj 
Cg  pairs.  As  a  result  the  values  of  both  C.^  and  Cpj  are 
unique  at  each  position  in  the  table  and  two  dimensional  ta- 
bles are  necessary  for  both  constituent  concentrations. 

Because  of  the  irregularity  of  tables  of  this  type,  the 
look-up  routine  that  is  required  is  somewhat  more  involved 
than  those  required  for  simpler  systems  and  the  transport 
coefficients  extracted  from  them  are  likely  to  be  slightly 
more  erroneous  than  those  extracted  from  the  more  orderly 
tables. 

As  with  simpler  mixed  associating  systems,  four  points 
are  located  in  the  table  for  interpolation,  and  these  four 
points  define  a  quadrilateral  which  contains  the  point 
'''At'J''  '-rt'J'''  "•'Bre  the  interpolation  is  to  take  place. 
However,  whereas  with  simpler  tables  it  is  always  possible 
to  select  points  which  allow  interpolation  (as  opposed  to 
extrapolation)  along  all  four  sides  of  the  quadrilateral, 
this  is  not  the  case  with  the  doubly  irregular  tables. 

Portions  of  each  type  of  mixed  association  table  are 
shown  in  figure  5.  In  the  first  two  types  of  tables  (AB  and 
AB^)  any  point,  (C^j(j),  Cg.p(j)),  will  lie  within  a  quadri- 
lateral in  such  a  way  that  (1)  Cg^(j)  will  fall  between  the 
Cgj  values  of  the  two  points  that  define  the  top  of  the 
quadrilateral,   (2)  C^j(j)     will   fall  between  the  Cg^  values 
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of  the  two  points  that  define  the  bottom  of  the  quadrilater- 
al' (3)  C^^(j)  will  fall  between  the  C^j  values  of  the  two 
points  that  define  the  left  side  of  the  quadrilateral  and 
(4)  PCAT(J)  will  fall  between  the  C,„  values  of  the  two 
points  that  define  the  right  side  of  the  quadrilateral. 
Furthermore,  these  four  conditions  can  always  be  satisfied 
(in  the  simpler  tables)  without  extending  any  edge  of  the 
quadrilateral  over  more  than  one  table  division  in  either 
constituent  concentration.  (eg.  in  figure  5-b  these  condi- 
tions are  met  by  selecting  table  points  (2,3),  (2,4),  (3,3) 
and  (3,2)  . 

Figure  5-c  depicts  a  situation  which  can  arise  while  in- 
terpolating from  the  double  irregular  table  of  an  A  B  sys- 
tem. The  point  C^^' J  '  ' '^BT  *  J  '  **  included  in  the  quadrilat- 
eral (2,3),  (2,4),  (3,4),  (3,3).  However,  Cg^(j)  is  less 
'!>»■>  Cg^(2,3)  so  that  the  first  condition  listed  above  is 
violated  and  an  extrapolation  will  be  necessary  along  this 

side  of   the  quadrilateral   with  respect   to   C„_.    Moreover, 

B  T ' 

there  are  no  four  points  which  satisfy  all  four  conditions 
without  extending  an  edge  of  the  quadrilateral  through  two 
table  divisions  in  one  constituent  concentration  or  the  oth- 
er. Although  the  irregularity  of  the  table  in  figure  5  has 
been  greatly  exaggerated  for  illustration  (in  a  typical  ta- 
ble the  deviation  from  r e c t angnl ar i ty  would  be  barely  noti- 
cable)  this  problem  is  real  and  will  occur  in  a  small  per- 
centage   of    the    interpolations.     However,    since    the 


72 
irregularity  is  very  small  in  relation  to  the  distance  be- 
tween adjacent  table  entries,  it  is  not  a  serious  problem 
and  a  sufficiently  accurate  "interpolation"  should  always  be 
possible. 

The  A^Bjii  search  routine  is  responsible  for  locating  an 
appropriate  "rectangular"  quadrilateral  for  interpolation  at 

*'At'J''  '-Bx'J'  (rhombic  quadrilaterals  are  not  considered  ). 
A  first  guess  is  made  by  locating  the  table  entries  in  the 
first  row  whose  Cg^  values  bracket  Cg^(j),  and  the  entries 
in  the  first  column  whose  C^^  values  bracket  C^x(j).  For 
instance,  if  Cg^(i,t-i)  <  c^^(j)  <  C^^(l,i.),  and  CAT(h-]..l) 
^  ''at''''  ^  C^^(h,l),  then  since  the  table  is  irregular,  but 
only  slightly  so,  ( C^^( j ) , Cg^( j ) )  should  be  included  either 
within  the  quadrilateral  (h-l,k-l),  (h-l,k),  (h,k),  (h,k-l) 
or  within  one  of  eight  quadrilaterals  surrounding  it. 

The  search  routine  then  sequentially  checks  each  edge  of 
the  test  quadrilateral  to  determine  whether  the  point 
'"at'-'''  ''Bt'J'  H"  °°  "**  inside  or  outside  of  the  quadri- 
lateral along  that  edge.    If  C^^Cj),  Cb^(j)  lies  outside  of 


In  addition  to  the  rectangular  quadr i 1  a t era  1 ( s )  several 
(4-12)   rhombic   quadrilaterals   may   contain   the   point 

''At'J''  *-Bx'J''  Although  one  or  more  of  these  rhomboids 
may  be  more  optimal  for  interpolation  than  the  rectangular 
quadrilateral  (the  most  ambitious  table  extraction  proce- 
dure might  be  to  average  the  interpolated  values  from  all 
(as  many  as  16)  of  the  quadrilaterals  containing  C,„(j), 
''Bt'J'  ^^^  "first  rectangular  quadrilateral"  method  should 
give  adequately  accurate  results,  is  much  simpler,  and 
would  consume  less  computer  time  than  procedures  that  con- 
sider other  quadrilaterals. 


Figure  5:   Types  of  Table  Organization 
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Regular  Table:  AB  systems  only 
C.rp  and  Cgrp  evenly  spaced.     ;,■ 

Irregular  Table:  AB   systems  only 

''AT  evenly  spaced,  Cg_,  spaced  unevenly. 

Double  Irregular  Table:  A  b„  systems 
'-^-j-  and  Cgj  both  spaced  unevenly. 
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the  quadrilateral  on  a  given  edge  then  the  test 
quadrilateral  is  moved,  horizontally  or  vertically,  one  ta- 
ble division  in  the  appropriate  direction.  Checking  and  re- 
positioning continue  until  the  correct  quadrilateral  is  lo- 
cated. 

The  following  algorithm  illustrates  the  procedure  used  by 
the  search  routine  to  locate  a  rectangular  quadrilateral,  in 
a°  *n^m  t^ble,  that  includes  the  point  ( C, ^ ( j ) , Cg j( j ) ) . 


Al gor  i  thm 

I)      Search  Routine. 

A)  Locate  table  entries  in  column  1  whose  C 


AT 


val- 


ues  bracket  C^j(j)  (ie.  C4j(h-l,l)  and  Cj^^(h,l) 


1.  h  =  2 

2.  If  C^^(h,l)  >.  C^t(j)  Then  Proceed 
Else  Begin 

a)  h  =  h  +  1 

b )  Go  to  2  . 

End  ■       ■   '    ' 

B)  Locate  table  entries  in  row  1  whose   C„_  values 

BT 

bracket   Cg^(j)   (ie.   Cg^d.k-l)   and   Cg^d.k): 

CBj(l,k-l)  <  CgjCj)  <.  CBj(l,k)). 
1  .   k  =  2 

2-   If  Cg^d^i)  >  Cgj(j)  Then  Proceed 

; El se  Begin 
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a)  k  =  k  +  1 

b)  Go  to  2 
End 

C)  Inspect  each  edge  of  the  quadrilateral  (h,k-l), 
(h-l,k-l),  (h-l,k),  (h,k)  and  adjust  h  and  k,  if 
necessary,  until  a  suitable  quadrilateral  is  lo- 
cated. 


f  '  >' 


If  '^At(J)  <  CiT(li.k-l)  + 
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(C3j(j)  -  Cg^di.k-D) 


Then  Proceed 
Else  Begin 


Cg^(h,k)  -  Cg^(h,k-1) 


a)  h  =  li  +  1 

b)  Go  to  3 
End 

If  <^At'J'  ^  C^.j.(h-l.i.)     + 


(Cg^Cj) 


Cg^(h-l,k-l)) 


C^j(h-l,k)  -  Cj^^(h-l,k-l) 


Then  Pro  ce  e  d 
Else  Begin 


BT 


(h-l,k)  -  Cg^(h-l,k-l) 


a)  h  =  h  -  1 

b )  Proceed. 
End 

If  St'J'  ^  Cg^Ch-l.k-l) 

"bt 


(C.T-Cj)  -  C._(h-l,k-l))  — 


(h,k-l) 


"-BT 


(h-l,k-l) 


'-AT 


Then  Proceed 
Else  Begin 


AT 


(h,k-l)  -  C^.f(h-l,k-l) 


a)  k  =  k  -  1 

b)  Go  to  C 
End 
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I^  ^Bt'J)   <  CBT(h-l,k) 


Then  Procee  d 
Else  Begin 


Cg-j.(h,k)  -  CgjCh-l.k) 
C._(h,k)  -  C.T.(h-l,k) 


a)  k  =  k  +  1 

b)  Go  to  C 
End 

End  Algorithin 

The   following   algorithm   illustrates   the   interpolation 

procedure  for  A  B   systems, 
n  m  ^ 


^^i  St  «'   c^T- 


Al gor  i  thm 

I)  Interpolate     the    values    of     the     transport     coefficient 

ilong    line     (h-l,k-l),     (h,k-l). 

A)    T(j,k-1)     =    T(h-l,k-l)     + 

Cjx'J)     -    CiT.(h-l,k-l) 
(T(h,k-1)     -    T(h-l,k-l)) 

C^x(h,k-1)     -    C^j(h-l,k-l) 


■*'     St'J'''^"!'     "    Cgx(h-l,k-l)     + 


Cat(J)     -    C,T,(h-l,k-l) 


(C       (h,k-l)     -    Cm.(h-l,k-l)) --- 


Cf^y{h..i.-l)-C^^{i,-X,i.--i.) 


II)    Interpolate  the  values  of  the  transport  coefficient 
«°d  Cgj  at  C^iCj)  along  line  (h-l ,  k)  ,  (  h,  k  ) 
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A)T(j,k)=T(li-l,k)  + 

C.T<J)  -  C.T(i.-l,k) 
(T(li,k)  -  T(li-l,k))— " 

C^j(li,k)  -  Cf^jlh-l.i.) 


B'  Cg^(j,k)  =  C3j(h-l,k)  + 

C.T<J'  -  C.^(ii-l,k) 
(Cg^di.k)  -  Cgj(li-l,k)) 

C^j(h,k)  -  Cj^j(li-l,k) 


III)   Interpolate    transport    coefficients    at    C._(i) 

AT  ^  J ' * 

Cbt(j). 

T(j)  =  T(j,k-1)  + 

Cbt(J)  -  C„T-(j,k-l) 
(T(j,k)  -  T(j,k-1))— 5- ^i 

Sl'J''''  "  CB^(j,k-l) 

A)  End. 
End  Al gor  i  thm 

2 . 6    REzINDEXING 

When  the  boundaries  are  shifted  in  a  sedimentation  trans- 
fer, the  entire  array  of  boxes,  the  simulation's  frame  of 
reference,  is  displaced  downward  in  the  cell.  Space  is  cre- 
ated between  the  initial  top  of  the  array,  r  _  and  the  first 
boundary  while  the  last  boundary  migrates  past  the  inital 
bottom,  r^^^^  (see  figure  6)  . 
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r,    Tj 


Figure  6:   Boundary  Movement 

If  this  process  were  all  owed  to  continue*  unchecked, 
through  a  large  number  of  sedimentation  transfers,  a  large 
part  of  the  array  would  wind  up  outside  the  cell. 

To  keep   the   simulation  arrays   from   sliding  out   of   the 

7 
cell,  new  boxes  are  usually  (but  not  always  )  added  to  the 

top  of   the   array   as   space   becomes   available,    A  constant 

number  of  boxes   is  maintained  by  eliminating  one  box  from 

the  bottom  of  the  array  for  every  box  added  to  the  top. 

When  the   space  between  r    and  r-   is  large   enough   (ie, 

when  r^  -  r   >  Ar),  a  new  box,  Ar  wide,   is  created  at  the 
1      o 

top  by  inserting  a  new  boundary.    The  new  boundary  becomes 
the  first  in  the  array  and  is  assigned  a  position 


Occasionally  it  is  desirable  to  use  a  very  short,  finely 
divided  array  (especially  for  simulating  systems  with 
large  species  that  generate  sharp  boundaries).  In  these 
cased  simulation  arrays  are  allowed  to  slide  and  re-index- 
ing is  bypa  s  se  d. 


new.l 


Ar 
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(94) 


In    most     cases     the     initial     solute    boundary     is     far    enougli    be- 
low    the     top     of     the     array     that     the     concentration     of     the     new 

box    can    safely    be     set    to    zero     (c  .=00) 

new.l    "  •"  '  • 

To  accommodate  a  new  box  at  the  top,  the  old  boxes  below 
must  be  re-indexed,  and  the  last  box  discarded.  To  avoid 
propagating  r^  and  c^  to  the  bottom  of  the  array,  reindexing 
must  be  carried  out  in  two  steps: 

1.  The  new  r  and  c  values  are  stored   with  the  appropri- 
ate new  indices  in  "temporary  holding  arrays"  (r 

new 

and  c    ) . 
new 

2.  The  r  and  c  array  members  are  assigned  the  values  of 

*^^^    'new  ^"^^    "nev    »rray  members. 
The  following  algorithm  adds  a  box  to  the  top  and  rein- 
dexes  the  simulation  arrays  (if  and  only  if  there  is  room). 


Al gor  i  thm 

I)      If  r(l)  -  , 
Then  Begin 


>  Ar 


A)  C^,„(l)  =0.0 

C)  Repeat  from  i  =  2  to  n+1 


1-  '^new'i'  =  "^(i-1) 

2-  C^ew'i)  =  C(i-l) 


i" 


.V.,. 


End  Repeat, 
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D)  Repeat  from  i  =  1  to  n+1 

1  .   r( i) =r    ( i) 

new  ^  ' 

2  .       C(i)  =  C    ( i) 
End  Repe  a  t ; 

End  .  ^'  - 

Else  Continue 
End  Algorithm 

2 . 7    REGULATI NG  BOX  SIZES 

Because  the  boundaries  sediment  at  different  rates,  the 
original  box  sizes  will  change  during  the  sedimentation  rou- 
tine. Some  boxes  may  become  much  wider  than  Ar ,  others  much 
narrower.  Boundary  spreading  can  be  expected  with  associat- 
ing systems  in  areas  where  the  gradient  profile  is  steep  and 
will  be  most  pronounced  when  the  concentration  there  corre- 
sponds to  a  rapidly  rising  region  of  the  S  vs.  C  curve  for 
the  system.  Compression  of  boxes  occurs  when  the  concentra- 
tion gradient  is  steep  and  the  sedimentation  coefficient  is 
decreasing  with  concentration  (due  to  hydrodynamic  depen- 
dence) . 

Boxes  that  are  too  wide  lead  to  serious  errors  in  inter- 
polation. Boxes  that  are  too  narrow  are  overly  sensitive  to 
small  errors  in  shifting  the  boundaries.  In  a  very  narrow 
box,  the  error  in  volume  that  would  result  from  a  slight  er- 
ror in  moving  one  of  its  boundaries  would  represent  a  sig- 
nificant fraction  of  the  boxes  total  volume  and  would  there- 
fore lead  to  a  sizeable  error  in  the  concentration. 
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In  order  to  avoid  these  problems  box  sizes  must  be 
regulated.  Thus,  following  each  sedimentation  event,  the 
sizes  of  all  of  the  boxes  are  examined,  and  any  boxes  that 
are  narrower  than  0.8  Ar  or  wider  than  2  Ar  are  eliminated. 

2.7.1    Box  S£lil  ti  n^ 

If  a  box  is  found  to  be  too  large  it  is  split  at  its  mid- 
point into  two  halves.  The  concentrations  in  the  upper  and 
lower  half-boxes,  C^.^  and  C^^^^,  are  calculated  using  ex- 
pressions that  incorporate  conservation  of  mass.  According 
to  the  Law  of  Conservation  of  Mass: 


"j'^j+l  -  '^V     ^  t^bot'-j+l  -  '-j)  ^  C,<,p(r 


-2 


'i' 


<95) 


where  c  is  the  concentration  of  the  original  box,  of  which 
^j_2  ^^^  ^1  ^^®  ^^®  upper  and  lower  boundaries,  and  r.  is 
the  midpo  int . 

Assuming  a  linear  gradient  within  box  j,  the  concentra- 
tions (at  the  midpoints)  of  the  upper  and  lower  box  halves 
are  related  to  one  another  as  follows. 


bot 


-top  *    gj'ir 


(96) 


5j      is     interpolated     between     gj +j^     at     r  j +j^     and     gj     at     rj      as 
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where  g   is  the  gradient  at  the  midpoint  of  box  j   and  Ar  i 
the  distance  between  the  midpoints  of  the  two  box  halves. 
Substituting  the  right  hand  side  of  equation  96  for  C. 
in  equation  95  gives 


bot 


'j<^^i  -  ^V 


(C 


top    -    ij^'X'j.l    -    ^j'^    C,<,p(rj    -    r^ 


Factoring    this    gives 


(98) 


Cj(r.^l    +    -jX-j+l    -    -j 


(C 


top    +    «jAr)(rj^i    +    ^jX-^j+l 


'^top'^j    ^    -jX^j    -    'j) 


-    r)     + 


(99) 


Then,     s  i  nee 


f ol 1 ow  s  : 


•j     ■"    <«j+l    -    «j)- 


(r. 


'P 


(r 


J+1 


^P 


(97) 


where     gl     and     i't  +  i     are     the     gradients     at     the     midpoints     be- 
tween   adjacent   ^oz    centers     (see    eq    20-21), 
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J+1  ~  ''j 


equation  99  can  be  simplified  to  give 

(101) 
Re-grouping  and  re-arranging  gives 


_  ^^jt'j.l  -  -j)  -  iiAr(r   1  +  -r  ) 

•top Z (102) 

^j+1  ^  ^j  ^  2-rj 


Since  2r  =  "^i+l  "*■  ^i'  ^^^  denominator  in  equation  102 
can  be  rewritten  as  2 ( r ^ ^^  +  r . ) .  Making  this  substitution 
and  simplifying  gives  "'■''' 


Fina 1 ly ,  since 


'top  °  °j  "  8j:7 <103) 

2<^j+l  +  Tj) 


'j+1  +  Tj  =  2rj  +  Ar   :/h  ' 

,1.1- 


and 


'j-i  '  ':   =  2^j. 


equation  (103)  can  be  rewritten 
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g.ArT     Arl 


(104) 


Equation   (104)   is  used  to  calculate   C^     which   is  then 

top 

used  in  equation  (96)  to  compote  C,   ..   xhe  arrays  are  then 

Dot  ^ 

re-indexed  to  include  the  new  boundary.   The  top  half  of  the 

split  box  becomes  box  j   with  a  concentration  of  c    (i)   = 

new  ^  ■'  ' 

Cf       .       The  bottom  half  becomes  box  j+1  with  c    (i  +  1)  =  C, 
top  •>  new'-''bot 

and  an  upper  boundary,  the  new  boundary,  at  r    (i+1).   The 

J  new  ^  ■'       ' 

indices  of  the  concentration  and  boundary  position  arrays 
are  incremented  by  one.  The  last  members  of  the  original 
arrays  are  not  copied  into  the  re-indexed  arrays  and  are 
thereby  discarded.  Figure  7  illustrates  box  splitting  and 
reindexing. 

The  following  algorithm  reindexes  the  r  and  c  arrays  when 
aboxis  split. 


Al gor i  thm 
I) 
II) 


5     (  i  )  =  C 
new^  J  I  v,^ 


op 


r(J  +  l)  =  C^^t 


before  box  splitting: 


)"'        j  j*l      j*2  n-1         n       ml 
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after  splitting   box  j; 


j"1       i  j*1      i*2     j.3 


n         n*1 


Figure  7:   Box  Splitting  and  Re-indexing 


m'   r^^^(j+l)  =  [(r(j)  +  r(j  +  l)]/2 
IV)    Repeat  from  i  =  (j+2)  to  (n+1) 

A)  r^^^(i)  =  r(i-l) 

B)  c^^^(i)  -  <:(i-l) 
End  repeat. 
Repeat  from  i  =  1  to  (n+1) 


A)  r(  i)  =  r    (  i) 
new  ^  ' 


V) 


B)  o(i)  =  c 
I 

End  r  epe  a  t . 
End  Algorithm 


,(i) 


2.7.2    Bo I  SquashinR 

When  a  box  is  found  to  be  too  narrow,  it  is  also  split  at 
its  midpoint  and  the  concentrations  in  the  two  halves  are 
computed  as  described  above  (equation  (104)  and  (96)).  In 
this  case  however,   the  upper  and  lower  halves  are  added  to 
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the  boxes  immediately  above  and  below  them,  and  the  original 
boundaries  of  the  undersized  box  are  dissolved,  ie.  three 
boxes  become  two,  (see  figure  8).  The  concentrations  of  the 
resulting  combined  boxes  are  computed  directly  from  the  con- 
servation of  mass  expressions: 


,u-i 


!Vp'<^J 


.!i!.!_!i=i!fL:.!izi!. 


<^1  -  ^J-^ 


(105) 


(Cv„.)(r? 


,(j)  = 


bot'^^j+1 


=  j+l<^^2 


(r2 
J+2 


tl06) 


The  arrays  are  then  re-indexed.  The  boundary  that  was 
added  at  the  center  of  the  original  undersized  box  becomes 
boundary  j  since  the  original  boundary  j,  the  upper  boundary 
of  the  undersized  box,  was  eliminated.  Since  the  original 
boundary  j+1  was  also  removed,  the  indices  of  all  array  mem- 
bers below  the  new  boundary  j  must  be  reduced  by  one.  Also, 
because  there  was  a  net  loss  of  one  boundary  (and  one  box) 
in  this  transaction,  a  new  boundary  must  be  added  to  the 
bottom  of  the  array.  The  new  box  is  usually  given  the  same 
width  and  concentration  as  the  box  immediately  above  it. 
Figure  8  illustrates  box  squashing  and  re-indexing. 
The  following  algorithm  re-indexes  r  and  c  arrays  following 
boxsquashing. 


before   box  squashing: 


i-1 


I       J* 


after  squashing   box  j: 


n.1 


n-1 


Figure  8:   Box  Squashing 
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Al gOT  i  thm 

^'     "^uew'J)  =  [r(j)  +  r(j  +  l)]/2 
II)    Repeat  from  i  =  j+1  to  n-1 


B)  "new'*)  =  c(i+l) 


III) 
IV) 
V) 
VI) 


End  Repe  at . 


"^new'"^!)  =  "^newf")  "^  ''^new'")  "  "^new*"-!" 


,(i) 


,(n-l) 


new ' - '     "new ' 
Repeat  from  i  =  1  to  n 


A)  r ( i)  =  r    (  j) 

new  ^  '  ' 

B)  c  (  i)  =  c    (  i) 

new  ^  ^ ' 


End  Repe  a  t , 


VII)   r(n+l)  =  r    (n+1) 
new  *"   ' 


End  Al gor  i  thm 
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Splitting  an  oversized  box  in  two  does  not  guarantee  that 
the  two  daughter  boxes  will  be  smaller  than  2Ar  .  If  the 
original  box  is  greater  than  4Ar,  the  the  two  new  boxes  will 
be  oversized  as  well.  Also,  eliminating  an  undersized  box 
by  dividing  its  volume  and  contents  between  its  neighbors 
could  cause  either  or  both  of  its  neighbors  to  become  over- 
sized where  they  had  not  been  before.  Thus  the  box  splitting 
routine  must  re-check  the  sizes  of  the  boxes  that  were  al- 
tered before  proceeding  further  in  the  array. 

2 . 8    TIME  INTERVALS 

2.8.1    Sedimentation  T r a n^X£ i  Time 

The  time  for  a  sedimentation  transfer,  At„,  is  usually 
chosen  so  that  each  sedimentation  transfer  will  create  a 
space  at  the  top  of  the  array  that  is  exactly  sufficient  for 
the  addition  of  one  box. 


Tjjg^d)  =  r(l)  +    Ai    =    r(l)exptS(l)a)2At  ]       (107) 


1     Frd)  +  Arl 

At   = -In (108) 

5(1)0)2   1^    r(l)   J 


S(l)  is  the  sedimentation  coefficient  at  the  first  bound- 
ary and  is  usually  taken  as  the  sedimentation  coefficient  of 
the  monomer  at  infinite  dilution  (the  concentration  at  the 


^Tr:-- 
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first  boundary  is  usually  zero).  When  two  or  more  constitu- 
ents are  present  the  sedimentation  coefficient  of  the  fast- 
est sedimenting  monomer  is  used  to  calculate  At_. 

2.8.1.1    Diffusion  Transfer  Time 

Complete  simulation  of  velocity  sedimentation  during  the 
interval  ^tg  includes  simulation  of  diffusion  for  an  equal 
time.  Thus  a  convenient  time  interval  for  diffusion  opera- 
tions, itp,  might  be  simply  Atg.  However,  Atg  is  usually 
too  long  for  a  single  diffusion  operation.  It  was  noted 
earlier  that  diffusion  times  must  be  short  because  the  simu- 
lation expressions  for  diffusion  do  not  account  for  gradient 
changes  during  the  diffusion  time  interval.  The  most  seri- 
ous error  resulting  from  a  long  diffusion  time  would  occur 
at  the  inital  sharp  solute  boundary.  There  the  steep  gradi- 
ent would  result  in  an  initial  rapid  diffusion  of  solute 
from  the  first  box  containing  solute  into  the  empty  box  im- 
mediately above  it.  In  a  real  cell  the  initial  diffusion 
would  rapidly  deplete  the  initially  steep  gradient  continu- 
ously decreasing  the  rate  of  diffusion  at  the  initial  sharp 
boundary  position.  In  the  discontinuous  simulation  opera- 
tion, however,  the  gradient  and  rate  of  diffusion  are  con- 
stant throughout  Atp.  if  At^is  too  long,  then,  an  inappro- 
priately large  amount  of  solute  would  be  transferred  in  the 
first  diffusion  operation  from  the  first  full  box  to  the 
empty  box  above  it.      ».*  i_..  Li--,-;.  '  '..J  .  ■   -  *>  t 


.  1  ■  .■.,■«.- 
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The  choice  of  a  maximum  value  for  At  is  based  on  the 
flow  of  solute  at  the  initial  sharp  boundary  in  a  rec tangu- 
li^    cell  with  equal  box  widths  (Ar)  (see  figure  9). 


Boundary:       j-1  j  j.l  j»2 

Box:  j-1  j  j.l 

Concentration:  O  O  C„ 


Figure  9:   Rectangular  Boxes,  Equal  Box  Widths 


The  mass  of  solute  that  will  pass  from  box  i+1  to  box  i 
in  At   is 


Am^  =  DAAtp 


i  +  1 


Ar 


=  DAAt^-- 
Ar 


(109) 


where   A   is   the   cross   sectional   area,   D   is   the   diffusion 
coefficient  and  (C,^^  -  C.)/Ar  is  the  gradient  at  boundary 


i+1  , 


The  mass  of  solute  that  will  flow  from  box  i  into  box  i+1 


in  At   is 


c .  -  c  •_, 

A",  =  DAAt-- —  =  0. 

Ar 
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(110) 


The  change  in  the  concentration  of  box  i,   the  first  box 
above  the  initial  sharp  boundary,  is  then 


Ac. 


Am^  -  Am^ 


AAr 


AAr 


(111) 


Substituting  the  right  side  of  (109)  for  Am   and  simplifying 
gives 


Ac  .  = 


"^^dCo 
'"Ar2 


(112) 


Ac  . 


aC_   where  a 


DAt, 


Ar-^ 


(113) 


a  is  the  fraction  of  solute  initially  in  box  i  +  1  that  is 
transfered  to  the  empty  box,  i,  in  the  first  diffusion  oper- 
ation. Obviously,  a  cannot  be  greater  than  0.5  since  if  it 
were,  then  the  sign  of  the  gradient  at  boundary  i+1  would  be 
reversed,  implying  that  some  of  the  solute  from  box  i  +  1  had 
travelled  into  box  i  against  a  positive  gradient.  An  arbi- 
trary maximum  value  for  a    has  thus  been  set  at  0.2. 


DAt, 


<  0.2 

Ar2  ~ 
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(114) 


0.2D 


^*D  < 


Ar  .2 
mm 


(115) 


max  ^'  *''*  largest  diffusion  coefficient  that  conld  occur 
in  the  system  -  usually  that  of  the  smallest  species  pres- 
ent, '^'nj^  is  the  smallest  allowable  box  width  (0.8  Ar,  ). 
^max  *^^  ^^min  *^®  used  to  ensure  that  a  will  not  be  greater 
than  0.2  under  any  circumstances  at  any  time  or  place  in  the 
cell  during  the  simulation. 

To  simulate  diffusion  for  a  time  equal  to  At-  it  is  nec- 
essary to  find  a  At   smaller  than  0 .2D     /Ar  -^  that  is  an 
D  max  '   min 

integral  factor  of  At_: 


-^'S  =  "D'^tD 


(116) 


The  simulated  transport  corresponding  to  a  total  time  of 
Atg  consists  of  n^  rounds  of  diffusion  for  At^  seconds  each 
followed  by  one  sedimentation  operation  corresponding  to  At 
seconds   of   sedimentation.    Atp  and   n^  can  be   computed  as 
f  ol  lows  : 
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0.2D 

max 


^'     ^*Dmax  =  --  — r  ("■" 

Ar  .  ^ 
min 


II  At-   II 


II)    n   =  11 II  +  1  (118) 

'itn     I  I 
Dmax ' < 


III) 


^^s 


'^^D  ^    -~-  <119' 


2 . 9    TI ME  AVERAGING 

In  two  constituent  systems  the  boundaries  of  the  two  con- 
stituent tracks  will  often  sediment  at  different  rates  so 
that  during  Atg  a  boundary  in  one  of  the  constituent  tracks 
will  move  through  a  range  of  concentrations  in  the  other 
constituent.  If,  as  is  usually  the  case,  the  local  sedimen- 
tation coefficient  of  each  constituent  depends  on  the  local 
concentration  of  the  other  constituent,  then  the  sedimenta- 
tion coefficients  at  the  boundaries  will  change  during  At 
Thus  new  boundary  positions  computed  from  the  initial  local 
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sedimentation  coefficients  will  be  systematically  in  error. 
This  error  can  be  largely  eliminated  by  computing  boundary 
sedimentation  coefficients  from  local  constituent  concentra- 
tions that  are  time-averaged  over  the  interval  At  The 
t ime- aver  aged  boundary  conditions  are  obtained  as  follows. 

I)  Initial  boundary  constituent  concentrations  are 
calculated  in  the  usual  way  (eq.  60-61)  and  are 
stored. 

II)  The  initial  conditions  in  I  are  used  to  compute 
sedimentation  coefficients  for  each  boundary. 

III)  The  sedimentation  coefficients  in  II  are  used  to 
make  a  "provisional  shift"  of  the  boundaries.  (The 
provisional  shift  provides  approximate  end-of-shift 
conditions) 

IV)  The  boundary  conditions  at  the  end  of  the  provi- 
sional shift  are  computed  and  stored. 

V)  Boundary  constituent  concentrations  from  I  and  IV 
are  averaged. 

VI)  The  time  averaged  boundary  conditions  in  V  are  used 
to  compute  time  averaged  sedimentation  coefficients 
at  each  boundary. 

VII)  The  boundaries  are  returned  to  their  initial  posi- 
tions and  new  boundary  positions  are  computed  using 
the  averaged  sedimentation  coefficients  in  VI. 


Chapter  III 
EXPERIMENTS  AND  SIMULATIONS 

3 . 1    DNCOOPERATIVE  AB   SYSTEMS 

We  begin  by  looking  at  uncooperative  AB^  systems  in  which 
the  binding  of  a  molecule  of  B  to  one  of  the  binding  sites 
on  A  is  independent  of  the  occupancy  of  the  other  sites  on 
A.  For  these  systems  the  mole-scale  step  wise  association 
constants,  K.  (where  i=l,2...n),  are  related  to  the  intrin- 
sic association  constant,  K    as  follows: 

[AB  ]  (n-i+1) 

K.  = =   K 

'    [B][AB._^]        ^    i 

where  n  is  the  number  of  binding  sites  (on  A  for  B)  and 

[AB.], 


E 


1  J 


^    [B] [AB.  ,] . 


the  association  constant  for  the  formation  of  a  particular 
species  of  AB.,  (AB^jj,  by  the  binding  of  a  molecule  of  B  to 
a  particular  binding  site,  j,  on  A. 

The  simulated  gradient  profile  (sehlieren  pattern;  (dc/ 
dr)  vs  r)  of  an  uncooperative  AB^  system  is  shown  in  figure 
10a.  In  this  particular  hypothetical  system  the  intrinsic 
association  constant,  K^ ,  is  1.0x10^  (fairly  tight),  the  mo- 
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lecular  weights  of  both  monomers  (A  and  B)  are  the  same 
(lOOKd),  and  the  frictional  ratios  (f/fo)  of  all  species  (A, 
B,  AB,  AB^)  are  equal.  The  initial  constituent  concentra- 
tions, C^^  and  Cgj,  are  both  5  mg/ml  or  5ilO~^M  (1:1  con- 
stituent mole  ratio) 

The  theoretical  gradient  in  figure  10a,  produced  by  the 
COXMIX  distorted-grid  simulation  program  (for  systems  under- 
going mixed  association  with  rapid  equilibration)  with  the 
SBABN  table  assembly  program  (for  AB  systems),'  represents 
the  situation  in  the  ul  tra cent r i f age  cell  after  2155.5  sec- 
onds of  velocity  sedimentation  at  60,000  rpm . ( 241 , 170ig  at 
r^  =  6  .00cm. )  . 

The  migrating  boundary  of  this  system  is  clearly  bimo- 
dal  .  The  trailing  boundary  migrates  at  the  sedimentation 
rate  of  free  monomer.  The  position  of  the  trailing  boundary 
(Rj)  is  the  same  as  that  predicted  for  free  monomer  with  no 
hydrodynamic  dependence.  (The  predicted  positions  of  indi- 
vidual species  with  minimum  and  maximum  hydrodynamic  depen- 
dence are  given  in  table  2).  The  individual  constituent 
gradients  (see  fig.  10)  indicate  that  the  trailing  boundary 
consists  of  excess  A  monomer  -  in  agreement  with  its  sedi- 
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The  simulation  and  table  assembly  programs  used  thro 
this  work  can  be  found  with  Dr.  David  J.  Cox,  Depa 
of  Biochemistry,  Kansas  State  University. 

In  general,  AB^  systems  with  n>l  give  bimodal  bound 
provided  the  association  is  not  too  weak  (see  fig 
Weakly  associating  AB  systems  with  n>l  have 
boundaries.  AB  (n  =  l)  iioundaries  are  not  bimodal 
11),  provided  that  the  initial  molar  ratio  of  con 
ent  s  i  s  1 : 1 . 


ughout 
r  tment 


a  r  1  e  s  , 
.  10)  . 
skewed 
(fig. 
s  t  i  t  u- 
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Figure  10:   Gradient  Profiles  of  Dncoope r a t i v e  AB   Systems 


lOOKd 


Wg  =  lOOKd 


''AT  =5.0  mg/ml 

[A]  =  5xlO~5 

[B]  =  5xlO"5 

[B]/[A]  =  1.0  " 

Frictional  ratios  (f/f^)  for  all  species  =  1.1 

All  self  and  cross  hy drodynam i c  constants  (k)  =  0.01 

Species  transport  coefficients  are  given  in  table  2. 

Rotor  speed  =  60,000  rpm. 

Initial  sharp  boundary  position  (r  )  =  fi.o  cm. 

At  =  2155  .51  sec. 


AB, 


•>•  ABj,  Kj  =  10- 


AB, 


100 


^      ::>•' 
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mentation  rate.  However,  the  concentration  of  constituents 
at  the  trailing  boundary  are  not  negligible  (particularly 
not  at  earlier  times),  and,  therefore,  if  the  trailing 
boundary  were  migrating  as  free  monomer  throughout  the  simu- 
lation its  final  position  would  be  expected  to  reflect  sig- 
nificant' hydrodynamic  effects. 

That  the  trailing  boundary  moves  at  an  average  rate 
greater  than  that  of  free  monomer  with  a  plausible  concen- 
tration dependence  is  probably  due  to  participation  of  ma- 
terial in  the  trailing  boundary  (excess  A)  in  aggregated 
species  (AB,  AB^)  at  earlier  times  in  the  experiment. 

The  position  of  the  leading  boundary,  R^^,  corresponds  to 
an  average  sedimentation  rate  greater  than  that  of  AB.  (at 
BT  o''  ^'^*  considerably  less  than  that  of  AB-  (at 


AT,  o 


ind  C 


'^AT,o  *°'^  ^BT,o)  .  ie. 

*minAB  =  «-54  '■<"' ■ '><\=6  .5  9<<R^.^^^^  =  6  .12  . 

Gradient  profiles  for  uncooperative  AB.   and  AB . 

3         '' 

are  also  shown  in  figure  10.  As  with  the  AB,  system  the 
profiles  of  the  AB^  and  AB^  systems  are  bimodal,  the  indi- 
vidual constituent  gradients  show  that  the  trailing  bounda- 
ries consist  almost  entirely  of  constituent  A  and  the  trail- 
ing boundaries  migrate  at  average  rates  somewhat  faster  than 
would  be  expected  for  concentration  dependent  sedimentation 
of  free  monomer.  Although  \IAB^)  >  Rl(AB3)  >  Rl(AB2)  lead- 
ing boundary  positions  foir  both  of  these  systems  still  fall 

be  tween  R    ,„  ._ j  r 

minAB  *°''  "minAB. 

....        .  -  I    .'"■•:    '-'.'. 


>.    'i". 
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Figure  11:   Gradient  Profiles  of  AB  Systems  of  Different 
S  tr ength  s 


W^  =  lOOKd 

"b  =  lOOKd 

'-j^j  =  5  .0  mg/ml  ■  " 

''BT  =  5  .0  mg/ml 

[A]  =  5x10"^  ,   ■ 

[B]  =  5xlO"5 
[B]/[A]  =1.0 

Frictional  ratios  (f/f^)  for  all  species  =  1.1 

All  self  and  cross  hy  dr  cidy  nam  i  c  constants  (k)  =  0.01 

Species  transport  coefficients  are  given  in  table  2. 

Rotor  speed  =  60,000  rpm. 

Initial  sharp  boundary  position  (r  )  =  6.0  cm. 

At  =  2155.51  sec 


-'SBF 


eoT 
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TABLE  2 


Species  Transport  Coefficients  and  Positions  (at  2155.51sec) 

'b= 


for  Systems  with  W^=Wg=100Kd . 


Spe  cie  s 

S 

0 

D 

* 

« 
m  in 

(il0^3  sec 

) 

dio''  CB^/sec) 

(cm) 

(cm) 

A 

7.0446 

6.3558 

6.371 

6.336 

B 

7.0446 

6.3558 

6.371 

6.336 

AB 

11  .183 

5.0446 

6.599 

6.542 

*«2 

14.653 

4  .4069 

6.797 

6.720 

*»3 

17.752 

4.0041 

6.978 

6.883 

AB, 

20.598 

3  .7168 

7.149 

7.036 

V  = 

0 

.73 

cm3/»8 

P  =■ 

1 

.0 

g/cm3 

11  = 

0 

.01 

po  i  se 

m  1  n 
wher  e 


rjjeip(S^u^At) 
r„exp(S„u^At/(l  +  k 


ave'^AT 


+  Cg.,.))) 


(I)  =  2000n  radians/sec 
At  =  2155.51  sec 


BT 


.  =  0.01  ml/mg 
=  5.0  mg  /ml 
=  5.0  mg/ml 


)•  A 


lOS 

3.1.1    Effect  of  Kj  on  Boundary  Shafie 

As  one  would  expect,  K^  is  an  important  determinant  of 
boundary  shape.  Gradient  profiles  of  several  uncooperative 
AB^  systems  with  different  K- s  are  shown  in  figure  12. 

With  weakly  associating  AB^  systems  (Kj<2.5xl0^)  bimodal 
behavior  is  not  observed.  At  K  =ixio^  nearly  all  of  the 
material  in  the  initial  plateau  is  monomer  (see  figure  13). 
The  resulting  boundary  (figure  12a)  is  nearly  symmetrical 
and  sediments  at  an  average  rate  just  slightly  greater  than 
that  of  free  monomer  (at  C^^^  and  C^j^)  .  As  Kj  is  in- 
creased, up  to  ca.  2.5x10^,  the  boundaries  become  more  no- 
ticeably skewed  (fig  12b)  and  an  obvious  trailing  shoulder 
develops  in  the  constituent  A  gradient,  paralleling  an  in- 
crease in  C^g  and  C^g  and  the  appearance  of  a  relative  ex- 
cess of  free  A  monomer  (ie.  relative  to  free  B  monomer)  in 
the  initial  plateau  (fig.  13).  A  trailing  shoulder  appears 
in  the  total  boundary  at  Kj=5il0*  and  at  K  =10^  (fig.  12c) 
the  boundary  is  fully  bimodal. 

The  position  of  the  leading  (reaction)  boundary^^  in- 
creases steadily  as  K^  is  increased  from  10^  to  about  lo'' . 
A  further  increase  in  Kj.  has  no  effect  on  the  shape  of  the 
boundary.  The  increase  in  reaction  boundary  sedimentation 
rate  with  Kj  js  related  to  an  increase  in  associated  spec- 
ies concentrations  in  the  initial  mixture.  Beyond  K-=lo' 
the   initial   species   concentrations   do  not   change   signifi- 
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The  unlmodal  boundaries  of  weak  systems  were  considered 
to  be  the  reaction  boundaries  of  these  systems. 


10« 


Figure  12:   Gradient  Profiles  Uncooperative  AB   Systems  of 
Different  Strengths 
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Figure  13:   Species  Concentrations  and  Average  Sedimentation 
Coefficients  for  Dncoopera ti ve  AB,  Systems 
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oantly.  The  initial  average  constituent  sedimentation  coef- 
ficients also  approach  maximum  values  at  E  _io'  (fig.  13b) 

With  systems  that  are  not  very  tight  (K-.<io^)  the  size  of 
the  leading  boundary  also  depends  on  Kj .  in  the  ABj  series 
(fig.  12)  leading  boundary  height  (H^^)  decreases  as  Kj  is 
increased  up  to  about  2.5x10^  where  E^  passes  through  a  min- 
imum (see  fig  18b).  E^  increases  slightly  with  Kj-  from  oa. 
2.5x10   to  ca.  10*  and  is  constant  at  higher  values. 

The  height  and  position  of  the  trailing  boundary,  when 
there  is  a  trailing  boundary  (ie.  when  K-.>5xl0^),  are  con- 
stant . 

The  changes  in  gradient  profiles  of  AB^  and  AB^  systems 
with  K  (figures  14  and  15)  follow  the  same  general  pattern 
as  in  the  AB   systems: 

I)  ABj  and  AB.  boundaries  are  bimodal  at  K,>ca. 
2.5x10  and  ca.  1.5x10  ,  respectively,  and  skewed 
at  1  ower  val  ue  s  . 

II)  Reaction  boundary  position,  R  increases  (sigmoi- 
dally  with  logjgKj  -  see  fig.  18b)  to  a  maximum  at 
Kj.~10  and  remains  constant  at  higher  values.  As 
with  AB^  systems  this  is  related  to  increases  in 
the  concentrations  of  associated  species  and  in  the 
average  constituent  sedimentation  coefficients 
(fig.  16  and  17) 

III)  Leading  boundary  size  decreases  as  K-  is  increased. 
However,   in  the  AB3   and  AB4   systems,   H^  does  not 
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pass  through  a  minimum  as  it  does  with  AB.  systems. 
Instead  it  decreases  steadily  to  K   =  ca.   10^  and 
remains  constant  at  higher  values. 
IV)    Trailing  boundary  height,  H.^,  and  position,  R^. 
constant. 
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Figure  14:   Gradient  Profiles  of  Uncooperative  AB 
of  Different  Strengths 
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Fi  gur e  1 5 : 


Gradient  Profiles  Uncooperative  AB.  Systems  of 
Different  Strengths 
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[B]/[A]  =  1 .0 

frictional  ratios  (f/f^)  for  all  species  =  1.1 

all  self  and  cross  hy drody nam i c  constants  (k)  =  0.01 

Species  transport  coefficients  are  given  in  table  2. 

rotor  speed  =  60,000  rpm. 

initial  sharp  boundary  position  (r  )  =  g.o  cm. 

At  =  2155  .51  sec 


a.  l.OilO^ 

b.  5.0x10^ 

c.  5.0x10'' 

d.  1.0x10* 


115 


"id 


116 


Figure  16:   Species  Concentrations  and  Average  Sedimentation 
Coefficients  for  Uncooperative  AB,  Systems 
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Species  transport  coefficients  are  given  in  table  2, 

a.  Species  concentrations  vs.  log   Kt  for 
uncooperative  AB^  systems 

b.  Average  constituent  sedimentation  coefficients 
vs.  ^"SjgKj  for  uncooperative  AB,  systems 
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Figure  17:   Species  Concentrations  and  Average  Sed i ment a t i or 
Coefficients  for  Uncooperative  AB.  Systems 
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3,1,2    Comj2_a  r^^ o n  o^  Un^ o_o£^ r^l v^  ABn  Bo un d a r^  Sh a£^^ 


Since  sediment ing  boundary  shapes  of  uncooperative  AB 
systems  vary  systematically  with  Kj  it  should  not  be  diffi- 
cult  to  distinguish  between,  say,  two  (uncooperative)  AB- 
systems  of  different  strengths  -  as  long  as  every th ing  else 
about  the  system  pertinent  to  its  transport  behavior  is 
known.  That  is,  if  K^  is  the  only  unknown  (which  it  rarely 
is)  then  modelling  will  give  an  unambiguous  solution.  A 
comparison  of  boundary  shapes  of  the  three  families  of  un- 
cooperative AB^  systems,  AB2 ,  AB3 ,  and  AB^ ,  suggests  that  it 
might  also  be  possible  to  obtain  an  unambiguous  solution 
when  both  s  t  o  i  ch  i  ome  try  and  K_  are  unknown.  First,  the 
three  st 0 i ch i ome tr i e s  give  rise  to  leading  boundaries  with 
distinctly  different  heights  at  values  of  K_.  greater  than 
ca,  10  (figure  18-a).  That  is,  the  leading  boundary  height 
of  a  system  with  one  of  these  three  s  t  o  i  chi  ome  tr  i  e  s  and 
Kj>ca.lO  is  different  than  the  leading  boundary  height  of 
any  system  with  one  of  the  other  two  s  t  o  i  ch  i  ome  tr  i  e  s  and 
K  >io  .  Further,  since  systems  with  any  of  these  stoichiom- 
etries  and  K<i.5ii0^  are  not  bimodal  (AB_  systems  are  bimo- 
dal  at  Kj>5xl0'^,  AB^  at  Kj>2.5xl0'*»  and  AB^  at  Kj>1.5xl0'** 
it  follows  that  systems  giving  rise  to  bimodal  boundaries 
have  K_ s  greater  than  10  and  therefore  have  "unique"  bound- 
ary shapes.  Thus,  for  example,  the  bimodal  boundary  of  an 
uncooperative  AB^  system  would  not  look  like  the  boundary  of 
any  AB^  or  AB^  system. 
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In  addition,  since  trailing  boundary  heights  differ  among 
systems  with  different  s t o i ch i ome tr i e s ,  but  are  constant,  or 
very  nearly  constant,  among  systems  with  the  same  stoichiom- 
etry  (H^(AB2)=  11  .0-11  .  7  (  mg  /  cm'' )  .  H^(AB3)=  13.0-13.8, 
Hj(AB^)=  14.2-14.8),  the  trailing  boundary  height  of  a  bimo- 
dal  AB  gradient  profile  could  allow  one  to  distinguish  be- 
tween various  stoichiometrics. 

For  weaker  systems  (not  bimodal)  boundary  size  by  itself 
cannot  be  used  to  distinguish  between  stoichiometrics  (eg. 
from  fig.  18a,  an  AB^  system  with  Kj  =  2.5xl0^.  an  AB^  system 
with  K  =4x10  ,  and  an  AB,  system  with  Kj=8xl0^  all  give  rise 
to  boundaries  of  identical  height.).  Similarly,  reaction 
boundary  position  alone  is  not  diagnostic  for  s t o i ch i ome t ry 
and  Kj .  Combinations  of  H,  and  R,,  however,  appear  to  be 
unique  for  AB^  systems  that  are  not  very  weak.  A  plot  of 
Hj^  vs.  R,  for  AB,  ,  AB,  ,  and  AB.  systems  shows  that  systems 
with  leading  boundary  heights  less  than  ca.  67(mg/cn'') 
(corresponding  to  AB^  systems  with  Kj>lilO^,  AB^  systems 
with  Kj.>i.5il0^,  and  AB^  systems  with  Kj>lxlO^'  have  unique 
combinations  of  boundary  height  and  position,  meaning  that 
they  have  unique  boundary  shapes  (fig  18-c).  Thus,  for  ex- 
ample, an  AB^  system  with  K-j>2xl0^  will  have  a  boundary 
shape  unlike  that  of  any  AB,  or  AB.  system. 

The  superimposed  gradients  of  AB^  and  AB.  (uncoopera- 
tive) systems  in  figure  19  clearly  demonstrate  the  differ- 
ences between  these  two  families  of  boundary  profiles. 


.   .1,  ' 
-.  .   A, 
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Figure  18:   Leading  Boundary  Heights  (H  )  and  Positions  ( R, ) 
of  Uncooperative  AB   Gradient  Profiles 
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Figure  19:   Comparison  of  Uncooperative  AB   and  AB^  Gradient 
Profiles  * 
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While  it  may  be  possible  to  distinguish  among  AB.  ,  AB,  , 
and  AB^  systems  the  apparent  tendency  of  the  H_  vs.  R, 
curves  to  lie  closer  together  for  higher  s t o i ch iome tr i e s 
(ABj  and  AB^)  than  they  do  for  lower  stoichiometrics  (AB, 
and  AB^)  suggests  that  the  boundary  shapes  of  uncooperative 
systems  of  higher  stoichiometrics  would  be  more  similar  to 
one  another  than  the  AB^  ,  ABj  ,  and  AB^  systems  are  to  one 
another  and  that  modeling  of  velocity  sedimentation  would 
probably  give  more  ambiguous  results  for  such  systems. 


3.1.3    Effect  of  Constituent  Mol^  Ratio 

All  of  the  gradient  profiles  discussed  so  far  have  been  for 
systems  at  equal  initial  constituent  concentrations 
(5xlO"^M)  or  a  constituent  mole  ratio  (MR=[B]/[A])  of  1.0. 
In  every  case  where  there  was  a  trailing  boundary  or  shoul- 
der the  individual  constituent  gradient  profiles  demonstrat- 
ed that  it  was  due  to  trailing  A  monomer,  probably  reflect- 
ing a  relative  excess  of  constituent  A  in  the  initial 
mixture.  Thus  one  might  expect  that  if  the  initial  concen- 
tration of  A  were  decreased  relative  to  that  of  B,  while 
keeping  the  total  constituent  concentration  constant,  the 
excess  in  A,  and  the  trailing  boundary  or  shoulder,  would  be 
reduced.  A  series  of  profiles  for  the  uncooperative  AB, 
system  with  Kj=io^  at  various  mole  ratios  (fig.  20)  shows 
that  this  is  exactly  what  happens.  As  MR  is  increased,  the 
trailing  boundary  of  A  shrinks,  closes  in  on  the  edge  of  the 
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leading  boundary  until  it  becomes  a  shoulder  at  the  trailing 

edge  of  the  leading  boundary  and  then  disappears  altogether. 

As  MR  is  increased  beyond  this  point  a  trailing  boundary  of 

B  develops  (first  as  a  shoulder  then  a  boundary) 

The  mole  ratio  at  which  there  is  no  excess  of  either  A  or 

B,   the   equivalence   mole   ratio   (MR   ),   depends   on   the 

eq 

strength  (Kj)  and  s t o i chi ome t ry  of  the  association.  In  gen- 
eral, tightly  associating  systems  have  equivalence  points  at 

higher  mole  ratios   than  weakly   associating   systems.    MR 

eq 

for  a  number  of  AB^  and  AB^  systems  were  determined  by  exam- 
ining  simulated  profiles   for   each   system   at   several   mole 

ratios.   The  equivalence  point  (MR   )  was  taken  as  the  mole 

eq 

ratio  at  which  the  trailing  edges  of  the  constituent  A  and  B 
boundaries  coincided.  A  plot  of  MR  vs.  log.nKr  is  shown 
in  figure  21.  Like  S^  and  Rj^  at  MR=1  ,  ME^  discriminates 
well  among  systems  with  the  same  s t o i ch i ome t ry  and  different 
values  of  K^. .  Moreover,  whereas  H,  and  R,  (at  MR=1)  are  es- 
sentially constant  among  tightly  associating  systems  (K-.>ca. 
5x10  ),  ^^gQ  continues  to  increase  beyond  Ky  =  10  and  may 
therefore  be  most  useful  for  modeling  relatively  tight  sys- 
tems. As  with  E^  and  Ej^  however,  MR^  can  have  the  same 
value  for  two  systems,  or  possibly  several  systems,  with 
different  stoichiometrics.  Thus  MR  by  itself  is  not  an 
unambiguous   indicator   of   both   s t o i ch iome try   and  K^ ,    How- 


ever,  combinations  of  MR 
\    (at  MR=1) ,  may  be. 


g   and  Ej     (at  MR=1)  or  of  MR^   and 


A  ^  i    . 


;V!    .. 
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Figure  20:   Gradient  Profiles  Uncooperative  AB.  Systems 
Different  Constituent  Hole  Ratios 
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Species  transport  coefficients  are  given  in  table  2. 
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Figure  21:   Equivalence  Mole  Ratio  vs.  log   Kt  for 
Uncooperative  AB^  and  AB4  Systems 
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Species  transport  coefficients  are  given  in  table  2. 

Rotor  speed  =  60,000  rpm. 

Initial  sharp  boundary  position  (r  )  =  6.0    cm. 

At  =  2155.51  sec 
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While  ME^^  could  be  helpful  in  some  instances  for  distin- 
guishing anong  various  model  systems,  there  are  some  prob- 
lems that  restrict  the  usefulness  of  this  type  of  analysis. 
First,  in  order  to  obtain  a  value  of  MR^^  for  a  real  system, 
it  would  be  necessary  to  run  the  system  in  the  ultracentri- 
fuge  at  several  mole  ratios  and  to  make  accurate  reproduc- 
tions or  tracings  of  the  schlieren  patterns  for  each  run. 
This  would  be  very  time  consuming.  Second,  even  if  the  num- 
ber of  possible  models  were  small,  many  costly  simulations 

would  be  required  to  determine  MR    for  each  of  them.   Fi- 

eq 

nally,  in  many  oases  it  will  not  be  possible  to  accurately 
determine  MR^  for  the  experimental  system.  The  MR  's  in 
fig.  21  were  determined  by  inspecting  the  simulated  constit- 


uent gradients  at  various  mole  ratios  with  MR    beine  taken 

eq      ® 

as  the  mole  ratio  at  which  the  trailing  edges  of  the  two 
constituent  gradients  coincided.  For  a  real  system  only  the 
total  gradient  profile  is  available  and  while  with  moderate- 
ly tight  and  tighter  systems  the  trailing  aspect  of  the  to- 
tal boundary  does  pass  through  a  detectable  minimum  at  MR 
(see  fig.  21)  this  is  not  the  case  for  weaker  systemi 
the  AB^  system  with  K.  of  10  ,  for  instance,  the  trailing 
edges  of  the  total  boundaries  are  essentially  identical  be- 
tween mole  ratios  of  1.0  and  1.5,  despite  the  obvious  dif- 
ferences, over  the  same  range  of  mole  ratios,  in  the  indi- 
vidual constituent  gradients  (fig.   22).     '\  ■ 


eq 
For 
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3.1.4    Systems  wiH  Non-Identical  Monomers 

So  far  we  have  looked  only  at  systems  with  identical  monom- 
ers (both  lOOKd) . 

Gradient  profiles  of  uncooperative  AB  systems  (K^=lo'. 
MR=1.0)  with  monomers  of  different  molecular  weights  are 
shown  in  figure  23 

(a:W^=140Kd,Wg=60Kd,  b : W^=6  0Kd, Wg=140Kd ) ^ ^ . 

The  boundaries  of  both  systems  with  unequal  monomers, 
like  boundaries  of  systems  of  comparable  strength  (K  =io^) 
with  identical  monomers,  are  bimodal  at  a  constituent  mole 
ratio  of  1.0.  However,  while  the  trailing  boundaries  (at 
MR=1)  of  systems  with  identical  (lOOKd)  monomers  and  of  the 
systems  with  W^  (140K)  larger  than  W^  (60K)  consist  of  ex- 
cess constituent  A,  it  is  constituent  B  that  is  in  excess 
and  present  in  the  trailing  boundary  (at  MR=1)  of  the  system 
wj'li  W^(i40K)  larger  than  Wg  (60K). 

As  with  systems  with  identical  monomers  (W,  =  Wj,)  the 
trailing  boundaries  of  both  of  these  systems  (*a>Wi,  and 
*g>W^)  move  at  about  the  rates  predicted  for  pure  monomer 
(the  monomer  in  excess).  (Transport  coefficients  and  pre- 
dicted final  positions  of  the  pure  species  in  these  systems 
are  given   in  table  3   and  4.)    Thus,   in  the  W_>w,   system. 


1 2 

The   sedimentation  and  diffusion   transfer   times   and  the 

total  time  simulated  were  shorter  here  than  in  the  previ- 
ous simulations.  At^  is  computed  (by  the  COXMIX  simula- 
tion program)  as  the  time  it  takes  for  the  pure  monomer 
with  the  largest  sedimentation  coefficient  to  move 
through  a  distance  equal  to  the  average  box  width,  (see 
p.  90) 
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Figure    22:       Gradient    Profiles    Uncooperative    AB       Systems    at 
Different    Constituent    Mole    Ratios 


Kj    =    1.0x10* 
W^    =    lOOKd 
Wg    =    lOOKd 

^AT    ■*■     '-BT    ^    ^^  •"    "g/ral 

[B1.J,    +     Ibkl.A]^    =    l.OxlO"'* 

Frictional  ratios  (f/f  )  for  all  species  =  1.1 

All  self  and  cross  hy drody nam i c  constants  (k)  =  0.01 

Species  transport  coefficients  are  given  in  table  2. 

Rotor  speed  =  60,000  rpm. 

Initial  sharp  boundary  position  (r  )  =  6.0  cm. 

At  =  2155  .51  sec 

[B]/[A] 

a.      1.00 
h,  1.08 

c.      1.22 


»s 


■      -  A 
•        B 
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Figure    23:       Gradient    Profiles    of    Uncooperative    AB,     Systf 
with    Non-Identical    Monomers 


Kj    =    l.Oxlo' 

Frictional  ratios  U/i^)     for  all  species  =  1.1 

All  self  and  cross  hydrodynamic  constants  (k)  =  O.Ol 

Rotor  speed  =  60,000  rpm. 

Initial  sharp  boundary  position  (r  )  =  g.o  cm. 

At  -  1722.50  sec 

a.  W^  =  140Kd 

Wg  =  60Kd  r 

*'AT  -    T  -^    mg/ml 
[A]  =  5.0x10"^  M 
Cgj  =3.0  mg/ml 


[B] 


5  .0x10 


[B]/[A]  =  1.0 

Species  transport  coefficients  are  given  in  table  3. 

"a  =  60Kd 

*B  =  140Kd  "     .  , 

^AT  °  ^ •"  "g/™! 
[A]  =  5.0x10"^  M 


BT 
[B] 


7  rO  mg  /ml 
5  .0x10"^  M 


[B]/[A]  =1.0 

Species  transport  coefficients  are  given  in  table  4. 
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A 
B 


"lb 
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in  excess. 


where  constituent  A 

^minA=6-19(«"'-'<RT=6-20<R„^j.A=6.21, 


and  in  the  *^>Wg  system,  where  constituent  B  is  in  excess, 
^ninB=RmaxB=«-19 

^°^  *A^*B  systems  the  shape  of  the  gradient  profile  seems 
to  depend  on  Kj.  in  about  the  same  way  that  it  does  with 
*j^=Wg  systems  (fig.  24):  Increasing  K-  decreases  the  height 
and  increases  the  average  sedimentation  rate  of  the  leading 
boundary.  However,  the  trailing  boundary  heights  of  the 
*^>Wg  systems  decrease  with  K  (fig.  24),  whereas  among 
*^=Wg  ( AB2  ,  unco  ope  r  a  t  i  V  e  )  systems  those  strong  enough  to 
give  a  bimodal  profile  under  the  conditions  simulated  had 
trailing  boundaries  of  the  same  size. 

A  possible  explanation  for  this  difference  might  be  that 
the  sedimentation  coefficients  of  the  monomer  in  excess  and 
the  aggregate  species  are  more  different  in  the  W,>w„  sys- 
tems (S^g/Sg=2 .232)  than  in  the  Wa=*B  systems 
'^AB'' ^A^-^ -^  *^'  so  that  W^>Wg  systems  give  rise  to  bimodal 
boundaries  at  lower  values  of  K  where  the  amount  of  excess 
trailing  monomer  varies  more  strongly  with  K. . 

The  dependence  of  boundary  shape  on  mole  ratio  is,  in 
general,  qualitatively  the  same  for  W^>Wg  and  yi^>Vf^  systems 
as  It  IS  for  ^^=Wg  systems.  As  MR  is  increased  from  values 
where  A  is  in  excess  the  trailing  boundary  or  shoulder  is 
depleted   of   excess   A   up   to   a   point   (equivalence)   beyond 
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TABLE  3 

Species  Transport  Coefficients  and  Positions  for  Systems 
with  W^=140Kd  >  W3=60Kd 


Species 


(il0^3  sec)      dio'  cm^/sec)  (cm)     (cm) 


A 

8.8161 

5  .6815 

6.371 

6.336 

B 

5.0114 

7.5356 

6.208 

6  .189 

AB 

11  .183 

5.0446 

6.474 

6.429 

^^2 

13  .320 

4.6222 

6.569 

6.515 

AB3 

15  .298 

4.3131 

6.658 

6.595 

AB, 

17  .155 

4 .0729 

6.742 

6.671 

V    =    0. 

.73 

cm3/»g 

p  =  1  .0  g/cm-* 
n  =  0.01  poise 


mi  n 

wher  e 


'o"P<S„„^At) 

r„exp(S„,.^At/(l  +  '^^^^(Cf^j    +  Cg^)  )  ) 


(i>  =  2000n  radians/sec 
At  =  1722.50  sec 


BT 


,  =  0.01  ml/mg 
=  7.0  mg/ml 
=  3.0  mg/ml 
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TABLE  4 

Species  Transport  Coefficients  and  Positions  for  Systems 
with  W3  =  140Kd  >  Wj^=60Kd 


Spe  cie  s 


(ilO^^  sec)     dio''  cm^/sec)    (cm) 


(  cm) 


A 

5  .0114 

B 

8.8161 

AB 

11  .183 

AB^ 

15  .929 

*«3 

20.046 

AB. 

23  .775 

7 

.5356 

6 

.208 

6 

.189 

5 

.6815 

6 

.371 

6 

.336 

5 

.0446 

6 

.474 

6 

.429 

4 

.2268 

6 

.686 

6 

620 

3 

7678 

6 

.876 

6 

7  92 

3 

.4597 

7 

.053 

6 

950 

V  =  0.73  cm^/'-S 
p  =1.0  g/cm^ 
r\    =  0.01  poise 


m  1  n 
wher  e 


r^exp  (S^u-'At) 

r„exp(S„,o2At/(l  +  k_  (c.,.  +  C„x)  )  ) 


ave^'-AT 


^BT' 


u  =  2000n  radians/sec 
At  =  1722.50  sec 
''ave  °  °-°^    ml/mg 
*^AT  =3.0  mg/ml 
^BT  -    T  -^    mg/ml 


■VI 


■-,,   *•" 
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Figure  24:   Gradient  Profiles  of  Uncooperative  AB,  Systems 
of  Different  Strengths  with  Non-Identical 


Monomers  C^Mg) 


140Kd 


Wg  =  60Kd 


AT 


7.0  mg/ml 


[A]  =  5  .Olio  '  M 

''BT  ~  ^  •"  mg/ml  ,  ;  . 

[B]  =  5  .OilO"5  M  .     ;' 

[B]/[A]  =  1 .0 

Frictional  ratios  (f/t^)     for  all  species  =  1.1 

All  self  and  cross  hydrodynamic  constants  (k)  =  0.01 

Rotor  speed  =  60,000  rpm . 

Initial  sharp  boundary  position  (r  )  =  6.0  cm. 

At  =  1722.50  sec 

Species  transport  coefficients  are  given  in  table  3. 


a.  5.0x10^ 

b.  5.0x10^ 
0.    l.OxlO^ 


-^5 
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which  a  trailing  shoulder  or  boundary  of  B  develops  (see 
figs.  25  and  26).  There  are,  however,  some  differences 
worth  noting.  First,  the  shape  of  the  W  >w  (Kt=10^)  gradi- 
ent, which  has  a  trailing  boundary  of  excess  B  at  MR=1.0,  is 
not  bimodal  when  the  system  is  in  excess  A  but  is  instead 
skewed  to  the  left  with  the  constituent  A  gradient  (fig. 
25).  The  absence  of  a  trailing  boundary  for  this  system 
with  A  in  excess  is  probably  a  result  of  the  relatively 
small  difference  between  the  sedimentation  coefficients  of  A 
and  the  aggregates  ( S^^/ s^=l  .26  8)  .  In  the  W^=Wb  and  V^yV^ 
systems  (with  Kj.  =  io^).  which  have  bimodal  boundaries  when  A 
is  in  excess  (as  well  as  when  B  is  in  excess),  the  sedimen- 
tation coefficient  of  A  is,  by  comparison  with  the  W  >w 
system,  smaller  in  relation  to  the  sedimentation  coeffi- 
cients of   the  associated  species   (for  W„>w     S   /  <!  =7  ■>■*■} 

B   A'   AB'  A  '■•^"■• 

for  W^=Wg, S^b/S^=1 .587) .  When  constituent  B  is  in  excess 
both  the  W^>Wb  and  Wb>W^  systems  (Kj=10^)  give  bimodal  pro- 
files with  trailing  boundaries  at  positions  corresponding 
approximately  to  the  predicted  positions  of  pure  B,  even 
though  the  difference  between  S^  and  S^g  in  the  Wg)!*^  sys- 
tem is  the  same  as  between  S^  and  S^^g  in  the  Wa>Wb  system, 
which  does  not  have  a  distinct  monomer  boundary  when  A  is  in 
excess.  The  relatively  greater  difference  in  the  W  >w,  sys- 
tem between  Sg  and  S^g   ( S^g  /Sb  =  1.81)   than  between  S^  and 

^AB     ii  the  W^>Wb  system  ( S^g  /S4=1.51)  may  be  responsible 
2  2 
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Figure  25:   Gradient  Profiles  of  Uncooperative  AB   Systems 
with  W^>Wg  at  Different  Constituent  Mole  Ratios 


Kj  =  1  .Olio' 
"a  =  140Kd 
"g  =  60Kd 


''AT  ■*■  '^BT  "^  10.0  mg/nl 


Frictional  ratios  ( f / f  )  for  all  species  «  1.1 

All  self  and  cross  hydrodynam ic  constants  (k)  =  0.01 

Rotor  speed  =  60,000  rpm. 

Initial  sharp  boundary  position  (r  )  =  6.0  cm. 

At  =  1722.50  sec 

Species  transport  coefficients  are  given  in  table  3. 

[B]/[A] 


0.80 
1  .00 
1  .50 
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Figure  26:   Gradient  Profiles  of  Dncoopera t i ve  AB   Systems 
with  Wg>v(^  at  Different  Constitnent  Mole  Ratios 


Kj  =  1  .0x10^  M 

"a  =  «OKd 

"g  =  140Kd 

^AT  *  ^BT  ^  10.0  mg/ml 

Frictional  ratios  ( f / f  )  for  all  species  =  1.1 

All  self  and  cross  hydrodynamic  constants  (k)  =  0.01 

Rotor  speed  =  60,000  rpm . 

Initial  sharp  boundary  position  (r  )  =  g.o  cm. 

At  =  1722.50  sec 

Species  transport  coefficients  are  given  in  table  4, 

[B]/tA] 


1  .00 

1  .50 

2  .20 

3  .00 


rt/i' 
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*B^*A  systems  (with  Kj=10  )  are  much  different  than  that  of 
the     system    with     identical     monomers     (K-=io^):     MR  ([B]/[A]) 

for  the  W^>flg  system  (ca.  0.8)  is  smaller  than  that  of  the 
*j^=Wg  system  (ca.  1.5)  which  is  smaller  than  that  of  the 
"gM^    system     (ca.     2.1). 

3 . 2  COMPLETELY    COOPERATIVE    SYSTEMS 

So      far     we     have      considered     only      uncooperative     AB         svstems 

n    *' 

where   the   step-wise   association   constants,   K,  ,   K, .    .    K 
are  described  in  terms  of  K   in  such  a  way  that  the  intrin- 
sic binding  of  a  molecule  of  B  to  a  vacant  site  on  A  is  in- 
dependent of  the  occupancy  of  other  sites  on  A.   We  will  now 

consider   completely   cooperative   AB    svstems   in   which   the 

n    ' 

formation  of  AB^  occurs  in  a  single  concerted  step  (A  +  nB 

<=>   ^^n^    without   the   formation   of   intermediate   species. 

For  these  systems  the  mole-scale  association  constant,  K 

nM' 

is  defined  as  follows 

"^    [A][B]° 


Simulated  gradient  profiles   of   cooperative  AB,   and  AB . 

2  4 

systems  with  various  values  of  K  „  are  shown  in  figures  27 

nM  ° 

and  28.  As  with  uncooperative  systems,  the  boundaries  of 
cooperative  systems  are  bimodal  when  association  is  strong 
enough  and  have  shoulders  or  are  skewed  when  the  association 
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is  weaker.  Also  similar  to  uncooperative  systems  are  the 
position  of  the  trailing  boundary,  which  is  nearly  constant 
and  about  equal  to  the  position  predicted  for  pure  monomer, 
and  the  direction  of  change  in  the  leading  boundary  position 
with  increasing  strength  of  association.  Unlike  the  leading 
boundaries  of  uncooperative  AB^  systems,  however,  which  move 
at  average  rates  that  are  considerably  slower  than  those  of 
the  pure  AB^  species  no  matter  how  tight  the  association, 
the  leading  boundaries  of  cooperative  AB   systems  migrate  at 


nM 

ly  slower  than  those  of  pure  AB^.  R^^  for  cooperative  AB^ 
and  AB^  systems  studied  had  maximum  values  of  6.71  cm.  and 
6.98  cm.,  respectively  (see  fig.  29-a),  which  are  compara- 
ble  to   the   theoretical   positions   of   the  pure  AB,   and  AB 

2  4 

species  (6.72  cm.   and  7.03  cm.-  see  table  2).    These  large 

values  for  Rj^,   which  occur  in  tightly  associating  systems, 

suggest  that  the  leading  boundaries  of  these  systems  consist 

largely  of  AB^.   This  may  also  be  inferred  from  the  initial 

species  concentrations  and  average  constituent  sedimentation 

coefficients   (see  figures  30  and  31).    In  AB   systems  C,„ 

2.  AB^ 

increases  sigmoidally  with  logioK2M  ''•^■^  ■  «*  large  values  of 
^2M'  '*  t'"  dominant  species  in  the  mixture.  Also  at  large 
values  of  K2M  the  average  sedimentation  coefficient  of  con- 
stituent B  is  not  much  smaller  than  that  of  pure  AB     Thus, 

for  instance,  at  K,„=io^-'-: 
2  M 

^B  0  =  14  .20xl0"-'^^sec.  ,  and 
^AB   0  =  14  .65xlO~^^sec. 


ISO 


A  0  '°  these   systems  is  considerably  smaller  than  S.o 

AB^.O' 

as  it  should  be,  since  there  is  considerable  excess  A  in  the 
mixture.    Similarly,   in  cooperative  AB^  systems,  when  asso- 


ciation  is  tight  the  predominant  species  is  AB.  and  Sj 


preaches  S 


AB, 


At  K4„=10 


20  . 


-13, 


^B_  o=18.64il0~-^''sec.  ,  and 
^AB   „=20 .60ilO"^^sec. 

There  is  a  general  difference  between  the  series  of  coop- 
erative profiles  in  figures  27  and  28  and  the  series  of  un- 
cooperative profiles  in  figures  27-30  in  the  manner  in  which 
bimodality  develops  as  the  strength  of  association  is  in- 
creased. In  the  cooperative  series  this  involves  formation 
of  a  leading  boundary  from  the  leading  edge  of  the  original 
single  boundary  (which  becomes  the  trailing  boundary), 
whereas,  in  the  uncooperative  series,  a  trailing  boundary 
develops  at  the  trailing  edge  of  the  single  boundary  (which 
becomes  the  leading  boundary).  Another  related  difference 
is  that  in  the  cooperative  series  the  height  of  the  leading 
boundary  increases  with  K^^  and  the  height  of  the  trailing 
boundary  decreases  with  K^^j,  while,  in  the  uncooperative  se- 
ries, the  leading  boundary  height  decreases  with  K  and  the 
trailing  boundary  height  is  constant. 

Changes  in  R^ ,  Hj^,  and  H^  with  K^^j,  for  cooperative  AB^ 
and  AB4  systems  are  illustrated  in  figure  29.  Of  these 
quantities,  Rj^,  leading  boundary  position,  appears  to  be  the 
least  ambiguous  in  discriminating  between  AB,  and  AB,  (coop- 
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Figure  27:   Gradient  Profiles  of  Completely  Cooperative  AB 
Systems  of  Different  Strengths 


"^  =  lOOKd 
Wg  =  lOOKd 

Cg  J  =  5.0  mg/ml 

[A]  =  5xlO"5  M 

[B]  =  5ilO"5  M 
[B]/[A]  =  1.0      ■ 

Frictional  ratios  (f/f^)  for  all  species  =  1.1 

All  self  and  cross  hydrody nam i c  constants  (k)  =  0.01 

Species  transport  coefficients  are  given  in  table  2. 

Rotor  speed  =  60,000  rpm . 

Initial  sharp  boundary  position  (r  )  =  fi.o  cm. 

At  =  2155  .51  sec 


d. 


^2M    («''' 

2.5xl0'7 

1.0x10^ 

l.Oxlo' 

1  .OxlO^l 

-.       ■    1 
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Figure  28: 


Gradient  Profiles  of  Completely  Cooperative  AB 
Systems  of  Different  Strengths 


W^  =  lOOKd 

Wg  =  lOOKd 

^AT  =5.0  mg/ml 

Cgj=5.0mg/ml 

[A]  =  5ilO~5  M 

[B]  =  5ilO~5  M 
[B]/[A]  =1.0 

Frictional  ratios  (f/f^)  for  all  species  =  1.1 

All  self  and  cross  hy drodynam i o  constants  (k)  =  0.01 

Species  transport  coefficients  are  given  in  table  2. 

Rotor  speed  =  60,000  rpm. 

Initial  sharp  boundary  position  (r  )  =  6.0    cm. 

At  =  2155  .51  sec 


^4M  («"'> 


a.  LOxlO^* 

b.  l.OxlO^^ 


c.    1.0x10 


d.    1.0x10 


18 


19 


154 


■     -A 
•     -  B 


^"•"••M..a»i 


/•• 


■^tWr*y»wn«tTfJ:t^" ■■■'" 


155 


Figure  29:   Boundary  Heights  (H   and  Ej)    and  Leading 

Boundary  Positions  TR.  )  for  Cooperative  AB 
Gradient  Profiles  " 


"^  =  lOOKd 
Wg  =  lOOKd 
''AT  "  '  •"  "'g/'i'l 
''BT  ~    ^  •"  mg/inl 

[A]  =  5xlO~5  M 

[B]  =  5ilO"5  M 
[B]/[A]  =1.0 

Frictional  ratios  (f/f  )  for  all  species  =  1.1 

o  '^ 

All  self  and  cross  hydrodynam i c  constants  (k)  =  0.01 

Species  transport  coefficients  are  given  in  table  2. 

Rotor  speed  =  60,000  rpm. 

Initial  sharp  boundary  position  (r  )  =  6.0  cm. 

At  =  2155  .51  sec 


'■       L    ^^'     '°*10'^nM    ^°^    completely    cooperative 
ABj    and    AB^    profiles. 

"  •    H^    vs.     los^oK^jj    for     completely    cooperative 
^^2    and    AB^    profiles.  ,,      .      ,    p_ 

c.     Rj^    vs.     loSio^nji    for    completely    cooperative 
'^^2    and    AB^    profiles. 
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Figure  30:   Species  Concentrations  and  Average  Sedimentation 
Coefficients  for  Completely  Cooperative  AB, 

Systems  ^ 


W^  =  lOOKd 
Wg  =  lOOKd 
^AT  =5.0  mg/ml 
''BT  =  5  .0  mg/ml 

[A]  =  5xlO"5  M 

[B]  =  5xl0~5  M 
tB]/[A]  =1.0 

Friotional  ratios  (f/f  )  for  all  species  =  1.1 
Species  transport  coefficients  are  given  in  table  2, 


a.  Species  concentrations  vs.  logmKow  for 

10  2M 

completely  cooperative  AB   systems 

b.  Average  constituent  sedimentation  coefficients 
vs.  ^°8^qK2jj  for  completely  cooperative 


AB 


2  sy  s  t  em  s 
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Co(mg/ml 


0 


So(x10"sec) 
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(a 
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Figure  31:   Species  Concentrations  and  Average  Sedimentation 
Coefficients  for  Completely  Cooperative  AB 
Systems  ^ 


W^  =  lOOKd  '■ 

"b  =  lOOKd 

''AT  =5.0  mg/ml 

'-gX'  =  5  .0  mg/nl  -_ 

[A]  =  5xlO"5  M 

[B]  =  5xlO~5  M 

[B]/[A]  =1.0  -  ; 

Frictional  ratios  (f/f  )  for  all  species  =  1.1 
Species  transport  coefficients  are  given  in  table  2, 


a.  Species  concentrations  vs.  log  -k    fnr 

10  4M 

completely  cooperative  AB   systems 

b.  Average  constituent  sedimentation  coefficients 
vs.  ^^S^qK^jj  for  completely  cooperative 


AB 


.  sy  s  t  em  s 
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Co(mg/ml) 


■  V  s 


So(x10"sec) 


■•..    V 


1«1 

erative)   systems.    At  all  values  of  K^„  where  R^  could  be 

measured,  the  leading  boundaries  of  the  AB^  systems  migrated 

more  rapidly  than  the  leading  boundaries  of  any  of  the  AB 

2 

systems.  AB^  systems,  however,  which  were  not  included  in 
this  study,  would  most  likely  have  values  of  R  overlapping 
those  of  ABj  and  AB^  systems. 

Plots  of  Hj^  vs.  R^  (figure  32-a)  demonstrate  a  signifi- 
cant lack  of  similarity  (with  respect  to  H^^  and  R^)  between 
the  shapes  of  all  cooperative  AB^  and  AB^  profiles  (for  sys- 
tems with  leading  boundaries)  as  well  as  between  the  shapes 
of  cooperative  and  uncooperative  systems.  Some  uncoopera- 
tive AB^  (U-4)  and  cooperative  AB,  (C-2)  profiles  have  simi- 
lar or  identical  leading  boundary  heights  and  positions 
(Hj^~39  mg/cm'*  and  Rj^=cti.  6.62  cm.).  These  U-4  and  C-2  sys- 
tems,   however,    have   very   different    trailing   boundary 

4) 
4 


heights.    The  U-4  systems   in  this  range   (K  ~4.0xl0'*)   give 


trailing  boundaries  with  heights  of  ca.  14.2  mg/cm*  while 
the  trailing  boundaries  of  the  C-2  systems  in  this  range 
'^2M"^  •^^■'^''^  have  heights  of  ca.  22.0  mg/cm*  (see  figures 
32-b  and  33).  Likewise,  uncooperative  AB  (U-2)  and  cooper- 
ative ABj  (C-2)  systems  with  similar  leading  boundary 
heights  and  positions  (E^~S4  mg/cm''  and  Rj^~6.66  cm.)  have 
different  trailing  boundary  heights  (U-2:  H  =11.7  mg/cm'', 
C-2:  H^.20.8  mg/cm'') 

Plots  of  H^  vs.  Hj^  and  Hj  vs.  R^  further  demonstrate  the 
dissimilarity  between   the   boundary   shapes   of   C-2   and   C-4 
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Figure  32:   Comparison  of  Boundary  Dimensions  for 

Uncooperative  and  Cooperative  AB   Systems 
"a  =  lOOKd 

Wg  =  lOOKd 
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=  5.0  mg/ml 


[A]  =  5x10 
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M 


[B]/[A]  =  1.0 

Frictional  ratios  (f/f  )  for  all  species  =  1.1 

o  *^ 

All  self  and  cross  hydrodynamic  constants  (k)  =  0.01 
Species  transport  coefficients  are  given  in  table  2. 
Rotor  speed  =  60,000  rpm.  »• 

Initial  sharp  boundary  position  (r  )  =  6,0  cm. 
At  =  2155  .51  sec 
Solid  Lines:  Cooperative  Systems 
Dashed  Lines:  Uncooperative  Systems 

a.  H,  vs.  R,  for  completely  cooperative 

and  uncooperative  AB,  and  AB,  profiles 
b ■  Hj  vs.  R,  for  completely  cooperative 

and  uncooperative  AB,  and  AB.  profiles 
"  •     H^  vs.  H,  for  completely  cooperative 

and  uncooperative  AB,  and  AB.  profiles      ' 

cooperative  AB^  and  AB^  profiles 
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Figure  33:   Gradient  Profiles  of  an  Uncooperative  AB   and 
Completely  Cooperative  AB^  System  with  Similar 
Leading  Boundaries 


Wj^  =  lOOKd 
Wg  =  lOOKd 
Cj^j  =  5.0  mg/ml 
Cgj  =5.0  mg/ml 

[A]  =  5x10"^  M 

[B]  =  5il0"5  M 
[B]/[A]  =1.0 

Frictional  ratios  (f/f  )  for  all  species  =  1.1 

0  ^ 

All  self  and  cross  hydrodynami c  constants  (k)  =  0.01 

Species  transport  coefficients  are  given  in  table  2, 

Rotor  speed  =  60,000  rpm. 

Initial  sharp  boundary  position  (r  )  =  6,0  cm. 

At  =  2155  .51  sec 
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systems  as  well  as  between  cooperative  and  uncooperative 
systems  (fig.  32-b  and  o). 

Gradient  profiles  of  C-2  and  C-4  systems  that  are  not 
very  weak  but  that  are  not  strong  enough  to  give  bimodal 
boundaries  are  also  dissimilar  as  seen  in  figure  34  in  which 
a  series  of  C-4  gradient  profiles  is  superimposed  onto  a  se- 
ries of  C-2  profiles.  The  most  obvious  difference  is  that 
the  C-2  boundaries  have  higher,  shorter  leading  shoulders 
than  C-4  boundaries. 

Similar  sets  of  superimposed  boundaries  for  U-2  and  C-2 
systems  (figure  35)  demonstrate  non- s imil ar i ty  between 
non-bimodal  0-2  and  C-2  boundaries.  n-2  and  C-2  boundaries 
that  are  skewed  are  skewed  in  opposite  directions.  U-2  sys- 
tems are  skewed  to  the  right  (downward  or  outward)  and  C-2 
boundaries  are  skewed  to  the  left  (upward  or  inward).  U-2 
and  C-2  boundaries  that  have  shoulders  have  them  on  opposite 
sides  of  the  main  boundary.  U-2  boundaries  have  trailing 
shoulders  and  C-2  boundaries  have  leading  shoulders.  U-4 
and  C-4  boundaries  differ  similarly  (figure  36). 
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Figure  34:   Comparison  of  Completely  Cooperative  AB,  and  AB^ 
Gradient  Profiles  * 
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Species  transport  coefficients  are  given  in  table  2. 
Rotor  speed  =  60,000  rpm . 
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Figure  35:   Comparison  of  Completely  Cooperative  and 
Unco 
"a  =  lOOKd 


Uncooperative  AB^  Gradient  Profiles 
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Figure  36:   Comparison  of  Completely  Cooperative  and 
Uncooperative  AB   Gradient  Profiles 
W^  =  lOOKd 

"b  =  lOOKd 
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Frictional  ratios  (f/f^)  for  all  species  =  1.1 

All  self  and  cross  hy drodynam i c  constants  (k)  =  0.01 

Species  transport  coefficients  are  given  in  table  2. 

Rotor  speed  =  60,000  rpm. 

Initial  sharp  boundary  position  (r  )  =  6.0  cm. 

At  =  2155  .51  sec 
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3.2.1    E f f e cl  of  Co n s tl t ue ai  Mole  Rati( 


For  completely  cooperative  AB   systems,  the  mole  ratio  at 

which  the  trailing  boundary  or  shoulder  is  smallest,  "MR   " 

eq  • 

seems  not  to  depend  upon  the  strength  of  association  (K 

nM ) 

Instead  *'Rg  appears  to  depend  only  on  the  s  t  o  i  ch  i  ome  try  of 
the  cooperative  system.  Cooperative  AB^  and  AB.  systems 
with  various  values  of  K^^  (and  identical  monomers)  had  min- 
imum trailing  boundaries'^  (or  shoulders)  at  mole  ratios  of 
2.0  and  4.0  respectively  (figures  37  -  39).  Thus,  in  gener- 
al, a  cooperative  AB^  system  (with  identical  monomers)  would 
be  expected  to  have  an  equivalence  point  at  a  mole  ratio  of 
equal  to  n. 

Since  uncooperative  AB^  systems  have  equivalence  mole 
ratios  of  less  than  n  (except  when  K.  is  very  large  -  see 
figure  21)  it  should  be  possible,  if  the  overall  stoichiom- 
etry  of  an  AB^  system  is  known  (ie.  if  n  is  known),  to  de- 
termine whether  the   system   is  uncooperative  or  completely 

cooperative  (if  it  is  one  or  the  other)   from  MR     alone. 

eq 

Also,  whether  or  not  n  is  known,  if  MR    is  not  a  whole  num- 

eq 

ber  greater  than  one,  then  one  can  probably  be  reasonably 
certain  that  the  system  in  question  is  not  completely  coop- 


13 


For  systems  with  bimodal  boundaries  at  their  equivalence 
points,  MR^  „4s  taken  as  the  mole  ration  at  which  H_  was 
smallest.  For  weaker  systems  that  were  not  bimodal  at 
^^eqj  ''"*  instead  had  a  trailing  shoulder,  MR  was  taken 
as  T;he  mole  ratio  at  which  the  trailing  sirkulder  was 
smallest,  approximated  by  superimposing  and  comparing 
profiles  at  different  mole  ratios. 


174 


Figure  37:   Gradient  Profiles  of  a  Completely  Cooperative 

AB^  Systems  at  Different  Constituent  Mole  Ratios 
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Figure  38:   Gradient  Profiles  of  a  Weak  Completely 

Cooperative  AB^  Systems  at  Different  Constituent 
Hole  Rat  i  os 
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Species  transport  coefficients  are  given  in  table  2.' 
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Figure  39:   Gradient  Profiles  of  Completely  Cooperative  AB 
Systems  at  Different  Constituent  Mole  Ratios 


^4M  =  1-Oxlo" 
W^  =  lOOKd 
Wg  =  lOOKd 

""AT  *  '^BT  "  ^^  •"  "g/"! 
[B]^  +  [A]^  =  1  .0x10"'' 

Frictional  ratios  ( f / f  )  for  all  species  =  1.1 
o  *^ 

All  self  and  cross  liy  drodynam  i  c  constants  (k)  =  0.01 

Species  transport  coefficients  are  given  in  table  2. 

Rotor  speed  =  60,000  rpm. 

Initial  sharp  boundary  position  (r  )  =  6.0  cm. 

At  =  2155  .51  sec 


[B]/[A] 

a.  3.00 

b.  4.00 

c.  5.00 


n.- 

a 

179 

so.' 
II.- 

■     -  A 
•     -B 

, ,  .-^ri«-'''                .It^^-  , 

s. 

M 

i.'"""' 

I.M                                                  1.30                                                 •-■/«                                                 I.M                                                 ^. 

-«*!l»5^ 


180 


owever,  and  MR^^  jj  an  integer 


erative.  If  n  is  not  known,  h 
greater  than  one,  then  the  system  in  question  could  be  ei- 
ther a  completely  cooperative  AB^  system  or  an  uncooperative 
AB   system  (where  m>n). 

The  changes  in  the  profiles  of  cooperative  systems  with 
mole  ratio  and  the  general  features  of  cooperative  profiles 
at  the  equivalence  point  are  distinctly  different,  qualita- 
tively, than  those  of  uncooperative  systems.  As  the  mole 
ratio  is  increased,  in  a  cooperative  system,  from  values  at 
which  there  is  an  excess  of  constituent  A  (and  most  of  the 
material  in  the  trailing  boundary  or  shoulder  is  constituent 
A)  the  constituent  B  gradient  begins  to  increase  its  contri- 
bution to  the  trailing  boundary  (or  shoulder)  before  all  of 
the  excess  constituent  A  is  eliminated  (ie.  before  the  con- 
tribution of  constituent  A  to  the  trailing  element  of  the 
boundary  stops  decreasing).  As  a  result,  the  trailing 
boundary  or  shoulder  of  a  cooperative  system  is  not  com- 
pletely eliminated  at  any  mole  ratio  as  it  is  with  uncooper- 
-g    .    Since   the  trailing  boundary  (or 


The  equivalence  point  occurs  at  the  mole  ratio  where  the 
increase  in  the  contribution  of  constituent  B  to  the 
trailing  boundary  (or  shoulder)  with  an  increase  in  MR  is 
equal  to  the  decrease  in  the  contribution  of  constituent 
A  to  the  trailing  element  with  the  same  increase  in  MR. 
For  systems  with  measurable  trailing  boundaries  this 
might  be  expressed  as: 


d(MR) 


dHj(A) 


d(MR) 
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shoulder)  does  not  disappear  in  cooperative  profiles  at  tie 
equivalence  point  as  it  does  in  uncooperative  profiles  it 
should  be  possible  to  determine  whether  a  system  is  com- 
pletely cooperative  or  uncooperative  (if  it  is  one  or  the 
other)  simply  by  examining  the  experimental  gradient  profile 
(schlieren  pattern)  at  the  equivalence  point.  So,  for  in- 
stance, when  MR  js  an  integer,  n,  greater  than  one  and  the 
e  q 

St  ox chiome try  is  not  known,  it  should  be  possible  to  deter- 
mine from  the  profile  at  the  equivalence  point  whether  the 
system  is  completely  cooperative  (AB  )  or  uncooperative  with 

a  higher  st o i ch i ome try  (AB   where  m>n). 

m 

3.2.2    Cooperative  Systems  with  Nonident  ical  Monomers 
Simulated  boundaries  of  C-4  systems  with  W,>Wu  (W,  =  140  Kd, 
Wg=60   Kd)   and   different   association   constants,   all   at 
MR=3.0,  are  shown  in  figure  40.   The  obvious  difference  be- 
tween these  profiles  and  any  others  considered  previously  is 


where  H  (a)  and  Hj(B)  are  the  heights  of  the  A  and  B 
stituent  gradients  at  the  peak  of  the  trailing  bound 
In  cooperative  systems  this  does  not  necessarily  occu 
the  point  where  the  trailing  edges  of  the  constit 
gradients  overlap.  In  fact  for  all  cases  studied 
trailing  edge  of  the  constituent  B  gradient  lies  a 
the  trailing  edge  of  cons t i t uent  A  gradient  at  the  eq 
alence  mole  ratio.  With  uncooperative  systems  essent 
ly  all  of  the  excess  A  is  eliminated  from  the  trai 
boundary  (increasing  MR)  and  the  trailing  edges  of 
constituent  gradients  meet  before  the  constituen 
boundary  begins  to  change,  and,  increasing  MR  beyond 
point  has  very  little  effect  on  the  trailing  edge  of 
constituent  A  boundary.  Thus,  for  uncooperative  syst 
MR  occurs  at  the  point  where  the  trailing  edges  of 
individual  constituent  profiles  overlap  the  most  (and 
total  boundary  is  least  skewed). 
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that  they  are  trimodal.   Constituent  gradients  show  that  the 

slowest   boundary   contains   constituent   B   only,   the   middle 

boundary   contains   mostly   constituent   A   and   the   leading 

boundary  contains  both  A  and  B  as  usual. 

As  the  constituent  mole  ratio  of  one  of  these  systems  is 

increased  (from  1.0)  the  fast  monomer  boundary  (A)  initially 

moving  at  about  the  rate  of  free  A  monomer  ( Ry.  (position  of 

trailing  boundary  containing  A)=6.35  cm.,   R  .  ,=«  34  cm.). 

mi  nA 

shrinks  and  migrates  at  progressively  faster  rates  until  it 
becomes  a  shoulder  on  the  trailing  edge  of  the  leading 
boundary  (figure  41).  Meanwhile,  the  slower  B  monomer 
boundary,  which  moves  at  a  rate  comparable  to  pure  B 
TB~^"^^  cm.;  ^maTB"^*^"^  cm.)  at  all  mole  ratios  consid- 
ered, steadily  becomes  larger. 

The  equivalence  point  for  these  systems  cannot  be  deter- 
mined accurately  from  either  the  total  gradient  or  constitu- 
ent gradient  profiles  because  differences  in  the  sizes  of 
the  total  trailing  element  at  mole  ratios  that  seem  to  be 
near  the  equivalence  point  (3.0<MR<5.0  in  figure  41)  are  am- 
biguous. Even  though  it  may  not  be  possible  to  determine 
the  equivalence  mole  ratio  of  this  type  of  system,  a  sam- 
pling of  gradient  profiles  at  various  mole  ratios  should  be 
helpful  in  distinguishing  between  completely  cooperative  and 
uncooperative  systems.  At  ME^^^  uncooperative  (AB^^)  sys- 
tems have  no  trailing  boundary  or  shoulder,  whereas  coopera- 
tive systems  have  at  least  a  trailing  shoulder  at  all  mole 
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Figure  40:   Gradient  Profiles  of  Completely  Cooperative  AB 
Systems  of  Different  Strengths  with  W,>w_ 


W^  =  140Kd 
Wg  =  60Kd 

''AT  ~  4.375  mg/ml 

[A]  =  3  .125x10"^  M 

^BT  ~  5.625  mg/ml 

[B]  =  9.375ilO"5  JI 
[B]/[A]  =  3 .0 

Frictional  ratios  (f/t^)     for  all  species  =  1.1 

All  self  and  cross  liy  drody  nam  i  c  constants  (k)  =  0.01 

Rotor  speed  =  60,000  rpm. 

Initial  sharp  boundary  position  (r  )  =  6.0  cm. 

At  =  1722.50  sec 

Species  transport  coefficients  are  given  in  table  3. 
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Figure  41:   Gradient  Profiles  of  Cooperative  AB   Systems 

with  Wj^>Wg  at  Different  Cons  t  i  t  nent  Mo  1  e  Ratios 


^4M  =  l.OxlO^' 
*A  =  140Kd 
Wg  =  SOKd 

''AT  ■^  '^BT  °  10  •"  ""g/"! 

Frictional  ratios  (ilf^)     for  all  species  =  1.1 

All  self  and  cross  hydrodynam i o  constants  (k)  =  0.01 

Rotor  speed  =  60,000  rpm. 

Initial  sharp  boundary  position  (r  )  =  g _o  cm. 

At  =  1722  .50  sec 

Species  transport  coefficients  are  given  in  table  3. 
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ratios  (except,  of  course,  when  K^j,  is  small  to  begin  with), 
and  at  least  some  cooperative  AB  systems  with  monomers  of 
significantly  different  sizes^^  have  trimodal  boundaries  at 
some  mole  ratios  which  uncooperative  systems  do  not  have  at 

any  mole  ratio. 

I 
'.J   .  ■ 
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Trimodal  boundaries  have  been  observed  for  AB  systems 
*'*''  *A"^'*°^  ^""^  Wg  =  60Kd  (figure  42),  as  well  as  for  sys- 
tems with  W^=60K  and  WB  =  140Kd   (figure  42),   but  not   for 

4   systems   with   monomers   that   differed   les 
(Wa=80K,  Wb  =  120K). 


m   size 
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Figure  42:   Gradient  Profile  of  a  Completely  Cooperative  AB, 
System  With  W„>W.  '' 


"a  =  60Kd 


Wg  =  140Kd 


''AT  =5.0  mg/ml 

[A]  =  8.333x10"^  M 
Cg.j.  =  5  .0  mg/ml 

[B]  =  3 .571x10"'  M 
[B]/[A]  =  0.429 

Frictional  ratios  (f/f  )  for  all  species  =  1.1 

All  self  and  cross  hydrodynamic  constants  (k)  =  0.01 

Rotor  speed  =  60,000  rpm. 

Initial  sharp  boundary  position  (r  )  =  6.0  cm. 

At  =  1722  .50  sec 

Species  transport  coefficients  are  given  in  table  4. 
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3 . 3    GENERAL  AB^  SYSTEMS 

It  has  previously  been  demonstrated  that  gradient  profiles 
of  uncooperative  and  cooperative  AB,  systems  are  generally 
not  alike.  Real  biological  AB^  systems,  of  course,  do  not 
all  necessarily  belong  to  one  of  these  two  classes  of  sys- 
tems. (In  fact  'completely  cooperative  systems  are  probably 
very  rare.)  In  the  full  range  of  possibilities,  uncoopera- 
tive systems  are  at  the  midpoint,  and  completely  cooperative 
systems  are  at  one  extreme  of  a  continuum  with  respect  to 
cooperativ ity .  At  the  other  end  of  the  continuum  are  anti- 
cooperative  systems,  the  extreme  of  which  would  be  an  AB 
system  in  which  binding  at  one  site  completely  prevents 
bonding  at  the  second  site.  This  case  is  indistinguishable 
from  a  simple  AB   system. 

Simulated  gradient  profiles  for  AB,  systems  with  various 
degrees  of  coope ra t i v i ty  and  ant i coope ra t i v i ty  were  studied 
to  see  if  boundary  shapes  differ  among  systems  with  differ- 
ent degrees  of  coope rat iv ity  as  they  do  among  completely  co- 
operative and  uncooperative  systems.  Plots  of  H  vs.  R 
(figure  43-a)  for  groups  of  systems  with  ^j^^/K' 2a  (an  index 
of  coopera ti V i ty^^ )  constant  and  K, 


For  uncooperative  AB^  systems  K,  „/K '  j  j,=  4  .0 ,  for  coopers 
tjve  systems  ^j  j|''^ '  2M^^  ■''  ^""^  *°'^  anti-cooperative  sys 
**'°*  ^im/K',„>4.0  .  wiiere 

^^    ^"  [AB] 


^"    [A][B] 


y    *>■ 
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appropriate  proportions  show  that,  at  least  when  the  associ- 
ations are  not  very  tight,  the  positions  and  heights  of  the 
leading  boundaries  are  together  unambiguously  diagnostic  for 
the  type  and  degree  of  coope  ra  t  i  v  i  ty  (Kj^jj/K-j^)-  For  exam- 
ple the  boundary  of  a  cooperative  system  with  K  „/K'   =2  5 

IM    2  H   * 

is  generally  different  from  that  of  a  more  cooperative  sys- 
tem (eg.  ^uj/K' 2j|=l  .0 )  or  of  a  less  cooperative  system  (eg. 
^1  jl/K '  2  jj-4  .0 )  as  long  as  the  association  is  not  very  strong. 
As  the  systems  become  tighter  the  curves  in  figure  43-a  be- 
gin to  run  together.  Plots  of  H^  vs.  R^  (fig.  43-b)  for  the 
same  groups  of  systems,  however,  show  that,  while  tight  sys- 
tems with  different  degrees  of  coope r a t i v i ty  may  have  lead- 
ing boundaries  with  the  same  heights  and  positions,  they 
will  probably  not  have  identical  trailing  boundary  heights. 
There  may  be  some  ambiguity  in  assigning  a  particular  model 
to  a  real  boundary  on  the  basis  of  its  trailing  boundary 
height  since  the  differences  between  H^s  of  systems  with 
different  values  of  K^jj/k'^j,  are  not  very  large,  but  it 
should,  at  least,  be  possible  to  distinguish  a  strongly  co- 


2M 
and 


[AB^] 
UB][B] 


2M 


=  KiM  ^  K' 


2M  - 


[ABj] 
[A][B]2 
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operative  system  from  an  uncooperative  or  anti-cooperative 
system. 

Since  the  values  of  K^jj  and  K'^^  are  unique  for  a  given 
degree  of  co  ope  ra  t  i  v  i  ty  (Kj^^/K'jj,)  and  association  strength 
'^1m''  '^  systems  with  different  degrees  of  coope ra t i v i ty 
and  association  strength  have  unique  boundary  shapes,   then, 

so   should   systems  with   different   values   of  K,  „   and  K'^u 

1 M  2  M  * 

The  independent  effects  of  K^^  and  K'jjj  on  boundary  shape 
are  illustrated  in  figure  44.  As  K  is  increased,  with 
2M  ^®^^  constant,  the  most  noticable  changes  in  boundary 
shape  are  an  increase  in  E^  and  a  decrease  in  H_.  Increas- 
^^6  ^'2M'  with  K^„  constant,  increases  the  leading  bounda- 
ry's sedimentation  rate.  The  leading  boundary  height  de- 
creases  as   K'^j,  is   increased   through   low  values   of  K'  „, 

then  passes  through  a  minimum  and  increases  as  K '  ,,  is  i  n- 

2  M 

creased  further.  The  effect  of  ^' ^n  on  H^  depends  on  the 
magnitude  of  ^^^.  At  low  values  of  K^j,,  the  trailing  bound- 
ary becomes  smaller  as  K'^j,  is  increased  and  at  large  values 
of  ^IM  '*  ^6<=o"6s  larger  as  K'jj,  is  increased. 

The   changes   in   H^^,   Hj,   and   R^   with   K'jji   for   several 

groups   of   systems  with   different  values  of  K,  „   (at   a  1:1 

1 M 

constituent  mole  ratio)  are  illustrated  in  figure  45. 
Plots  of  n^  vs.  Rj^  '^'2M  varied,  Kj^j,  constant)  for  the  same 
groups  of  systems  indicate  that  boundaries  of  systems  with 
different  values  of  Kj^jj  are  not  alike  at  low  and  intermedi- 
ate values  of  K'jj,,  but  are  similar  (with  respect  to  E^    and 
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Figure  43:   Comparison  of  Boundary  Dimensions  for  AB, 

Systems  with  Different  Degrees  of  Cooper  a t iv i ty 


W^  =  lOOKd 
"b  =  lOOKd 
''AT  ~  ^  •"  mg/ml 
''BT  ~  ^  •"  "g/i"! 

[A]  =  5x10"^  M 

[B]  =  5xl0"5  M 
[B]/[A]  =1.0 

Frictional  ratios  d/f^)     for  all  species  =  1.1 

All  self  and  cross  hydrodynamic  constants  (k)  =  0.01 

Species  transport  coefficients  are  given  in  table  2. 

Rotor  speed  =  60,000  rpm. 

Initial  sharp  boundary  position  (r  )  =  g.o  cm. 

At  =  2155  .51  sec 

^-     Hj^    vs.     Rj^       for    families    of    AB2     systems 
*'*''    ^1M''^'2M    constant 


!'•    H^    vs.     Rj_       for    families    of    ABj    systems 
"i'Ji    ^Ijl/K'jji    constant 


K|m/K2m 


194 


po.o- 


R^^cn\) 


X 


K|m/K; 


Rjicm) 


195 


Figure  44:   Gradient  Profiles  of  AB   Systems  w  i  tli  K    and 
^'2M  Varied  Independently 


W^  =  lOOKd 


"b  =  lOOKd 


^j^j    =    5.0  mg/ml 

BT  ~  ^  '^    mg/ml 

[A]  =  5ilO~5  M 

[B]  =  5xlO"5  M 
[B]/[A]  =1.0 

Frictional  ratios  (f/f^)  for  all  species  =  1.1 

All  self  and  cross  hydrodynamic  constants  (k)  =  0.01 

Species  transport  coefficients  are  given  in  table  2. 

Rotor  speed  =  60,000  rpm. 

Initial  sharp  boundary  position  (r  )  =  g.o  cm. 

At  =  2155  .51  sec 


»•  K^jj  varied,  K'^jj  =  5.0x10'' 
''■  ^IM  varied,  K'jj,  =  1.0x10^ 
" '     '^'2M  varied,  Kj^j,  =  5.0x10'* 
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Kj^)  at  higher  values  of  K'jjj,  where  these  curves  converge 
(see  figure  46-a).  Plots  of  H^  vs.  R^  (figure  46-b)  demon- 
strate additional  differences  in  boundary  shape  that  seem  to 
extend  to  somewhat  tighter  systems.  These  differences  among 
the  tighter  systems  are,  however,  fairly  small  and  may  not 
be  practically  useful. 

Species  concentration  vs.  logj^K'^M  curves  (figure  47) 
show  that  at  large  values  of  K' ^^  there  is  essentially  no 
difference  in  the  compositions  of  systems  with  different 
values  of  ^^.w  Average  constituent  sedimentation  coeffi- 
cients are,  therefore,  also  similar  at  large  K'  s  (figure 
48).  In  view  of  these  similarities  it  is  not  too  surprising 
that  the  sedimenting  boundary  profiles  of  these  systems  are 
also  alike. 

In  some  cases,  where  systems  with  significantly  different 

values   of  ^^ff     and  K'2j||  have  very   similar  boundaries   at   a 

constituent  mole  ratio  of  1:1,  they  may  have  more  noticeably 

different  boundaries   at  another   ratio.    For   example,   at   a 

1:1  mole  ratio  the  boundary  of  the  system  with  K,„=2xlo'  and 

IM 

^'2)1  =  2.5x10^  is  nearly  identical  to  that  of  the  system  with 
Kj^y=5xl0'  and  K'2j,=  3xlO^  (see  figure  49-a),  but,  at  a  con- 
stituent ratio,  [B]:[A],  of  2:1  the  same  two  systems  produce 
obviously  different  boundaries  (fig.  49-b). 

At  a  2:1  constituent  ratio  the  effects,  on  boundary 
shape,  of  changes  in  K'^j,  (K^j,  constant)  are  qualitatively 
different   than  those   at  a  1:1  mole  ratio   (see  figure  50). 
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Figure  45:   Changes  ia  Boundary  Dimensions  with  K'     for  AB 


Systems  with  K^^    Constant 


2H 


W^  =  lOOKd 
"b  =  lOOKd 

*'BT  ~  ^  •"  ■"g/'ol 

[A]  =  5x10"^  M 

[B]  =  5xlO"5  M 
[B]/[A]  =1.0 

Frictional  ratios  (f/(  )  for  all  species  =  1.1 

All  self  and  cross  hydrodynamic  constants  (k)  =  0.01 

Species  transport  coefficients  are  given  in  table  2. 

Rotor  speed  =  60,000  rpm. 

Initial  sharp  boundary  position  (r  )  =  6.0    cm. 

At  =  2155.51  sec 

a-     Hj^    vs,     l''8io^'2M    ^°'    families    of 
systems    with    K^^^    constant 


b.     H^    vs.     loSio^'2M    ^°'^    families    of 
systems    with    K^j^    constant 


"■     R^    ^'-     l°8io^'2M    ^°^     families    of 
systems    with    K^^j    constant 
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Fignre  46:   Comparison  of  Boundary  Dimensions  for  AB, 


Systems  with  Different  values  of  K' 


2U 


Wj^  =  lOOKd  • 

Wg  =  lOOKd  -   . 

C^j  =5.0  mg/ml 

Cg J  =  5  .0  mg/ml 

[A]  =  5x10"^  M 

[B]  =  5x10"'  M 
[B]/[A]  =  1 .0 

Frictional  ratios  ( f / f  )  for  all  species  =  1.1 

All  self  and  cross  liy  drody  nam  i  c  constants  (k)  =  0.01 

Species  transport  coefficients  are  given  in  table  2. 

Rotor  speed  ==  60,000  rpm. 

Initial  sharp  boundary  position  (r  )  =  g.o  cm. 

At  =  2155  .51  sec 


systems  w  i  th  K' 


2H 


constant 


*>•  ^j    vs.  Rj^   for  families  of  ABj 
systems  with  K'^j,  constant 
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Figure  47:   Species  Concentrations  for  AB   Systems 


'^  =  lOOKd 


Wg  =  lOOKd 


AT 


5 .0  mg/ml 


[A]  =  5il0"5  M 

[B]  =  5xl0~5  H 
[B]/[A]  =  1.0 


a.  Species  concentrations  vs.  log   f" 

lU   2M 

systems  with  K, „  =  5.0x10* 

IN 


for 


b.  Species  concentrations  vs.  log   K'«u  for 
systems  with  K.j,  =  2.0x10^ 


c.  Species  concentrations  vs.  ^"^-in^'  ju    for 
systems  with  K    =  5.0x10^ 


Co(mg/ml) 
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Figure  48:   Average  Constituent  Sedimentation  Coefficients 
for  AB^  Systems 


W^  =  lOOKd 
*B  =  lOOKd 
''AT  ~  ^  •"  ""g/""! 
*'BT  ~  ^  •"  ""g/nl 
[A]  =  5il0~5 


[B]  =  5x10 


-5 


[B]/[A]  =  1 .0 

Frictional  ratios  (f/f  )  for  all  species  =  1.1 

Species  transport  coefficients  are  given  in  table  2, 

a.  Sbar^  vs.  K'-ji  for  systems  with 
different  values  of  K. „ 


b.  SbaTg  vs.  E'jM  for  systems  with 

different  values  of  K,,, 
IM 
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Figure  49:   Gradient  Profiles  of  Two  AB,  Systems 
Different  Constituent  Ratios 


it  Two 


Wj^  =  lOOKd 

Wg  -  lOOKd 

Frictional  ratios  (f/f^)  for  all  species  =  1.1 

All  self  and  cross  liy  drody  nam  i  c  constants  (k)  =  0.01 

Species  transport  coefficients  are  given  in  table  2. 

Rotor  speed  =  60,000  rpm. 

Initial  sharp  boundary  position  (r  )  =  6,0  cm. 

At  =  2155  .51  sec 

a.  C    =5.0  mg/ml 


BT 
[A] 


=  5.0  mg/ml 


5x10 


[B]  =  5x10  ^ 
[B]/[A]  =  1.0 


AT 

r 
BT 

[A]  =  3  .33x10" 


[B] 


6  .67x10 


-5 


[B]/[A]  -  2.0 
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At  a  2:1  mole  ratio  H^^  increases  steadily  as  K'2j|  is  in- 
creased, rather  than  passing  through  a  ninimum  as  it  does  at 
a  1:1  ratio,  and  H^  decreases  with  increasing  K'2in  at  a  2:1 
constituent  ratio  whereas,  at  a  constituent  ratio  of  1:1,  it 
decreases  if  K  „  js  small  (<  ca.  3x10^)  but  increases  if 
■^Ijl  is  larger. 

Plots  of  E^  vs.  R^  and  H^  vs.  R^  for  AB2  systems  at  a  2:1 
constituent  ratio  (K'^j,  varied,  Kj^j,  constant;  figure  51) 
suggest  that,  in  general,  boundaries  of  Afi^  systems  with 
different  association  constants  may  be  more  readily  distin- 
guishable, when  K'  is  large,  at  a  2:1  ratio  than  at  a  1:1 
ratio. 


.-,  l> 
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Figure  50:   Changes  in  Boundary  Dimensions  With  K'     for  AB 


Systems  with  K^^    Constant 


2H 


Wj^  =  lOOKd 

Wg  =  lOOKd 

^AT  ^  3.33  mg/ml 

^BT  ~  6-67  mg/nl 

[A]  =  3  .33ilO"5  M 

[B]  =  6.67x10"'  M 
[B]/[A]  =2.0 

Frictional  ratios  (f/f^)     for  all  species  =  1.1 

All  self  and  cross  hydrodynamic  constants  (k)  =  0.01 

Species  transport  coefficients  are  given  in  table  2. 

Rotor  speed  =  60,000  rpm. 

Initial  sharp  boundary  position  (r  )  =  6.0    cm. 

At  =  2155  .51  sec 

^-  H^  ■'*•  '°8ioK'2M  ^°'     families  of 

systems  with  K^^    constant  at  [B]/[A]  =  2.0 


b .  H^  vs.  logio^'2M  ^"^  families  of 

systems    with    K^j,    constant    at     [B]/[A]     =    2.0 

<=  •     \    vs.     lo8ioK'2M    for    families    of 

systems    with    K^^    constant    at     [B]/[A]    =    2.0 
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Figure  51:   Comparison  of  Boundary  Dimensions  for  AB 


Systems  with  Different  Values  of  K  „ 
Constituent  Mole  Ratio  ^ 


It  a  2:1 


W^  =  lOOKd 


B 


=  lOOKd 


^AT  "^  3.33  mg/ml 

*'BT  =  6  -67  mg/ml  "'  ; 

[A]  =  3 .33xlO~5  M 

[B]  =  6.67xl0~5  M 
[B]/[A]  =  2  .0 

Frictional  ratios  (f/f^)  for  all  species  =  1.1 

All  self  and  cross  hy drodynam i c  constants  (k)  =  0.01 

Species  transport  coefficients  are  given  in  table  2. 

Rotor  speed  =  60,000  rpm. 

Initial  sharp  boundary  position  (r  )  =  fi.o  cm. 

At  =  2155  .51  sec 


a-  H^  vs.  Rj^       for  families  of  AB^ 

systems  with  K' ^^    constant  at  [B]/[A]  =  2.0 

t"  •  Hj  vs.  Rj^   for  families  of  AB^ 
'2M 
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3 . 4    EFFE CT  OF  AGGREGATE  FRICTIONAL  RAH OS 

Since  the   sedimentation  and  diffusion  coefficients  of  each 
species  ij  in  a  mixed  associating  system  (with  monomers  of  a 

given  volume  (Mv/N))  depend  on  the  molecules  frictional  ra- 

1  7 
tio     '^'^o'ii'   *°  does  the  shape  of  the  boundary  produced 

by  the  system  in  a  velocity  sedimentation  experiment. 

The   frictional   ratio  of   each  monomer,   A  and  B,   can  be 

computed   from   its   diffusion   coefficient   and   the   minimum 

frictional  coefficient. 

RT      f 

DNf     f 
o     0 

or   from   the  Stokes   radius   (R^  =  t/6m\)     obtained  from  gel 
permeation  and  R   (K„  =  (  3Mv/ 4tiN)  ^^  ^  )  . 


The   sedimentation  and  diffusion   coefficients   of   a  mol- 
ecule are  related  to  its  frictional  ratio  as  follows: 
M(l  -  vp)                 RT 
S  = D  = 

'^^<f''fo>  Nf„(f/f„) 


where  f / f ^  is  the  ratio  of  the  frictional  coefficient  of 
the  molecule  itself  to  that  of  an  unhydrated  spherical 
particle  of  equal  volume  (with  radius  R  \ 


where 
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f/6nil 
o/6jiTi 


The  frictional  ratios  of  the  aggregates,  however,  cannot 
be  obtained  by  these  methods  since  the  individual  aggregates 
cannot  be  isolated  from  other  species  in  a  mixed  associating 
system.  It  may,  therefore,  be  necessary,  when  modeling  ve- 
locity sedimentation  of  mixed  associating  systems,  to  con- 
sider systems  with  a  variety  of  plausible  combinations  of 
aggregate  frictional  ratios. 

Simulated  boundaries  for  uncooperative  AB,  systems 
(Kj=10  )   in   which   the   frictional   ratios   of   AB   and   AB, 

o^  fiB  """^  '^'^o'aB  ^^^  varied,  but  kept  equal  to  one 
another,  are  shown  in  figure  52,  As  one  would  expect,  in- 
creasing the  aggregate  frictional  ratios  sharply  decreases 
the  average  sedimentation  rate  of  the  reaction  boundary,  but 
does  not  affect  the  sedimentation  rate  of  the  trailing  mo- 
nomer boundary.  The  heights  of  both  the  leading  and  trail- 
ing boundaries  are  also  increased  slightly.   If  ( f / f  )    and 

^''  ^o^AB  ^^®  made  large  enough,  the  boundary  becomes  unimo- 
dal,  with  a  trailing  shoulder  that  becomes  less  prominent  as 
the  frictional  ratios  are  increased  further.  Thus  systems 
in  which  the  aggregates  are  highly  asymmetrical  and/or  hyd- 
rated  (ie,  large  f / t ^)  are  likely  not  to  generate  bimodal 
boundaries  unless  the  association  is  relatively  tight. 
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Figure  52:   Gradient  Profiles  of  Hnooope ra t i v e  AB,  Systems 
with  Different  Aggregate  Frictional  Ratios 


100  Kd 


"b  =  100  Kd 


'"AT  =5.0  mg/ml 
''BT  ~  ^  '^  mg/ml 
[A]  =  5x10"'  M 


[B]  =  5ilO" 


H 


[B]/[A]  =1.0 

All  self  and  cross  liy  drody  nam  i  c  constants  (k)  =  0.01  ml/i 

Species  transport  coefficients  are  given  in  table  2. 

Rotor  speed  =  60,000  rpm. 

Initial  sharp  boundary  position  (r  )  =  6,0  cm. 

At  =  2155  .51  sec 
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While  it  may  be  that,  among  AB.  systems  with  the  same  ag- 
gregate frictional  ratios,   systems  with  different  values  of 

^m  and  K'2j(  generally  generate  boundaries  of  different  di- 
mensions or  shapes  (unless  the  association  is  very  tight  - 
see  p.  192),  it  does  not  seem  unlikely  that  systems  with 
different  aggregate  frictional  ratios,  as  well  as  different 
association  constants,  would,  in  some  cases,  generate  iden- 
tical or  very  similar  boundaries.  Thus,  when  fitting  model 
system  boundaries  to  the  observed  boundaries  of  real  sys- 
tems, satisfactory  fits  might  be  obtained  with  two  or  more 
different  systems  (ie.  systems  with  different  combinations 
of  frictional  ratios  and  association  constants).  Different 
systems  that  have  similar  boundaries  at  one  constituent  mole 
ratio,  however,  might  not  have'  such  similar  boundary  shapes 
at  another  mole  ratio.  Some  incorrect  models  might,  there- 
fore, be  eliminated  by  modeling  at  different  mole  ratios. 

3  . 5    EFFECT  OF  HYDRODYNAMIC  DEPENDENCE 

Hydrodynamic  concentration  dependence  for  sedimentation 
coefficients  is  described  in  the  simulation  model  by  the 
following  express  ions. 


^B  =Sb,o/(1  *    ''BBCbT  +  ^baCat' 


wher  e  "5  . 


sedimentation  coefficients,   C    and  C 
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■"A,  0  ^""^  ^B  0  ^^^     *^*  ideal  local  average  constituent 

ire  tlie  local  con- 
stituent concentrations,  k^^  a„d  t^^  are  the  self  hydrody- 
namic  constants  and  k^g  and  kg^  are  the  cross  hydrodynamic 
constant  s . 

In  all  of  the  simulated  experiments  discussed  so  far  all 
four  of  the  hydrodynamic  constants  had  values  of  .01  (.01  is 
about  average  for  real  systems  in  which  actual  values  may 
vary  from  a  low  of  about  .005  to  an  upper  limit  of  about 
.020).  The  effects  of  hydrodynamic  dependence  on  boundary 
shapes  of  AB^  systems  were  studied  by  varying  hydrodynamic 
constants  in  three  different  ways. 

I)  The  effect  of  general  hydrodynamic  dependence  was 
investigated  by  varying  all  of  the  constants  to- 
gether (ie.  keeping  all  the  hydrodynamic  constants 
equal ) . 

II)  The  effect  of  general  cross  hydrodynamic  dependence 
was   studied  by  varying  k     and  k  .   together   (ie. 


AB-'^BA' 


)  with  k,,  and  kop  constant  at  .01, 


'AA 


"■BB 


III)   The  specific  effects  of  each  of  the  four  dependen- 
cies  were  studied  by  varying  k^^_  kgg^  k^j  and  kg^ 
one   at  a   time  while  keeping  the  other   three  con- 
stant a t  . 01  . 
Simulated  profiles  of  uncooperative  AB,  systems  (K,=lo'. 
^ f^=^^  =  i.0  0t.i,      Cj^x~''BT~^  •"  mg/ml)  with  various  degrees  of  hy- 
drodynamic dependence   (all   k's   equal)   are   shown  in  figure 


219 

53-a.  As  one  might  expect,  an  increase  in  overall  hydrody- 
nam i c  dependence  results  in  a  sharpening  of  both  the  leading 
and  trailing  boundaries  as  well  as  a  decrease  in  their  aver- 
age sedimentation  rates.  Because  the  leading  boundaries' 
sedimentation  rate  is  more  strongly  affected  than  that  of 
the  trailing  boundary,  increasing  overall  hydrodynamic  de- 
pendence brings  the  leading  and  trailing  boundaries  closer 
together.  For  weaker  systems  this  means  that  the  distin- 
guishing bimodal  characteristics  of  the  sedimenting  boundary 
may  be  obscured  by  large  hydrodynamic  effects  (Figure 
53-a,K  =10  ).  However,  if  velocity  sedimentation  is  allowed 
to  continue  for  a  longer  time,  bimodalty,  or  related  charac- 
teristics, that  are  obscured  by  hydrodynamic  dependence, 
might  eventually  appear. 

In  general,  if  a  boundary  is  leaning  toward  being  bimodal 
(ie.  if  it  has  four  curvature  inflections),  then  it  is  capa- 
ble of  becoming  bimodal  at  a  later  time,  once  it  has  sedi- 
mented  through  a  distance  large  enough  to  separate  the  two 
components.  However,  in  a  real  situation,  if  bimodality  is 
to  be  observed  then  it  must  develop  before  the  boundary  be- 
comes worthlessly  distorted  by  end-of-cell  effects. 

Hydrodynamic  dependence  seems  to  obscure  bimodality,  as 
already  noted,  because  it  causes  the  two  components  of  the 
boundary  to  separate  more  slowly.  However,  if  there  is  a 
large  hydrodynamic  effect,  the  leading  and  trailing  bounda- 
ries,  or   their   rudiments,   will   be   relatively   sharp,   and 


'"^ 
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Figure  53:   Effect  of  Overall  Hydrody nam i c  Dependence  on 
Boundaries  of  Uncooperative  AB-  Systems 


W^  =  lOOKd 
"b  =  lOOKd 
C^^  =5.0  mg/ml 

''BT  ~  '  •"  "S/""! 

[A]  =  5x10"^  M 

[B]  =  5x10"'  M 
[B]/[A]  =  1.0 

Frictional  ratios  ( f / f  )  for  all  species  =  1.1 

Species  transport  coefficients  are  given  in  table  2, 

Rotor  speed  =  60,000  rpm. 

Initial  sharp  boundary  position  (r  )  =  6.0  cm. 

At  =  2155  .51  sec 


a.  Kj  =  1  .0x10^  M  ^ 


-----I 


b.  Kj  =  5  .0x10'*  M~l 
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should,  therefore,  become  resolved  at  relatively  small  sepa- 
rations. Comparisons  of  boundaries  of  otherwise  identical 
systems  with  different  k's,  which  were  allowed  to  sediment 
until  their  leading  boundaries  had  migrated  through  equal 
distances  (ie.  longer  times  for  higher  values  of  k )  demon- 
strate that  reasonable  and  even  unreasonable  hydrodynamic 
effects  do  not  obscure  bimodality.  In  fact,  in  some  in- 
stances (figure  54),  systems  with  higher  dependencies  are 
more  likely  to  have  developed  a  bimodal  boundary  by  the  time 
their  leading  boundaries  have  reached  a  given  position  than 
are  systems  with  lower  dependency.  Thus,  in  at  least  some 
cases,  hydrodynamic  effects  (boundary  sharpening)  can  actu- 
ally help  to  bring  out  bimodal  tendencies  in  a  boundary. 

Cross  hydrodynamic  dependence  (figure  55)  has  the  same 
effect,  qualitatively,  on  the  leading  boundary,  as  general 
dependence  (ie.  the  leading  boundary  becomes  sharper  and 
moves  at  a  slower  rate).  The  effect  on  the  leading  boundary 
of  changing  k^^  and  kg^^  from,  for  instance,  .010  to  .015, 
while  being  the  same  qualitatively,  is,  of  course,  consider- 
ably less  drastic  than  the  effect  of  increasing  all  of  the 
hydrodynamic  constants  by  the  same  amount.  Increasing  cross 
hydrodynamic  dependence  also  increases  the  trailing  bounda- 
ry's height  (as  does  general  hydrodynamic  dependence)  but 
has  very  little  if  any  effect  on  either  its  migration  rate 
or  sharpness.  That  there  is  no  effect  on  the  trailing 
boundary's  migration  rate  or  width  is,   of  course,  not  sur- 
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Figure  54:   Effect  of  Overall  Hydrodynam i c  Dependence  on 
Bimodality  of  Uncooperative  AB   Profiles 


boundaries  of  systems  with  different 
hydrodynamic  dependences  after  migrating 
through  equal  distances. 


Kj  =  3  .5x10* 
W^  =  lOOKd 
Wg  =  lOOKd 

''AT  "  '  •"  ""g/""! 

''BT  =  5  .0  mg/ml 

[A]  -  5x10"'  M 

[B]  =  5xlO"'  M 

[B]/[A]  =1.0  ,    . 

Frictional  ratios  ( f / f  )  for  all  species  =  1.1 

Species  transport  coefficients  are  given  in  table  2. 

Rotor  speed  =  60,000  rpm. 

Initial  sharp  boundary  position  (r  )  =  6.0  cm. 

At  =  2155  .51  sec 
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prising,  since  the  trailing  boundary  consists  of  only  one 
constituent  (A)  and  cross  hy drodynam i c  effects  should  there- 
fore be  absent  there. 

Since  leading  boundaries  of  systems  with  large  cross  hy- 
drodynamic  effects  move  more  slowly,  in  relation  to  their 
trailing  elements,  than  do  those  of  systems  with  smaller 
cross  effects,  bimodal  behavior  might  be  obscured  in  rela- 
tively weak  systems  by  cross  hydrodynamic  dependence  alone, 
particularly  since  it  does  not  sharpen  the  trailing  bounda- 
ry. .   ■(-_ 

The  independent  effects  of  each  of  the  two  self  and  two 
cross  hydrodynamic  constants  (k^^  and  kgg,  and  k^g  and  kg^) 
are  illustrated  in  figures  56  and  57. 

The  effects  of  the  four  types  of  hydrodynamic  dependence  on 
boundary  shape  are  not  each  distinct.  In  general  the  ef- 
fects of  cross  or   self  hydrodynamic  dependence   on  ?.   (je 

A 

the  effects  related  to  k..  and  k.p)  are  very  similar  as  are 


the  effects  of  cross  and  self  hydrodynamic  dependence  on  "S- 

6 

(ie.  the  effects  related  to  k^^  and  kg^)  Both  self  and  cross 
dependence  of  ^^  are  associated  with  higher  trailing  bounda- 
ries and  lower,  narrower  leading  boundaries,  while  cross  and 
self  dependence  of  5g  are  associated  with  lower,  broader 
trailing  boundaries  and  taller,  sharper  leading  boundaries 
as  well  as  slower  leading  boundary  migration  rates.  Self 
dependence  of  S^  C^aa'  also  slows  down  and  sharpens  the 
trailing  boundary.   None  of  the  other  three  dependences  af- 


Figure  55:   Effects  of  Cross  Hydrodynam i c  dependence  on 
Boundaries  of  Uncooperative  AB»  Systems 
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^f^    =  lOOKd 
Wg  =  lOOKd 

''AT  ^  ^  ■"  ""g/""! 
Cg-j.  =  5  .0  mg  /ml 

[A]  =  5x10"^  M 

[B]  =  5x10"^  M 
[B]/[A]  =1.0 

Frictional  ratios  (f/f  )  for  all  species  =  1.1 

Species  transport  coefficients  are  given  in  table  2. 

Rotor  speed  =  60,000  rpm. 

Initial  sharp  boundary  position  (r  )  =  6.0    cm. 

At  =  2155  .51  sec 
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Figure  56:   Effects  of  Individual  Self  Hydrodynamic 

Dependences  on  Boundaries  of  Uncooperative  AB 
Sy  St  ems 


W^  =  lOOKd 
Wg  =  lOOKd 

""AT  ^  ^  •"  "g/"! 

C 
BT 

[A]  =  5x10"^  M 

[B]  =  5x10"^  M 
[B]/[A]  =1.0 

Frictional  ratios  ( f / f  )  for  all  species  =  1.1 

Species  transport  coefficients  are  given  in  table  2, 

Rotor  speed  =  60,000  rpm. 

Initial  sharp  boundary  position  (r  )  =  fi.o  cm. 

At  =  2155  .51  sec 


^*  ^AA  v^^i®<i»  *11  other  ks  =  .01 


b.  kgg  varied,  all  other  ks  =  ,01 
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feet  either  the  trailing  boundary's  migration  rate  or  its 
width,  nor  should  they,  since  the  trailing  boundary  contains 
constituent  A  only. 


Since  hydr ody nam i c  dependence  affects  the  position  and 
height  of  the  leading  boundary  it  is  possible  that  some  sys- 
tems with  different  combinations  of  association  constants 
that  also  have  different  combinations  of  self  and  cross  hy- 
drodynamic  constants  would  give  rise  to  identical  bounda- 
ries. However,  whereas  changes  in  association  constants 
that  slow  down  the  leading  boundary  generally  either  widen 
the  leading  boundary  or  do  not  affect  its  width  at  all,  hy- 
drodynamic  dependence,  in  addition  to  slowing  down  the  lead- 
ing boundary,  also  decreases  its  width.  It  might  therefore 
be   possible   in  modeling   real   mixed   systems   to   separate, 

somewhat,  the  effects  of  hydrodynamio  dependence  and  associ- 

1  8 
ation   .    The   situation  is   likely   to  be   complicated,   how- 
ever, by  a  lack  of  knowledge  about  the  frictional  ratios  of 
the  aggregates,  since  these  may  also  affect  the  width  of  the 
leading  boundary. 


18 


The  self  hydrodynamic  constants  k  ^nd  kno  would  usually 
be  obtained  more  easily  from  modeling  studies  with  the 
individual  constituents. 
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Figure  57:   Effects  of  Individual  Cross  Hydrodynam i c 

Dependences  on  Boundaries  of  Uncooperative  AB, 

Systems        ■  ,  ^ 


Kj  =  1  .0x10^ 
"^  =  lOOKd 
Wg  =  lOOKd 

''AT  °  ^  ■"  ""g/"! 
Cg'j'=5,0mg/ml 

[A]  =  5x10"^  M 

[B]  =  5ilO"5  M 
[B]/[A]  =1.0 

Frictional  ratios  ( f / f  )  for  all  species  =  1.1 

Species  transport  coefficients  are  given  in  table  2. 

Rotor  speed  =  60,000  rpm. 

Initial  sliarp  boundary  position  (r  )  =  6.0  cm. 

At  =  2155  .51  sec 


^*  ^AB  ^^ried,  all  other  ks  =  .01 


^-  ^BA  varied,  all  other  ks  =  .01 
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3 . 6    EFFE CT  OF  SELF  AS SO CI ATION 

Occasionally  one  or  both  of  the  constituents  of  a  mixed  as- 
sociating system  will  self  associate.  The  effects  of  self 
association  of  A  and  B  to  dimer^ 

2A  <  =  >  A„ 


[Aj] 

^1^  ~  Tkv 


and 


2B  <=>  B 


K     =  [III 
^2"    [B]^ 


on  the  sedimentation  of  profiles  of  AB,  systems  were  studied 
using  the  COXMIX  distorted  grid  simulation  program  for  mixed 
associating  systems  and  a  modification  of  the  SBABN  table 
assembly  program  accommodating  self  association  of  either  or 
both  monomers  to  the  dimer. 

Simulated  profiles  of  uncooperative  AB   systems  (K,=lo'), 
at  a  constituent  ratio  of  1:1,  with  various  values  of  K 
and  K. 

Increasing  K     j^  these  systems,  at  this  mole  ratio,  has 
a  strong  effect  on  the  shape  of  the  trailing  boundary.   Even 


^2 


^3       are     shown    in    figures    58    and    59. 
2 
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at  relatively  small  values  of  K.  (5x10^).  the  trailing 
boundary  is  considerably  lower  and  broader  than  it  is  in  the 
absence   of   self  association,   probably  due   to  spreading  of 

the  constituent  A  boundary.   Also,   as  K    is  increased,  the 

*2 
average  sedimentation  rate  of  the  trailing  boundary  increas- 
es, probably  because  of  the  contribution  of  a  dimer  gradient 
to  the  trailing  boundary.  The  leading  boundary's  height  and 
position  are  not  affected.  Apparently  only  the  excess  A  in 
the  trailing  boundary  is  significantly  involved  in  self  as- 
sociation. / 

**  ^A    is  increased  and  the  difference  between  the  aver- 
2 

age  migration  rates  of  the  leading  and  trailing  boundaries 
becomes  smaller,  the  bimodal  appearance  of  the  boundary  is 
lost.  Hence,  self  a.ssociation  may  obscure  the  distinctive 
features  of  an  AB^  boundary,  particularly  if  the  mixed  asso- 
ciation is  weak. 

Introducing  self  association  of  B  to  the  system,  at  the 
same  constituent  concentrations,  has  relatively  little  ef- 
fect  on   boundary   shape,   except   at   large   values   of   K 

4  ^2 

(5x10  ),  where  the  size  of  the  trailing  boundary  is  in- 
creased and  the  average  sedimentation  rate  and  size  of  the 
trailing  are  decreased  (figure  59)  Apparently,  because 
there  is  no  B  in  the  trailing  boundary,  self  association  of 
B  cannot  affect  the  boundary  shape  unless  it  is  strong 
enough  to  significantly  alter  the  availability  of  B  monomer 
for  the  mixed  association. 
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Figure  58:   Effect  of  Self  Association  of  A  on  the  Boundary 
)f  an  Uncooperative  AB   System  at  a  1:1 


Constituent  Ratio 


Kj.  =  1  .0x10^  M  ^ 
W^  =  lOOKd 


lOOKd 


''AT  '^    ^  •'^    mg/ml 
Cg^  =  5.0  mg/ml 


[A] 


5x10"^  M 


[B]  =  5x10"^  M 

[B]/[A]  =  1 .0 

Frictional  ratios  ( f / f  )  for  all  species  =  1.1 
0  ^ 

All  self  and  cross  hydrodynamic  constants  (k)  =  0.01 

Species  transport  coefficients  are  given  in  table  2. 

Rotor  speed  =  60,000  rpm. 

Initial  sharp  boundary  position  (r  )  =  6.0  cm. 

At  =  2155  .51  se  o 

a.     no     self    association 

>>■     K         =    5.0x10^    M-1 

*2 
=•     K.        =    2.5x10"    M-1 


,,  1 
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Figure  59:   Effect  of  Self  Association  of  B  on  the  Boundary 
of  an  Uncooperative  AB   System  at  a  1:1 
Const  i  tuent  Ratio 


Kj  =  1.0x10^  M~^ 
Wj^  =  100  Kd 
Wg  =  100  Kd 
C^j  =  5.0  mg/ml 

^BT  ~  ^"^  ng/nl 
[A]  =  5xlO~5  M 


[B]/[A]  =  1.0 

Frictional  ratios  ( f / f  )  for  all  species  =  1.1 

All  self  and  cross  hy drody nam i c  constants  (k)  =  0.01 

Species  transport  coefficients  are  given  in  table  2. 

Rotor  speed  =  60,000  rpm. 

Initial  sharp  boundary  position  (r  )  =  6,0  cm. 

At  =  2155  .51  sec 

a.  no  self  association 

b.  K    =  5.0x10''  M-1 

"2 
<:•  K    =  s.Oxlo'  M-1 
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If  the  const! tuent  concentrations  are  changed,  making  B 

the  constituent  in  excess,  the  effect  of  self  association  of 

B  is  much  more  obvious  at  low  values  of  K    (SxlO^)-   With  B 

in  excess  the  effects,  on  the  boundary,  of  an  increase   in 

Kr   are  about  the  same  as  the  effects  of  increasing  K*   when 

2  ^2 

A  is  in  excess   (ie.   the  trailing  boundary  becomes  shorter 

and  broader  and  moves  at  a  faster  rate  so  that  it  merges 
with  the  leading  boundary  -see  figure  61)  At  these  constitu- 
ent concentrations  increasing  K.    has  no  noticeable  effect 

^2 
on  the  boundary  shape  with  values  of  K.   as  large  as  5x10   '» 

*2 
ie.  The  effect  is  the  same  as  increasing  Kq   when  A  is  in 

excess  (figure  60). 

Since,  when  one  or  the  other  of  the  constituents  is  in 
excess,  the  average  migration  rate  of  the  trailing  boundary 
is  affected  only  by  self  association  of  the  constituent  in 
excess,  and,  since  the  average  migration  rate  of  trailing 
boundary  is  usually  constant  (among  otherwise  identical  sys- 
tems with  different  mixed-association  constants  and  aggre- 
gate frictional  ratios),  when  one  of  the  constituents  is  in 
excess  the  trailing  boundary  postion  or  its  average  migra- 
tion rate  might  be  useful  in  modeling  studies  to  assess  self 
association.  Usually,  however,  the  self  association  con- 
stants for  A  and  B  can  be  obtained  from  studies  with  the  in- 
dividual constituents. 
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Figure  60:   Effect  of  Self  Association  of  A  on  the  Boundary 
of  an  Uncooperative  AB   System  at  a  2:1 
Const  itnent  Ratio 


Kj  =  1  .0x10^  M"l 

^ j^   =  100  Kd 

Wg  =  100  Kd 

''AT  "    3.333  mg/ml 

Cgj  =  6  .667  mg/ml 

[A]  =  3.333x10"^  M 

[B]  =  6  .667x10"^  M 
[B]/[A]  =  2.0 

Frictional  ratios  (f/f^)  for  all  species  =  1.1 

All  self  and  cross  hydrodynam i c  constants  (k)  =  0.01 

Species  transport  coefficients  are  given  in  table  2. 

Rotor  speed  =  60,000  rpm. 

Initial  sharp  boundary  position  (r  )  =  6,0  cm. 

At  =  2155.51  sec 

a.  no  self  association 

b-  K    =  5.0x10*  M-1 
*2 
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Figure  61:   Effect  of  Self  Association  of  B  on  the  Boundary 
of  an  Uncooperative  AB,  System  at  a  2:1 
Constituent  Ratio 


Kj  =  l.Oilo'  M"l 
W^  =  100  Kd    ■ 
Wg  =  100  Kd 
''AT  "  3.333  mg/ml 
Cg^  =  6.667  mg/ml 

[A]  =  3 .333xl0~5  M 

[B]  =  6.667il0"5  M     ■ 
[B]/[A]  =2.0 

Frictional  ratios  (f/f  )  for  all  species  =  1.1 

All  self  and  cross  hydrodynamic  constants  (k)  =  0.01 

Species  transport  coefficients  are  given  in  table  2, 

Rotor  speed  =  60,000  rpm. 

Initial  sharp  boundary  position  (r  )  =  6,0  cm. 

At  =  2155  .51  sec 

a.  no  self  association 

b.  K    =  5.0x10^  M-1 

^2 


c.  K„ 


2.5X10''  M-1 
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3 . 7    EFFE CT  OF  NON-INTERACTING  CONTAMINANTS 

Occasionally  one  or  both  of  the  constituent  solute  samples 
used  to  prepare  a  mixture  for  the  ul  tr  a  cen  t  r  i  f  uge  cell  may- 
be contaminated  with  inactive  monomer  and/or  fixed  dimer. 
Inactive  (crippled)  monomer  may  be  generated  during  extrac- 
tion by  procedures  that  can  alter  the  conformation  or  chemi- 
cal specificity  of  the  protein's  binding  site(s).  Fixed  di- 
mer and  other  fixed  self  aggregates  usually  result  from 
formation  of  disulfide  bonds. 

Since  the  presence  of  non- int e r a c t i ng  species,  in  a  sam- 
ple consisting  of  given  total  amounts  of  A  and  B,  lowers  the 
concentration  of  reactive  constituents  in  the  mixture  (ie. 
the  effective  constituent  concentrations),  it  follows  that 
the  initial  distribution  of  reactive  constituents  among  the 
monomer  and  aggregate  species  would  be  different  for  prepa- 
rations with  different  proportions  of  inactive  species.  It 
should  not  be  surprising,  therefore,  that  the  presence  of 
either  fixed  dimer  or  crippled  monomer  have  noticeable  ef- 
fects on  the  shape  of  a  system's  gradient  profile. 

Gradient  profiles  of  AB«  systems  with  various  amounts  of 
non- i nt e r a c t i ng  species  (monomer  or  dimer)  were  studied  us- 
ing the  distorted  grid  simulation  program  for  mixed  associ- 
ating systems  modified  to  include  a  third  channel  of  non-in- 
teracting contaminant.  The  modified  program  accounts  for 
cross  hydrodynamic  dependences  between  the  contaminant  and 
the  reactive  constituents  as  well  as  self  hydrodynamic  de- 
pendence for  the  contaminant. 
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3,7.1    £i.^^  d  P^m^  r 

Simulated  profiles  for  preparations  of  an  uncooperative  AB. 
system  with  different  amounts  of  fixed  dimer  in  the  constit- 
uent  A  component  and  the  constituent  B  component  of  the  mix- 
ture are  shown  in  figures  62  and  63. 

With  the  transfer  of  part  of  the  total  constituent  A  con- 
centration to  fixed  A^(fig,  62),  the  trail ing  boundary, 
which  consists  of  free  A,  becomes  smaller  as  excess  reactive 
A  is  depleted,  ^' 

Transferring  part  of  constituent  B  to  fixed  dimer(fig. 
63),  on  the  other  hand,  increases  the  trailing  boundary's 
size  -  probably  by  freeing  more  of  constituent  A. 


The  height  of  the  leading  boundary  is  also  changed  in  op- 
posite directions  by  the  presence  of  fixed  A_  and  fixed  B« . 
As  the  fraction  of  fixed  A  is  increased  the  resulting  in- 
crease in  the  size  of  the  non-interacting  dimer  boundary, 
which  sediments  just  inside  or  below  the  trailing  edge  of 
the  leading  reactive  constituent  boundary,  increases  the 
height  of  the  total  leading  boundary.  When  the  fraction  of 
fixed  B„  is  increased  the  gain  in  the  leading  boundary 
height  due  to  the  dimer  boundary  is  outweighed  by  a  reduc- 
tion in  size  due  to  a  "transfer"  of  constituent  A  from  the 
reaction  boundary  to  the  trailing  boundary.  Depleting  reac- 
tive B  increases  excess  reactive  A,  and  the  net  effect  is  a 
slight  decrease  in  leading  boundary  height. 
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Figure    62:       Effect    of    Fixed    A^    on    the    Boundary    of    an 
Uncooperative    AB,     System. 
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Figure  63:   Effect  of  Fixed  B,  on  the  Boundary  of  an 
Uncooperative  AB   System. 
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Contamination  of  either  A  or  B  with  fixed  dimer  decreases 
the  average  sedimentation  rate  of  the  leading  boundary.  The 
individual  constituent  and  fixed  dimer  boundaries  show  that 
this  decrease  is  due  to  the  contribution  of  the  slower  dimer 
boundary  to  the  total  leading  boundary  rather  than  a  change 
in  the  sedimentation  rates  of  the  leading  reactive  constitu- 
ent boundar  i  e  s . 

The  unusual  appearance  of  a  negative  gradient  in  the 
non- i nt e r a ct i ng  dimer  profiles  of  these  systems  probably 
arises  because  of  the  large  A  and  B  gradients  in  that  area, 
which  would  result  in  a  sharp  decrease  in  the  effective  sed- 
imentation coefficient  of  the  dimer,  causing  a  p i  1  e- up  of 
dimer  behind  the  reaction  boundary  (ie.  a  Johns t on-Ogs t on 
effect ) . 


3.7.2    C  rlBEl  e  d  M  o  n^m£  r 

Simulated  profiles  of  preparations  of  an  uncooperative  AB« 
system  with  various  amounts  of  inactive  A  and  B  monomer  are 
shown  in  figures  64  and  65  With  constituent  A  initially  in 
excess  (C  ^=Cgnr.=  5.0  mg/ml),  as  the  fraction  of  crippled  A 
monomer  is  increased,  the  crippled  A  boundary,  which  sedi- 
ments at  the  same  rate  as  the  noncrippled  A  monomer  bounda- 
ry, becomes  progressively  larger,  while  the  reactive  A  mo- 
nomer boundary  becomes  progressively  smaller.  The  net 
result  is  no  effect  on  the  trailing  boundary  size.  Increas- 
ing the  fraction  of  crippled  B  monomer  (with  A  initially  in 
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excess)  increases  the  amount  of  excess  A.  This,  along  with 
the  increase  in  the  contribution,  to  the  trailing  boundary, 
of  the  crippled  B  monomer  boundary,  results  in  a  sharp  in- 
crease in  the  total  Size  of  the  trailing  boundary. 

The  presence  of  either  crippled  A  or  B  reduces  the  size 
and  average  sedimentation  rate  of  the  leading  boundary  - 
probably  by  shifting  the  equilibrium  of  non-crippled  compo- 
nents  toward  the  smaller  species  (AB,  A  and  B). 

Since  the  presence  of  crippled  species  affects  boundary 
sizes  and  positions,  failure  to  recognize  or  consider  their 
presence  could  lead  to  faulty  interpretations  of  the  ob- 
served schlieren  patterns  and  make  modelling  very  difficult. 
Fortunately  the  presence  of  fixed  dimer  can  usually  be  de- 
tected in  the  schlieren  patterns  of  the  individual  constitu- 
ent samples  and  the  fraction  of  dimer  can  be  determined  ei- 
ther directly  from  area  measurements  or  indirectly  by 
modelling  and  simulation.  When  the  contaminant  is  crippled 
or  inactive  monomer,  however,  its  presence  is  less  likely  to 
be  recognized  beforehand,  as  it  has  no  effect  on  the 
schlieren  patterns  of  the  monomer  samples.  If  the  monomers 
are  undergoing  reversible  self  association  then  it  might  be 
possible  to  determine  the  presence  and  quantity  of  crippled 
monomer  from  modelling  studies  with  the  individual  constitu- 
ent preparations  at  different  concentrations.  If  the  amount 
of  crippled  monomer  cannot  be  assessed  prior  to  modeling  the 
mixed  system   then  erroneous  models  might  be  put   forth  for 
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Figure  64:   Effect  of  Crippled  A  Monomer  on  the  Boundary  of 
an  Uncooperative  AB.  System. 
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a.   no  crippled  A  monomer 
0.5  mg/ml  (of  total  C   )  present  as  crippled  A  monomer 
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Figure  65:   Effect  of  Crippled  B  Monomer  on  the  Boundary  of 
an  Uncooperative  AB   System. 
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the  underlying  reactive  system.  However,  since  the  concen- 
trations of  crippled  monomers  vary  linearly  with  the  total 
concentrations  and  the  distribution  of  active  constituents 
among  the  various  species  does  not,  it  seems  unlikely  that 
an  erroneous  reaction  model  (not  accounting  for  crippled  mo- 
nomer) would  give  adequate  matches  to  the  real  data  at  more 
than  a  few  select  constituent  ratios. 


Chapter  IV 
DISCUSSION 

In  general,  we  can  conclude  that,  under  many  circumstances, 
it  will  be  possible  to  distinguish  between  various  plausible 
models  for  an  AB-  system  by  comparing  simulated  velocity 
sedimentation  profiles  of  the  suspected  models  with  that  of 
the  real  system.  Specifically,  when  everything  about  a  sys- 
tem that  affects  its  behavior  during  velocity  sedimentation 
(ie.  boundary  shape),  except  for  the  values  of  the  associ- 
ation constants,  are  known,  then  it  should  be  possible  to 
discover  the  values  of  the  association  constants  by  trial 
and  error  modelling  and  simulating. 

It  has  also  been  shown  that,  in  some  instances,  it  will 
be  possible  to  differentiate  between  AB  systems  with  dif- 
ferent overall  stoichiometrics  (AB,,  AB,  and  AB , , ) .  Specif- 
ically, if  the  (AB  )  system  is  uncooperative,  and  the  stoi- 
chiometry  and  intrinsic  association  constant  (^j  )  are  the 
only  unknowns,  then  it  should  be  possible  to  identify  the 
correct  st oi chi ome try  and  K   simultaneously.    In  addition. 
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The  important  variables  that  contribute  to  the  boundary 
shape  for  a  rapidly  relaxing  AB,  system  are:  1.  the  asso- 
ciation constants,  K  and  Kj ,  2.  the  molecular  weights 
and  partial  specific  volumes  of  the  monomers,  3.  the 
frictional  ratios  of  each  species,  4,  the  hydrodynamic 
dependencies  of  the  sedimentation  coefficients  and  5.  the 
concentrations  of  inactive  species  in  the  preparation 
used  in  the  ul  tr  a  cen  t  r  i  f  uge  run. 
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it  has  been  shown  that  cooperative  AB.  and  AB^  systems  have 
boundaries  that  are  generally  different  from  one  another,  as 
well  as  from  the  boundaries  of  un  cooperative  AB«  and  AB^ 
systems.  Actually,  strictly  speaking,  we  have  shown  only 
that  these  distinctions  can  be  made  when  the  the  systems 
have  a  particular  set  of  properties  (W.=W|j  =  100Kd.,  all 
f  /  f  =1.10,  all  hydrodynamic  constants.  It.. ,=0,01,  and  no  in- 
active contaminants  or  self  association.),  and  it  is  not  at 
all  unlikely  that  systems  with  a  different  set  of  properties 
would  be  less  distinguishable.  Lower  molecular  weights  and 
larger  frictional  ratios  would  be  particularly  likely  to 
make  modelling  more  difficult  since  both  would  tend  to  ob- 
scure the  distinctive  bimodal  features  of  a  migrating  bound- 
ary . 

Of  course,  as  a  general  rule,  the  more  that  is  already 
known  about  the  system  before  attempting  to  model  velocity 
sedimentation  the  more  likely  it  is  that  such  modeling  stud- 
ies  will  resolve  the  residual  ambiguity  and  identify  or  con- 
f i rm  a  single  model  as  the  correct  one. 

In  particular,  it  is  likely  to  be  necessary  to  carry  out 
thorough  modelling  studies  with  the  individual  constituents 
before  attempting  to  obtain  an  unambiguous  result  for  mixed 
sy  s  t  em  s  . 

Usually  there  will  be  some  uncertainty  about  several  of 
the  variables  other  than  the  association  constants,  that  in- 
fluence  sedimentation   behavior,   the   most   troublesome   of 
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which,  for  our  purposes,  would  probably  be  the  aggregate 
frictional  ratios,  the  concentrations  of  inactive  species, 
particularly  monomer,  and  the  cross  hydrodynamic  dependen- 
cies. When  these  things  are  not  known  accurately  it  may  be 
necessary  to  settle,  in  the  end,  for  a  range  of  possible 
models  for  the  system,  since  these  variables  affect  the  same 
structural  features  of  the  gradient  profile  that  reflect  the 
values  of  the  association  constants  (the  heights,  widths, 
and  migration  rates  of  the  boundaries).  It  is  therefore  im- 
portant that  every  effort  be  made  to  define  plausible  ranges 
that  are  as  narrow  as  possible  for  all  of  the  variables,  in- 
cluding the  association  constants,  before  attempting  to  mod- 
el the  velocity  sedimentation  behavior  of  the  mixed  system. 
Thorough  velocity  sedimentation-modelling  studies  of  the  in- 
dividual components  should  be  completed  first.  The  trans- 
port properties  of  the  monomers,  self-association  constants, 
and  the  concentrations  of  inactive  species  will  usually  be 
easier  to  obtain  from  modeling  studies  with  the  monomers 
than  they  will  from  studies  with  the  mixed  system. 

In  practical  situations,  additional  information  is  usual- 
ly available  from  other  techniques,  particularly  from  such 
thermodynamic  equilibrium  techniques  as  sedimentation  equi- 
librium, light  scattering  and  osmometry.  A  correct  physical 
description  of  a  real  system  must  be  consistent  with  the  re- 
sults of  both  transport  and  equilibrium  experiments.  Nei- 
ther category  of  data  alone  will  usually  suffice  to  define 
any  but  the  simplest  real  system. 
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The  present  work  has  shown,  however,  that  transport 
experiments  can  radically  narrow  the  range  of  acceptable 
models  for  a  real  system.  Fitting  a  real  sedimentation  ve- 
locity profile  precisely  in  every  respect  -  the  number,  ve- 
locity, height  and  width  of  subsidiary  boundaries  -  is  quite 
difficult  to  do.  Fitting  is  particularly  difficult  when  an 
acceptable  model  must  be  consistent  with  the  behavior  of  the 
individual  constituents  and  must  describe  the  mixed  system 
at  several  constituent  ratios.  That  difficulty  confirms  the 
power  of  the  simulation  technique  to  distinguish  correct 
from  incorrect  models. 
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The  shape  of  the  migrating  boundary  produced  by  a  revers- 
ibly  associating  system  of  macromolecules  during  velocity 
sedimentation  in  the  ultracentrifuge  is  generally  thought  to 
be  a  fairly  distinctive  indicator  of  the  reaction  scheme 
(ie.  stoichiometrics  and  association  constants)  of  the  sys- 
tem. Computer  simulation  of  velocity  sedimentation  may  be 
useful  in  distinguishing  between  correct  and  incorrect  reac- 
tion models  when  other  techniques  (eg.  sedimentation  equi- 
librium, light  scattering,  and  osmotic  pressure)  give  ambig- 
uous re  s  ul  ts  . 

In  the  present  work,  velocity  sedimentation  in  the  ultra- 
centrifuge  was  simulated  for  a  variety  of  AB  mixed-associ- 
ating systems  using  the  distorted-grid  model  of  Cox,  The 
predicted  sedimenting  boundary  profiles  of  these  systems 
were  compared  in  order  to  evaluate  their  ability  to  differ- 
entiate among  dissimilar  systems.  We  have  found  that  1) 
among  AB^  systems  boundary  shape  may  be  an  unambiguous  indi- 
cator of  the  values  of  the  association  constants  K  and  K^ , 
and  that  2)  among  uncooperative  AB  systems  boundary  shape 
may  be  an  unambiguous  indicator  of  the  s t o i ch i ome try  and  as- 
sociation constant  (K  )  of  the  system.  In  general,  our  re- 
sults indicate  that  modelling  studies  of  velocity  sedimenta- 
tion data  can  be  a  quite  effective  approach  for 
distinguishing  between  various  plausible  models  for  a  real 
interacting  system. 


